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ABSTRACT 

A better understanding of the turbine flow field is needed in order for 
turbine designers to increase the efficiency and to improve the performance of 
turbines. This has been hampered by the fact the present knowledge of the axial 
flow turbine flow field is not adequate. Several of the major areas that need more 
investigation include turbine secondary flow, turbine nozzle and rotor wakes and 
rotor-stator interaction. In view of this need, the three-dimensional steady flow 
field in a turbine nozzle has been investigated experimentally with an emphasis on 
the nozzle secondary flow and the nozzle wakes. The two-dimensional unsteady 
flow field in a turbine rotor has also been investigated experimentally at midspan 
with an emphasis on the interaction of the nozzle wake with the rotor flow field. 
The nozzle wake decay upstream and in the rotor passage and the rotor wake decay 
characteristics were also measured and analyzed. 

The nozzle flow was measured at several axial locations. These 
measurements were carried out using a five hole probe, a two-component laser 
doppler velocimeter (LDV) and a single sensor hot wire. The static pressure on the 
nozzle endwall and vane surfaces were also measured. Based on these 
measurements, the nozzle secondary flow is found to be weak at midchord. At the 
nozzle exit, two secondary flow loss cores (casing and hub passage vortices) 
dominate the flow field. The casing passage vortex is larger in area than the hub 
passage vortex, while the hub passage vortex has a higher maximum loss than the 
casing passage vortex, which is contrary to the results of linear cascades (where the 
casing and hub passage vortices are identical.) Radial inward flow was observed 
over the whole passage, which was more pronounced in the wake. The nozzle wake 
decay was compared to the decay of other turbomachinery blade wakes. The nozzle 
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wake is shown to decay much faster than a compressor cascade wake, an annular 
turbine cascade wake or a turbine nozzle wake with a large nozzle-rotor spacing. 
This is due to the presence of a rotor in close vicinity. The radial variation of 
nozzle wake total velocity defects shows larger velocity defects in the hub and 
casing secondary flow regions. These large defects results from the interaction of 
the passage vortices and the wake which cause deeper and wider wakes. 

The steady and the unsteady flow field in the rotor was also measured at 
midspan using a two-component LDV. Measurements were acquired at 37 axial 
locations from just upstream of the rotor to one chord downstream of the rotor. The 
propagation of the nozzle wake through the rotor passage was captured by this 
method. As the nozzle wake enters the rotor passage it becomes bowed because the 
convection velocity at midpitch is higher than that near the rotor leading edge. As 
the nozzle wake travels further through the passage, it becomes distorted with the 
region of the nozzle wake near the rotor suction surface moving faster than the 
region near the pressure surface. The region of the nozzle wake near the rotor 
suction side continues to move faster than the region near the rotor pressure side, 
until the nozzle wake has turned more than 30 degrees from its orientation at the 
rotor inlet. It is now spread out along the rotor pressure surface from midchord to 
the trailing edge. This in contrast to other measurements of the nozzle wake in the 
turbine rotor passage, where at the rotor exit, the nozzle wake still spans the rotor 
passage in the pitchwise direction, extending from the pressure to suction surface. 
The reason for the difference is that the turbine rotor employed has a much larger 
ratio of the suction surface velocity to pressure surface velocity than the other 
researchers' turbines, which causes the region of the nozzle wake near the suction 
surface to travel much more rapidly than the region of the nozzle wake near the 
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rotor pressure surface. Thus the nozzle wake rotates about the region of the nozzle 
wake near the pressure surface, ending up spread out along the pressure surface. 

As the nozzle wake travels through the rotor flow field, the magnitude of 
the velocity defect grows until close to midchord, after which it decreases. This is 
contrary to the common theory that the nozzle wake would continuously decay as it 
travels through the rotor passage. Also, the nozzle wake total relative unsteadiness 
does not decrease steadily as the nozzle wake travels through the rotor passage. 
Instead, it increases and decreases as the nozzle wake travels through the rotor 
passage. 

High values of relative total unresolved unsteadiness were observed near the 
rotor leading edge, which result from an increase in the production of unresolved 
unsteadiness due to large mean flow gradients near the leading edge. High values 
of relative total unresolved unsteadiness were also observed near the rotor pressure 
surface. This increase in unresolved unsteadiness is due to the interaction of the 
nozzle wake with the pressure surface boundary layer, along with the concave 
curvature effects, which destabilize the flow. 

The rotor wake decay characteristics were also analyzed. Correlations are 
presented which match the decay of the various wake properties. The rotor wake 
velocity defect decays rapidly in the trailing edge region, becoming less rapid in 
the near and the far wake regions. The rotor wake velocity defect decays according 
to the inverse of the square route of streamwise distance downstream of the rotor. 
The decay of the maximum unresolved unsteadiness and maximum unresolved 
velocity cross correlation is very rapid in the trailing edge region and this trend 
slows in the far wake region. While the interaction of the nozzle wake with the 
rotor wake does not influence the decay rate of the various wake properties, it does 
change the magnitude of the properties. 



TABLE OF CONTENTS 


ABSTRACT iv 


LIST OF FIGURES xi 


LIST OF TABLES xxi 


NOMENCLATURE xxii 


ACKNOWLEDGEMENTS xxvii 


CHAPTER 1: INTRODUCTION 1 

1.1 Turbine Flow Field 2 

1.1.1 Secondary Flow 3 

1.1.2 Turbine Blade Wakes 1 3 

1.1.3 Rotor-Stator Interaction 15 

1.2 Objective 18 

1 .3 Method of Investigation 19 

CHAPTER 2: EXPERIMENTAL FACILTTTY AND MEASUREMENT 
TECHNIQUE 21 

2.1 Facility Description 21 

2.2 Design Features and Overall Performance Parameters 24 

2.3 Modifications to the AFTRF 30 

2.4 Static Pressure Instrumentation 31 

2.5 Five Hole Probe 32 

2.6 Hot Wire 32 

2.7 Laser Doppler Velocimeter 32 

2.8 LDV Data Analysis 36 

2.9 LDV Measurement Errors 43 



Vll 

CHAPTER 3: RADIAL DISTRIBUTION OF FLOW PROPERTIES AT 
DESIGN AND OFF DESIGN CONDITIONS 45 

3.1 Nozzle Inlet Flow Field 45 

3.2 Overall Performance 47 

3.3 Rotor Exit Measurements 51 

3.4 Off-Design Measurements 58 

CHAPTER 4: NOZZLE PASSAGE FLOW FIELD 63 

4. 1 Vane Surface and Endwall Static Pressure Distribution 63 

4.2. Flow Field Near Midchord (X/C = 0.56) 77 

4.3 Flow Field Near Trailing Edge (X/C = 0.935) 83 

4.3.1 Stagnation and Static Pressures 84 

4.3.2 Velocity and Flow Angles 86 

4.3.3 Secondary Flow Vectors 90 

4.3.4 Streamwise Vorticity 96 

CHAPTER 5: NOZZLE EXIT FLOW FIELD 101 

5.1 Total Pressure Loss and Static Pressure Drop 102 

5.2 Velocity and Flow Angles 112 

5.3 Secondary Flow Velocity 120 

5.4 Streamwise Vorticity 125 

5.5 Nozzle Wake Characteristics 134 

5.5.1 Nozzle Wake Static Pressure 134 

5.5.2 Total Velocity 136 

5.5.3 Nozzle Wake Decay Characteristics 139 

5.6 Mass-Averaged Properties 147 

CHAPTER 6: STEADY AND UNSTEADY FLOW FIELD AT ROTOR 
MIDSPAN 156 

6. 1 Measurement Procedure 1 57 

6.2 Cycle Average Properties 162 

6.2.1 Relative Total Velocity and Total Unresolved 

Unsteadiness 164 

6.2.2 Axial Velocity and Relative Flow Angle 1 74 

6.2.3 Unresolved and Periodic Unsteady Velocity 

Correlations 174 

6.3 Nozzle Wake Propagation through the Rotor (Rotor Time 

Resolved Flow Field) 185 



V1U 


6.3.1 Relative Total Unresolved Unsteadiness 186 

6.3.2 Relative and Absolute Total Velocity 194 

6.3.3 Relative Flow Angle 198 

6.3.4 Unresolved Velocity Cross Correlations 204 

6.3.5 Axial and Tangential Components of Unresolved 

Velocity Correlations 210 

6.3.6 Periodic Velocity Correlations 220 

6.3.7 Unsteady Velocity Vectors 227 

6.4 Mass-Averaged Properties 235 

6.5 Nozzle Relative Flow Field (Time- Averaged Properties) 253 

CHAPTER 7: ROTOR FLOW FIELD WAKE CHARACTERISTICS 259 

7.1 Nozzle Wake Characteristics in Rotor Passage 259 

7.2 Rotor Wake Characteristics 268 

7.2.1 Cycle- Averaged Properties 271 

7.2.2 Cycle-Averaged Wake Decay Properties 278 

7.2.3 Rotor Wake Profiles at Individual Nozzle/Rotor 

Locations 297 

CHAPTER 8: SUMMARY AND CONCLUSIONS 312 

8.1. Nozzle Flow Field 3 1 3 

8.2 Rotor Flow Field 315 

8.3. Rotor Wake Properties and Characteristics 319 

CHAPTER 9: RECOMMENDATIONS FOR FUTURE RESEARCH 322 


APPENDIX A: ERROR ANALYSIS 324 


APPENDIX B: FIVE-HOLE PROBE AND STATIC PRESSURE 
UNCERTAINTY ANALYSIS 328 

B.l Sources of Error 328 

B.2 Uncertainty Analysis 330 


APPENDIX C: HOT WIRE UNCERTAINTY ANALYSIS 


337 



ix 


C.l Probe Calibration 337 

C.2 Measurement Errors 339 

APPENDIX D: LDV ERRORS 342 

REFERENCES 346 



X 


LIST OF FIGURES 

Figure 1.1. Classical Secondary Flow Model of Hawthorne (1955) 4 

Figure 1.2. Endwall Flow Models by Klein (1966) and Langston (1980) 6 

Figure 1.3. Comer Vortex 8 

Figure 1 .4. Cascade Endwall Flow Structure (from Sharma and Butler, 

1987) 10 

Figure 1.5. Endwall Three-Dimensional Separation and Reattachment 

Lines (from Langson et al., 1977) 11 

Figure 1.6. Nozzle Wake Development and Propagation 14 

Figure 2. 1 . Schematic and Flow Path of AFTRF 22 

Figure 2.2. Nozzle Vane Profile 28 

Figure 2.3. Isometric View of the Rotor 29 

Figure 2.4. Decomposition of Instantaneous Velocity 40 

Figure 2.5. Velocity Decomposition for Phase-Locked Averaged Rotor 

Blade 44 

Figure 3.1. Axisymmetric Flow Field 1 Chord Upstream of Nozzle 46 

Figure 3.2. Turbulence Intensity 1 Chord Upstream of Nozzle 48 


Figure 3.3. Total and Static Pressure Profiles 1 Chord Upstream of Nozzle 


.49 



Figure 3.4. Velocity Profiles at 1 Chord Upstream of Nozzle 50 

Figure 3.5. Performance Curve 52 

Figure 3.6. Total and Static Pressure Profiles 2 Chords Downstream of 

Rotor 53 

Figure 3.7. Velocity Profiles 2 Chords Downstream of Rotor 54 

Figure 3.8. Radial and Tangential Flow Angles 2 Chords Downstream of 

Rotor 55 

Figure 3.9. Off-Design Velocity Profiles at 2 Chords Downstream of Rotor 59 

Figure 3.10. Off-Design Radial and Tangential Flow Angle Profiels at 2 

Chords Downstream of Rotor 60 

Figure 3.1 1. Velocity Triangles 61 

Figure 4. 1 . Nozzle Flow Field Measurement Location 64 

Figure 4.2. Experimental Static Pressure Distribution on Nozzle Surface 65 

Figure 4.3. Nozzle Vane Profiles at 5 Radial Locations 67 

Figure 4.4. Comparison of Predicted and Measured Static Pressure Profiles 68 

Figure 4.5 Nozzle Pressure Coefficient based on Exit Total Velocity 71 

Figure 4.6a. Nozzle Casing Static Pressure Coefficient 72 

Figure 4.6b. Nozzle Hub Static Pressure Coefficient 73 



XII 

Figure 4.7a. Nozzle Pressure Surface Static Pressure Coefficient 74 

Figure 4.7b. Nozzle Suction Surface Static Pressure Coefficient 75 

Figure 4.8. Comparison of LDV and FHP Total Velocity Data at X/C = 0.56 78 

Figure 4.9. Total Velocity at X/C = 0.56 79 

Figure 4.10. Yaw Angle at X/C = 0.56 80 

Figure 4.11. Total Turbulence Intensity (Tu) at X/C = 0.56 81 

Figure 4.12. Total Pressure Loss Coefficient at X/C = 0.935 85 

Figure 4.13. Static Pressure Coefficient at X/C = 0.935 87 

Figure 4.14. Total Velocity at X/C = 0.935 88 

Figure 4.15. Yaw Angle at X/C = 0.935 89 

Figure 4.16. Radial Flow Angle at X/C = 0.935 91 

Figure 4.17. Secondary Flow Velocity Vectors at XJC = 0.935 92 

Figure 4.18a. Secondary Flow Velocity Vectors at X/C = 0.935 94 

Figure 4.18b. Secondary Flow Velocity Vectors at X/C = 0.935 95 

Figure 4.19. Streamwise Vorticity at X/C = 0.935 98 

Figure 5.1. Total Pressure Coefficient at 1.025 103 

Figure 5.2. Total Pressure Coefficient at 1.09 104 



Xlll 

Figure 5.3 Maximum Total Pressure Loss For Hub and Casing Secondary 

Flow Regions 105 

Figure 5.4. Maximum Width of Hub and Casing Secondary Flow Regions 106 

Figure 5.5. Circumferentially Mass- Averaged Pressures at X/C = 1.025 108 

Figure 5.6. Static Pressure Drop Coefficient at X/C = 1.025 1 10 

Figure 5.7. Static Pressure Drop Coefficient at X/C = 1 .09 Ill 

Figure 5.8. Circumferentially Mass- Averaged Velocities at X/C = 1.025 113 

Figure 5.9. Yaw Angle at X/C = 1.025 1 14 

Figure 5.10. Yaw Angle at X/C = 1.09 115 

Figure 5.11. Circumferentially Mass- Averaged Angles at X/C = 1.025 117 

Figure 5.12. Radial Flow Angle at X/C = 1.025 118 

Figure 5.13. Radial Flow Angles at X/C =1.09 119 

Figure 5.14. Secondary Flow Velocity Vectors at X/C = 1.025 121 

Figure 5.15. Secondary Flow Velocity Vectors at X/C = 1.09 122 

Figure 5.16. Secondary Vorticity at X/C = 1.025 126 

Figure 5.17. Secondary Vorticity at X/C = 1.09 127 

Figure 5.18. Pitchwise Averaged Distribution at X/C = 1.025 and 1.09 130 



XIV 


Figure 5.19. Nozzle Exit Secondary Flow Vorticity 132 

Figure 5.20. Static Pressure Variation at X/C = 1.025 135 

Figure 5.21. Total Velocity Profiles at X/C = 1.025 137 

Figure 5.22. Total Velocity Profiles at X/C = 1.025 138 

Figure 5.23. Decay of Total Velocity Defect with Streamwise Distance 

(Midspan) 140 

Figure 5.24. Total Velocity Profiles at Midspan 142 

Figure 5.25. Variation of Wake Edge Velocity with Downstream Distance 143 

Figure 5.26. Decay of Radial Velocity 145 

Figure 5.27. Radial Variation of Total Velocity Defect 146 

Figure 5.28. Variation of Static Pressure and Semi-Wake Width with 

Streamwise Distance 148 

Figure 5.29. Mass- Averaged Total Pressure Loss 149 

Figure 5.30 Mass-Averaged Static Pressure Coefficient 151 

Figure 5.31. Predicted and Measured Mass- Averaged Total Velocity 152 

Figure 5.32. Predicted and Measured Mass-Averaged Yaw Angles 153 

Figure 5.33. Area-Averaged Total Pressure Loss 155 

Figure 6.1. Rotor LDV Measurement Locations 158 



XV 


Figure 6.2. Rotor Inlet Velocity Triangles 163 

Figure 6.3. Cycle-Averaged Relative Total Velocity 165 

Figure 6.4. Cycle-Averaged Relative Total Velocity 167 

Figure 6.5. Cycle-Averaged Relative Total Velocity Vectors 168 

Figure 6.6. Cycle- Averaged Relative Total Unresolved Unsteadiness 169 

Figure 6.7. Cycle-Averaged Relative Total Unresolved Unsteadiness 171 

Figure 6.8. Cycle- Averaged Total Unresolved Unsteadiness 172 

Figure 6.9. Cycle-Averaged Axial Velocity 175 

Figure 6.10. Cycle-Averaged Relative Flow Angle 176 

Figure 6.1 1. Cycle Averaged Axial Velocity 177 

Figure 6.12. Cycle-Averaged Relative Flow Angles 178 

Figure 6.13. Cycle-Averaged Unresolved Velocity Correlation 179 

Figure 6.14. Cycle-Averaged Unresolved Velocity Correlation 180 

Figure 6.15. Cycle- Averaged Periodic Velocity Correlation 182 

Figure 6.16. Cycle-Averaged Unresolved Velocity Correlation 183 

Figure 6.17. Cycle-Averaged Unresolved Velocity Correlation 184 



xvi 

Figure 6.18. Relative Total Unresolved Unsteadiness at the Six 

Nozzle/Rotor Positions 1 87 

Figure 6.19. Blade to Blade Profiles of Total Unresolved Unsteadiness 193 

Figure 6.20. Relative Total Velocity at the Six Nozzle/Rotor Locations 195 

Figure 6.21. Blade to Blade Profiles of Relative Total Velocity 199 

Figure 6.22. Relative Flow Angle at the Six Nozzle/Rotor Positions 200 

Figure 6.23. Relative Flow Angle Upstream of Rotor (Five Hole Probe 
Measurements) 203 

Figure 6.24. Blade to Blade Profiles of Relative Flow Angle 205 

Figure 6.25. Unresolved Velocity Correlation at the Six Nozzle/Rotor 

Positions 206 

Figure 6.26. Blade to Blade Profiles of Unresolved Velocity Correlation 21 1 

Figure 6.27. Unresolved Velocity Correlation at the Six Nozzle/Rotor 

Positions 212 

Figure 6.28. Unresolved Velocity Correlation at the Six Nozzle/Rotor 

Positions 215 

Figure 6.29. Periodic Velocity Correlation at the Six Nozzle/Rotor Positions 221 

Figure 6.30. Periodic Velocity Correlation at the Six Nozzle/Rotor Positions 224 

Figure 6.3 1. Periodic Velocity Correlation at the Six Nozzle/Rotor Positions 228 

Figure 6.32. Unsteady Velocity Vectors at the Six Nozzle/Rotor Positions 231 



xvu 


Figure 6.33a. Mass-Averaged, Cycle-Averaged Relative Total Velocity 236 

Figure 6.33b. Mass-Averaged Relative Total Velocity 237 

Figure 6.34. Mass-Averaged, Cycle-Averaged Absolute Total Velocity 239 

Figure 6.35a. Mass-Averaged, Cycle-Averaged Relative Flow Angle 240 

Figure 6.35b. Mass-Averaged Relative Flow Angle 242 

Figure 6.36a. Mass-Averaged, Cycle-Averaged Axial Velocity 243 

Figure 6.36b. Mass-Averaged Axial Velocity 244 

Figure 6.37. Mass-Averaged, Cycle-Averaged Absolute Flow Angle 246 

Figure 6.38a. Mass-Averaged, Cycle-Averaged Relative Total Unresolved 
Unsteadiness 248 

Figure 6.38b. Mass-Averaged Relative Total Unresolved Unsteadiness 249 

Figure 6.39a. Mass-Averaged, Cycle-Averaged Unresolved Velocity 

Correlation 251 

Figure 6.39b. Mass-Averaged Unresolved Velocity Correlation 252 

Figure 6.40. Time-Averaged Total Unresolved Unsteadiness 254 

Figure 6.41. Time-Averaged Absolute Total Velocity Vectors 255 

Figure 6.42. Time-Averaged Axial Velocity 256 


Figure 6.43. Time- Averaged Unresolved Velocity Correlation 


258 



xvm 


Figure 7.1. Nozzle Wake Total Velocity Profiles Upstream of Rotor 260 

Figure 7.2 Nozzle Wake Total Velocity Profiles in Rotor Passage 261 

Figure 7.3. Variation of Peak Nozzle Wake Velocity Defect in Rotor 

Passage 263 

Figure 7.4. Nozzle Wake Total Unresolved Unsteadiness in Rotor Passage 265 

Figure 7.5. Variation of Peak Nozzle Wake Total Unresolved Unsteadiness 

in Rotor Passage 267 

Figure 7.6. Nozzle Wake Unresolved Velocity Correlation in Rotor Passage 269 

Figure 7.7. Variation of Peak Nozzle Wake Unresolved Velocity 

Correlation in Rotor Passage 270 

Figure 7.8. Rotor Wake Cycle-Averaged Relative Velocity Profiles 272 

Figure 7.9. Rotor Wake Cycle-Averaged Total Unsteadiness Profiles 274 

Figure 7.10. Rotor Wake Cycle- Averaged Unresolved Streamwise-Normal 
Velocity Correlation Profiles 276 

Figure 7.11. Rotor Wake Cycle- Averaged Relative Flow Angle () Profiles 277 

Figure 7.12. Decay of Cycle- Averaged Rotor Wake Velocity Defect 279 

Figure 7.13. Variation of Cycle- Averaged Semi-Wake Width with 

Streamwise Distance 282 

Figure 7.14. Decay of Rotor Wake Cycle- Averaged Maximum Relative 

Total Unresolved Unsteadinesses 284 



xix 

Figure 7.15. Decay of Rotor Wake Cycle- Averaged Maximum Unresolved 


Velocity Correlation 287 

Figure 7.16. Decay of Rotor Wake Cycle- Averaged Maximum Unresolved 

and Periodic Velocity Correlations 289 

Figure 7.17. Decay of Rotor Wake Cycle- Averaged Maximum Unresolved 

and Periodic Velocity Correlation 290 

Figure 7.18. Decay of Rotor Wake Cycle- Averaged Maximum Unresolved 

and Periodic Velocity Correlation 291 

Figure 7.19. Variation of Rotor Wake Cycle- Averaged Momentum 

Thickness with Streamwise Distance 294 

Figure 7.20. Variation of Rotor Wake Cycle- Averaged Shape Factor with 
Streamwise Distance 295 

Figure 7.21. Rotor Wake Relative Velocity Profiles 298 

Figure 7.22. Rotor Wake Total Unresolved Unsteadiness 300 

Figure 7.23. Rotor Wake Unresolved Velocity Correlation 301 

Figure 7.24. Rotor Wake Relative Flow Angle 302 

Figure 7.25. Rotor Wake Relative Velocity Profiles at 304 

Figure 7.26. Rotor Wake Total Unresolved Unsteadiness at 305 

Figure 7.27. Decay of Rotor Wake Velocity Defect with Streamwise 

Distance for Each Nozzle/Rotor Location 307 



XX 

Figure 7.28. Variation of Rotor Wake Semi-Wake Width with Streamwise 
Distance for Each Nozzle/Rotor Location 308 

Figure 7.29. Decay of Maximum Relative Total Unresolved Unsteadiness 

with Streamwise Distance for Each Nozzle/Rotor Location 309 

Figure 7.30. Decay of Maximum Unresolved Velocity Correlation with 
Streamwise Distance for Each Nozzle/Rotor Location 310 



XXI 


LIST OF TABLES 

Table 2.1. Design Performance Parameters 25 

Table 2.2. The Design Features of the AFTRF 26 

Table2.3. Experimental Operating Conditions 27 

Table 6.1. Axial Measurement Locations in the Rotor 159 



XXII 


NOMENCLATURE 


B Bias error 

C Nozzle axial chord at midspan 

C L Nozzle axial chord 

Cp Nozzle static pressure coefficient (p — P Q| )/ 0.5p V 2 
Cp 2 Nozzle static pressure coefficient (p - P 0] )/o.5pV 2 2 
C p ^ Pitch angle calibration coefficient for five hole probe 
C p ^ Yaw angle calibration coefficient for five hole probe 

Static pressure calibration coefficient for five hole probe 
C Pm Total pressure calibration coefficient for five hole probe 
H Spanwise distance (R-R h )/(R, -R h ) 

H Shape factor 

j Individual measurement window location in the rotor passage 
k Kinetic energy 

k Nozzle/rotor position 

L Semi- wake width 

L s Maximum pitchwise width of the passage vortex from the wake center 
M Mach number 

rh Mass flow rate 

n Total number of measurements in each measurement window 

NB Number of rotor blades 

NRP Number of nozzle/rotor positions 

NW Number of measurement windows 

NWB Number of measurement windows per blade 

P Precision error 



XX1X1 


p 

PS 

^atm 

R 

Re 

S 

S 

S 

S r 

ss 

Tu 

Tu r 

Tu t 


Static pressure on nozzle surface 
Pressure surface 
Atmospheric pressure 
Radius 

Reynolds number 

Percentage pitchwise distance in nozzle passage (Percentage pitchwise 
distance from nozzle pressure surface to nozzle suction surface) 

Percentage pitchwise distance from the wake center at midspan (positive on 
pressure side, negative on suction side) 

Precision error 

Percentage pitchwise distance in rotor passage (Percentage pitchwise 
distance from rotor pressure surface to rotor suction surface) 

Suction surface 


Turbulence intensity, — — 


2 + v* 2 ) 


x 100% 


Relative total unresolved unsteadiness. 


^(tT 2 +v' 2 )/2 

W 


x 100% 


J{ 1T 2 + v ,2 )/2 

Total unresolved unsteadiness, — x 100% 


U Uncertainty 

U m Blade speed at midspan 

u' Fluctuating (unresolved) velocity in axial direction 
u Periodic velocity in axial direction 

V Absolute total velocity 

v' Fluctuating (unresolved) velocity in tangential direction 

v' Fluctuating (unresolved) velocity in axial direction 
v' Fluctuating (unresolved) velocity in tangential direction 



XXIV 


v Periodic velocity in tangential direction 

V c Defect in absolute velocity at the wake center normalized by U m 
(V r ) ref Radial velocity at midspan passage, midspan 
V sec Secondary Velocity 
W Relative total velocity 

W c Defect in relative velocity at the wake center normalized by U m 
X Axial distance from nozzle leading edge at midspan 

X L Axial distance from nozzle leading edge 

X r Axial distance from rotor leading edge at midspan 

Y Percentage pitchwise distance in the rotor wake with wake center equal to 
zero 

Z Axial distance from nozzle or rotor trailing edge normalized by nozzle or 
rotor axial chord, respectively 

a Absolute tangential flow angle or yaw angle (measured from axial 
direction) 

a p Primary flow angle 

a 0 Nozzle vane outlet angle 

P Relative tangential flow angle or yaw angle (measured from axial direction) 
P 0 Rotor blade outlet angle 

y Radial flow angle 

8’ Displacement thickness 

e Dissipation 

0 Momentum thickness 

C, Total pressure loss coefficient (P ol - P o )/0^pV x 2 , 

T| Efficiency 

p Density 



XXV 


'F Total pressure loss coefficient (p o1 - P 0j )/< 0.5 pU* 

'Floss Total pressure loss coefficient (P 0] - P 0 )/o.5 pU^ 

% Static pressure drop coefficient (p, - p)/05pU^ 

^ Static pressure difference in the wake (p^ - p^ )/ (p™, - 

to Vorticity 

<o s Streamwise vorticity (Equation 4.5, normalized by C/V x i) 
co Sb Streamwise vorticity at midspan 

®sec Secondary vorticity (Equation 5.1, normalized by C/V x j) 
Subscripts 

1 Inlet to nozzle 

2 Exit to nozzle 

3 Exit of rotor 

b blade 

h Hub 

i Instantaneous 

1 local 

m Midspan 

max Maximum 

min Minimum 

n Normal 

r Rotor 

ref Reference conditions 
rev Revolution 

s Streamwise 

s Stator 



sec Secondary flow 

o Total 

r,0,x Radial, tangential, axial directions 
t Tip 

te Trailing edge 


Superscripts 

— Circumferentially mass-averaged properties (Chapters 3, 4 and 5) 

— Ensemble-averaged properties (Chapters 2, 6 and 7) 

= Axea-averaged properties (Chapters 3, 4 and 5) 

= Time-averaged properties (Chapters 2, 6 and 7) 

r\ Circumferentially mass-averaged properties (Chapters 6 and 7) 
Cycle-averaged properties 
Fluctuating quantity 
Periodic quantity 
A Aperiodic quantity 

u Phase-locked averaged 

* Nondimensionalized variable 



ACKNOWLEDGEMENTS 


The research work reported here is supported by NASA Lewis Research Center through 
Grant NAG-3-555, with R. Boyle, C. Civinskas and J. Schwab as technical monitors. Comments 
and criticisms and support by these individuals and L. J. Bober, R. Gaugler and L. Povinelli of 
NASA Lewis Research Center is very much appreciated. The funding for the facility was provided 
by DoD through Instrumentation Grant (DAAL-03-86-G-0013), NASA Grant (NAG-3-555) and the 
Pennsylvania State University. The NASA personnel mentioned earlier and T. Doliglski and R. E. 
Singleton of Army Research Office are acknowledged for their assistance. The LDV System used 
in this investigation was acquired with the funds provided by NSF through Equipment Grant 
MSME-85-06626. 

The assistance by D. Ristic, S. Khalatov and J. Fetterolf in this experimental program is 
acknowledged. 




CHAPTER 1 
INTRODUCTION 


The major goals of the axial flow turbine designer are to increase the 
aerodynamic, thermal and mechanical performance of the turbine, and to lower 
manufacturing effort and cost. To increase the turbine performance, a better 
understanding of the flow field is needed. This has been hampered by the fact the 
present knowledge of the axial flow turbine flow field, especially the rotor flow 
field is not adequate. The flow field is complex being three-dimensional and 
unsteady, with the presence of laminar, transitional and turbulent regions near the 
blade surface. Some of the three-dimensional effects present include a three- 
dimensional boundary layer which can be laminar, transitional or turbulent, shock 
wave boundary layer interaction, radial density (or temperature) gradient, non- 
uniform entry flow and temperature field, steady and unsteady flow, and leakage 
and secondary flows. The three-dimensional viscous flows and turbulence effects 
are mainly caused by the three-dimensional boundary layers on blades and wakes, 
annulus wall and hub wall boundary layers, shock boundary layer interaction, and 
secondary flows in annulus wall and hub wall boundary layers. The presence of a 
horseshoe vortex at the leading edge combined with thick blades and high turning 
makes the flow field truly complex. There are many other basic problems related to 
turbines that remain unresolved, including rotor-stator interaction and its effect on 
the flow field, unsteady heat transfer and flow induced vibration. 

The other goal of the turbine designer is to lower manufacturing effort and 
cost, which means a reduction in the size, weight and number of blades. This 
requirement can only be met by reducing the blade row spacing and the number of 
stages, which involves increasing the aerodynamic stage loading. Reducing the 
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blade weight also results in decreasing the blade height. Thus the aspect ratio (the 
ratio of the blade height to the blade chord) of the blade decreases causing an 
increase in three dimensional viscous effects in the turbine flow field. This 
reduction in size and weight of the turbine also results in a decreasing gap between 
the nozzle and rotor, thus causing an increase in rotor-stator interaction effects. 

Thus a better understanding of the turbine flow field is essential for more efficient 
design of turbines. 

1.1 Turbine Flow Field 

The turbine flow field is very complex as is described above. Thus it is 
useful to describe the axial flow turbine flow field before proceeding. The flow 
field in the nozzle is described first. At the center of the blade passage, away from 
the blade surfaces and endwalls, no significant viscous effects are taking place. 

This region is basically inviscid and is called the inviscid core region. Along the 
blade surfaces, a three-dimensional boundary layer develops, which transforms into 
a three-dimensional wake downstream of the nozzle. The turning of the endwall 
boundary layers inside the nozzle induces three dimensional viscous flow effects 
called secondary flows. After the flow exits the nozzle it enters the rotor. The 
interaction of the nozzle wakes with the rotor and the interaction of the nozzle and 
rotor potential flow fields is called rotor-stator interaction. The rotor flow field is 
similar to the nozzle, except that it has the additional effects of the rotor rotation 
and the nozzle wake propagation through the rotor passage. 

Though the turbine inviscid flow field is understood well, the turbine 
secondary flow and wakes along with the turbine rotor-stator interaction are not 
completely understood as of yet. A complete description of these turbine flow field 
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phenomena and why they need to be investigated further is provided below. 

1.1.1 Secondary Flow 

Secondary flow in turbines is associated with the endwall flow and arises 
due to the turning of the endwall boundary layer by the blade row. This flow is 
highly three-dimensional and contains numerous vortices, including the passage 
vortex, the horseshoe vortex, the trailing filament and trailing shed vortices and the 
comer vortex. 

Figure 1.1 describes the classical secondary flow vortex system of 
Hawthorne (1955). It includes the passage, the trailing filament and trailing shed 
vortices. A comprehensive review and analysis of classical secondary flow theory 
is given by Horlock and Lakshminarayana (1973). The passage vortex is a result of 
the flow turning in the blade passage. The deflection of the flow in the inviscid 
region generates a cross-channel pressure gradient, from pressure to suction side, 
which is balanced by the centrifugal force in this region. But in the endwall 
boundary layer, the velocity of the flow decreases, and since the pressure in the 
boundary layer is the same as that of the inviscid core flow, the radius of curvature 
of the streamlines must decrease near the endwalls in order for the centrifugal and 
pressure force to be in equilibrium. Thus a cross flow is generated from the 
pressure to suction surface which rolls up into a vortex core to form a passage 
vortex. There are two passage vortices in each passage, one at each endwall, which 
rotate opposite to each other. The strength of the passage vortex is a function of the 
main flow turning, the thickness of the incoming end-wall boundary layer, and the 
velocity gradient. Presence of the passage vortex causes over-turning of the flow 
near the endwalls and undertuming of the flow away from the wall. 
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Figure 1.1. Classical Secondary Flow Model of Hawthorne (1955) 
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The trailing filament vortices and the trailing shed vorticity compose the 
vortex sheet at the trailing edge. The trailing filament vortices are generated due to 
the stretching of the inlet vortex filaments when passing through the cascade 
withdifferent velocities between the suction side and the pressure side. The trailing 
shed vorticity is caused by the spanwise change of blade circulation. The trailing 
shed vortices are shed all along the blade span, whereas the trailing filament 
vortices are confined to the endwall region. Together, these two vortices form a 
vortex sheet with the same sense of rotation, but opposite to that of the passage 
vortex. 

The horseshoe vortex is caused by the rolling up of the inlet boundary layer 
as it approaches the leading edge of the blade. It is called this because of its shape. 
The physical presence of the blade causes the flow approaching the blade to 
decelerate. A direct consequence of this flow deceleration is that the incoming 
endwall boundary layer has to negotiate an adverse pressure gradient. Thus, the 
inlet boundary layer (and the associated vorticity normal to the freestream) grows 
rapidly and separates to form a horseshoe vortex. One leg of the horseshoe vortex 
wraps itself around the pressure side of the blade, while the other leg wraps itself 
around the suction side. Although the horseshoe vortex is a well known 
phenomena, its existence in the turbine flow field has only been realized recently. 
According to Sieverding (1985b), Klein (1966) was the first person to identify the 
horseshoe vortex in a turbine passage, calling it a stagnation point vortex, (see top 
of Figure 1.2). But it was Langston et al. (1977) who first presented a detailed 
analysis of the development of the horseshoe vortex. The bottom of Figure 1 .2, 
from Langston (1980), shows the horseshoe vortex. 

The comer vortex is the vortex that rotates in the opposite sense to the 
passage vortex and is located right in the end wall/suction side comer. It is caused 
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KLEIN'S MODEL 



LANGSTON'S MODEL 


Figure 1.2. Endwall Flow Models by Klein (1966) and Langston (1980) 
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by the fact that the overturning of the flow near the endwall is so extreme that the 
flow interferes almost at a right angle with the blade suction surface. Sieverding 
(1985b) states that this interaction is similar to the one causing the boundary layer 
ahead of the leading edge to roll up into a horseshoe vortex. Because it is so small, 
it is hard to see, but its existence is often shown by a characteristic reduction of the 
overturning near the endwall, as shown in Figure 1.3a. Figure 1.3b shows endwall 
limiting streamlines indicating the generation of the comer vortex. The overturning 
at the endwall is such that the limiting streamlines interfere at almost a nght angle 
with the blade suction surface near the position of maximum surface curvature. 

The interaction of all these vortices causes the flow in the nozzle passage to 
be very complex and highly three-dimensional. Many researchers have conducted 
experiments to try to understand the nozzle flow field. Until the late 1970 s, most 
of this research was limited to surveys of flow conditions upstream and 
downstream of the turbine cascade, which is called the "control volume approach . 
Since no measurements were made within the blade passage, researchers had to 
infer from the inlet and exit measurements the processes occurring inside the 
cascade passages. Numerous loss correlations were developed from experiments 
using the control volume approach, but since these correlations were not based on 
the actual physics of the flow field, there were major discrepancies between the 
various correlations. There was no sound basis on which to decide what correlation 
to use for the performance of a particular turbine stage. This was strongly criticized 

by several researchers including Sieverding (1985a). 

Researchers thus realized that only by obtaining a better physical 
understanding of the evolution of secondary flow both inside and outside of the 
blade passage could accurate predictions of the losses occurring in the turbine be 
calculated. Several research programs were started in the mid 1970's to take 




Figure 1.3. Comer Vortex: 

(a) Reduction of Overturning Near Endwall Due to Comer Vortex (from 
Gregory-Smith and Graves, 1983) 

(b) Endwall Limiting Streamlines Indicating Generation of Comer Vortex 


(fromBelik, 1975) 
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detailed measurements of the endwall flow development. Most of this work was 
conducted in linear cascades. Sieverding (1985b) gives an excellent review of the 
research conducted in turbine cascades in order to provide an understanding of the 
basic aspects of secondary flows in turbine cascades. Sharma and Butler (1987) 
also present a detailed description of the evolution of secondary flow in a turbine 
cascade. From these and other researchers, a description of the flow field in a linear 
turbine cascade is presented as follows. 

The incoming boundary layer at both endwalls separates at the leading edge 
of the nozzle, forming a horseshoe vortex with two counter rotating legs, the 
pressure side leg and the suction side leg. This is shown in Figure 1.4. As the 
pressure side leg enters the vane passage, it is convected toward the suction surface 
of the vane, due to the pressure gradient between the pressure and suction surfaces. 
It meets the suction surface near the minimum static pressure point, lifts off the 
endwall and grows rapidly as it travels along the suction surface. As the pressure 
side leg moves across the passage, it merges with the passage vortex, while 
entraining low energy boundary layer fluid from the endwall and blade surfaces. 
Most of the fluid from the inlet endwall boundary layer is trapped in the passage 
and horseshoe vortex. Since all of the fluid particles from the inlet boundary layer 
have either become part of the passage and horseshoe vortices or been convected 
toward the suction side, a new boundary layer starts at the endwall downstream of 
the separation line defined by the passage vortex (see Figure 1.5). The suction leg 
of the horseshoe vortex follows the contour of the blade suction surface, until it 
reaches the point where the pressure side leg of the horseshoe vortex meets the 
suction side. At this point, the suction side leg of the horseshoe vortex leaves the 
endwall and orbits around the outer edge of the passage vortex, as it moves 
downstream. It is hard to detect this downstream of the vane, since it loses intensity 
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FLOW FROM 
INNER REGION 



Figure 1.4. Cascade Endwall Flow Structure (from Sharma and Butler, 1987) 
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in contact with the passage vortex. 

The separation lines on the endwall are shown in Figure 1.5 (from 
Sieverding (1985b), based on the endwall ink traces of Langston et al. (1977)). The 
flow field is divided into distinct regions through the three-dimensional separation 
lines S and the reattachment lines R (stagnation streamline) with the separation 
saddle point A at their intersection. The horseshoe and the passage vortices form 

behind the separation line S, with the pressure side leg of the horseshoe vortex and 
the passage vortex forming behind S p (starting from A) and the suction side leg of 

the horseshoe vortex behind S s . 

While this is a good description of the flow in a linear cascade, it is not an 
accurate one of the flow in a real turbine nozzle. Important turbine parameters such 
as a radial pressure gradient that can be the same order of magnitude as the 
pitchwise pressure gradient and radial variation in blade loading, which are present 
in an actual nozzle, do not exist in a linear cascade. Also, Boletis and Sieverding 
(1984) show that the presence of a downstream rotor causes local changes to the 
outlet flow of an annular cascade. Though there have been a few investigations 
carried out in annular cascades, most of them have been in configurations where 
there is no rotor behind the nozzle such as Yamamoto and Yanagi (1986), Binder 
and Romey (1982) and Sieverding, Van Hove and Boletis (1984). Even fewer 
measurements have been completed in nozzles with a rotor downstream and they 
were all carried out with a nozzle-rotor spacing that is much larger than is found in 
modem axial flow turbines, such as the experiments done by Richards and Johnson 
(1988) and by Joslyn and Dring (1990). Thus more work needs to be done in 
annular turbine nozzles with a downstream rotor and a realistic nozzle-rotor 


spacing. 
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1.1.2 Turbine Blade Wakes 

An understanding of turbine nozzle and rotor wakes are also important for 
the efficient design of axial flow turbomachinery used in land, space, naval and 
aircraft applications. A major cause of noise and vibration characteristics of 
turbomachinery is caused by wakes. Turbine wakes represent a source of loss in 
efficiency, since the mixing of the wakes with the freestream dissipates energy. 

The three-dimensional characteristics of the wake, the decay characteristics and the 
path that it follows is important in the design of the following blade rows. This 
information is essential for both the prediction of the aerodynamic and mechanical 
performance of a turbine and for building quieter turbomachines. An understanding 
of the wake development and its decay is also essential because of the role it plays 
in the rotor-stator interaction (see below). 

Both the nozzle and rotor wakes are three-dimensional. In the nozzle wake, 
shown in Figure 1.6, the radial component results from an imbalance between the 
pressure and centrifugal forces within the shear layer. The tangential velocity 
decreases as the surface of the nozzle is approached, and thus the centrifugal force 
decreases, while the radial pressure gradient remains constant. This imbalance 
between the centrifugal and the pressure forces sets up a radial inward flow, which 
continues even in the wake. On the other hand, the radial velocity is outward in a 
turbine rotor wake due to the centrifugal force. Although there has been extensive 
research done on compressor wakes (see Reynolds (1978), and Ravindranath 
(1979), for a complete literature survey on compressor wake research), not much 
work has been done on turbine wakes. Mee et al. (1990) measured the wake of a 
linear cascade at midspan at six axial locations. Sitaram and Govardhan (1986) 
reported the wake of a turbine rotor cascade blade at midspan at seven axial 
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Figure 1.6. Nozzle Wake Development and Propagation 
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locations. Goldman and Seasholtz (1992) measured the wake in an annular turbine 
cascade for three axial locations at midspan and Dring et al. (1987) reported wake 
data at one axial location in a 1 1/2 stage axial flow turbine with a nozzle-rotor 
spacing of 50% nozzle axial chord. Only one of these (Dring et al., 1987) had a 
realistic three-dimensional configuration with a rotor downstream. Dring et al. s 
rotor was located at 50% nozzle axial chord, which is much higher than usually 
employed in industry. Furthermore, the wake was measured at only one axial 
location. The only reported measurement of a wake from a rotating turbine rotor is 
also from Dring et al. (1987), and the rotor wake was measured at only one axial 
location, too. Thus, there is a need for a comprehensive survey of the wake at 
different axial and spanwise locations for a turbine nozzle with a realistic axial 
distance between the nozzle and the rotor, and also a need for a survey of turbine 
rotor wakes. 

113 Rotor-Stator Interaction 


The understanding of rotor-stator interaction is also important for the 
efficient design of axial flow turbines. The flow around the blades of a turbine is 
highly unsteady. This unsteadiness is caused by the aerodynamic interaction of the 
flow between the nozzle and rotor which is called rotor-stator interaction. Rotor- 
stator interaction can affect the aerodynamic, structural and thermal performance of 
a turbine. According to Dring et al. (1982), rotor-stator interaction can be divided 
into two parts. These are potential flow and wake interactions. The potential flow 
gradients extend both upstream and downstream of the blade and they decay 
exponentially with a length scale of the order of the pitch or chord of the blade 
row. If the axial gap between the blade rows is less than a chord (which it is in a 
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typical axial flow turbine), then the potential influence can cause unsteadiness in 
both upstream and downstream of the blades. However, the wake is convected 
downstream and has a far field rate of decay much lower than that of the potential 
flow. The wake will still be felt several chords downstream. But in most modem 
axial flow turbines, which have a rotor-stator spacing close to 20 percent of a blade 
chord, both the potential and wake effects will occur together. And in the future as 
gas turbine designers try to reduce weight, and thus the rotor-stator spacing gets 
smaller, these effects will become more prominent. 

Even though unsteady flow plays a major role in axial flow turbines, 
turbines are designed using three-dimensional steady flow calculation methods 
(Sharma et al., 1992). Empirical correlations are used to account for the effect of 
the unsteadiness. Because actual models of the loss generating mechanisms in 
unsteady flow turbomachinery do not exist, these correlations are based on results 
from stationary cascade data and do not represent the actual fluid mechanics in the 
flow field (Sharma et al., 1985). Thus these correlations must be multiplied by 
some factor to obtain a good estimate of the actual losses that occur in turbines. 
Although these correlations have worked well in the design of existing turbines, 
they do not represent the true physics of the flow field and are only useful in the 
areas from which they were obtained, namely design point predictions and turbines 
which are similar to existing designs (Hathaway, 1986). Thus, a more thorough 
knowledge of unsteady flow interactions is needed in order to increase both the 
design and off-design performance of existing turbines, and to design turbines 
which are considerably different than existing turbines. To obtain this knowledge, 
good time accurate data from inside the rotor is needed. This knowledge, in turn, 
can be used to model the unsteady flow mechanisms that are not currently in 
existing design codes. 
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Many researchers have investigated rotor-stator interaction. Greitzer (1985), 
Gallus (1987) and Sharma et. al. (1992) provide comprehensive surveys of 
unsteady flows in turbomachinery. Unfortunately, most of these investigations 
have been in isolated airfoils or axial flow compressors. And there should be 
significant differences in the unsteady turbine flow field as compared to the 
unsteady compressor or isolated airfoil flow field (Sharma et al., 1985). This is 
because the flow field is accelerating in a turbine where as in a compressor the 
flow field is decelerating and also because the flow turning is much larger in a 
turbine than in a compressor thus causing stronger secondary flow in a turbine than 
in a compressor. Also there is the added effect of the horseshoe vortex in a turbine 
that does not occur in a compressor. 

Only four groups have measured the unsteady flow field inside a turbine 
rotor; the groups at UTRC (Sharma et al., 1985), Cambridge (Hodson, 1984), 
Allison/Calspan (Rao, Delaney and Dunn, 1992) and at DFVLR (Binder et al. 
,1985, Binder, 1985 and Binder et al. ,1987) . The Cambridge group measured the 
flow field at midspan (nozzle-rotor spacing is 50% of nozzle axial chord), the 
Allison/Calspan group measured the blade surface pressures at midspan (nozzle- 
rotor spacing is 22.5% of nozzle axial chord), the UTRC group measured the blade 
surface pressures and the rotor exit flow field (nozzle-rotor spacing is 65% of 
nozzle axial chord), and the researchers at DFVLR used a L2F velocimeter to 
measure two components of velocity and turbulence intensity (axial and 
circumferential) at four axial planes inside of the rotor (nozzle-rotor spacing is 
61% of nozzle axial chord). Three out of the four research groups had a nozzle- 
rotor spacing of 50% of nozzle axial chord or greater which is not a realistic 
spacing for modem turbine designs. At this spacing the potential flow interactions 
will be very small and the nozzle wake will have decayed significantly by the time 
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it enters the rotor passage. Typically, modem turbines have a nozzle-rotor spacing 
of around 25% or less nozzle axial chord. The only group to take measurements 
with a realistic nozzle-rotor spacing (Allison/Calspan), measured only the blade 
surface pressures at midspan (and did not measure the flow field between the 
blades). Thus there is a need for flow field measurements in a turbine rotor with a 
realistic nozzle-rotor spacing that will include both potential flow and wake 
interactions. 

1.2 Objectives 

The objectives of this research are as follows: 

(1) To understand the three-dimensional flow field in an axial flow turbine 
nozzle with a focus on secondary flow. 

(2) To understand the nozzle wake characteristics and to show how the 
presence of a rotor in close vicinity affects the decay of a nozzle wake. 

(3) To understand the two dimensional unsteady rotor flow field at midspan. 

(4) To understand the effects of rotor- stator interaction on the flow field in 
an axial flow turbine. 

(5) To capture the nozzle wake trajectory through the rotor passage. 

(6) To determine the influence of the rotor on the decay of the nozzle wake 
in the rotor passage. 

(7) To understand the rotor wake characteristics and to show the decay and 
profile characteristics. 

(8) To understand how the interaction of the nozzle wake with the rotor 
wake affects the rotor wake decay characteristics 

(9) To generate accurate data for CFD code validation. 



The ultimate goal of this research is that it will be used to build more efficient 
turbines. 


1.3 Method of Investigation 

The axial turbine flow field has been investigated experimentally. The 
nozzle flow field was measured with a two-dimensional LDV, a five hole probe 
and a hot wire. The rotor flow field was measured using a two-dimensional LDV. 
These measurements were performed in a systematic manner. The turbine inlet 
asymmetry was determined first. Next, the performance of the turbine was 
measured in order to pick an operating point that was as close to the design point as 
possible. Then measurements were carried out inside the nozzle flow field. The 
nozzle surface and endwall static pressure were measured, and the nozzle flow field 
was investigated at two axial locations. The first location was at midchord using 
the LDV and the second just upstream of the trailing edge using a five hole probe. 
The LDV was used provide both the mean velocity and the turbulence intensity, 
while the five hole probe was used to obtain the total pressure loss in the passage as 
well as the three component of steady velocity. The data from the five hole probe 
can also be used to calculate the secondary flow vectors and the streamwise 
vorticity. Downstream of the nozzle, the flow field was measured also using a five 
hole probe. Using the measurements in the nozzle passage and nozzle exit, the 
evolution of the secondary flow vortices, and their importance to the flow field can 
be determined. The characteristics of the nozzle wake, including its decay rate and 
the growth of its width can also be determined. The nozzle exit flow field is also 
used to provide the flow field at the rotor inlet and to compare with the LDV 



20 


measurements upstream of the rotor. 

The rotor flow field at midspan was measured using a two-component 
LDV. The LDV was used to measure the steady and unsteady components of the 
axial and tangential velocity of the rotor flow field. The interaction of the nozzle 
flow field and the rotor flow field can be measured, including the propagation of 
the no zzl e wake through the rotor passage and its influence on the rotor flow field 
The rotor wake was also measured in detail at midspan in order to determine its 
characteristics, including the decay of the velocity defect and unsteadiness. 
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CHAPTER 2 

EXPERIMENTAL FACILITY AND MEASUREMENT TECHNIQUE 
2.1 Facility Description 

The Axial Flow Turbine Research Facility (AFTRF) of The Pennsylvania 
State University is an open circuit facility 91.4 cm in diameter and a hub to tip 
radius ratio of 0.73, with an advanced axial turbine blading configuration. The 
facility consists of a large bellmouth inlet, a turbulence generating grid section, 
followed by a test section with a nozzle vane row and a rotor as shown in Figure 
2.1. There are 23 nozzle guide vanes and 29 rotor blades followed by outlet guide 
vanes. The bellmouth inlet is housed in an enclosure (not shown) covered with 
wire mesh and a thin layer of rubber foam to filter the air prior to entry to the inlet. 

A variable through flow is provided by two auxiliary, adjustable pitch, axial 
flow fans and an aerodynamically designed throttle. The two fans in the series 
produce a pressure rise of 74.7 mm Hg (40 " of water) with a mass flow of 10 m^ 
per second under nominal operating conditions. The power generated by the 
experimental turbine rotor assembly is absorbed by an eddy-current brake which is 
capable of absorbing up to 60.6 kw (90 Hp). The speed of the rotor can be varied 
between 175 and 1695 rpm with the "dynamic" control system and can be held 
constant to ± 1 rpm, with normal fluctuations in line voltage. The eddy current 
brake is cooled by a closed loop chilled water cooling system. 

The rotor and nozzle vane passages are instrumented with fast response 
instrumentation to measure steady (time averaged), unsteady pressures and wall 
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shear stresses. A window for LDV measurements of the flow field upstream of the 
nozzle, in the nozzle passage, in the spacing between the rotor and the nozzle, in 
the rotor passage and downstream of the rotor passage is included. 

The facility is equipped with two traversing mechanisms. One of the probe 
traverse units is mounted directly behind the rotor disk and has provisions for the 
radial and circumferential traverses in the rotating frame. It is controlled by a 
stepping motor driven by a computer indexer at tangential increments of 0.019 
degrees/step to allow accurate measurements of the rotor passage and downstream 
flow field. The intra blade radial and tangential surveys are accomplished by a 
computer controlled indexer and stepping motor in 0.05 mm steps and 0.0234 
degree steps in the radial and circumferential directions respectively. The second 
probe traversing unit is mounted on the outside casing of the turbine, and it 
provides radial and circumferential traversing at any axial location from one chord 
upstream of the nozzle to two chords downstream of the rotor. A probe mounted on 
this traversing unit can be positioned to an accuracy of 0.05 mm in the radial 
direction, 0.03 degrees in the circumferential direction and 0.5 mm in the axial 
direction, allowing for very detailed measurements of the absolute flow field in the 
turbine. 

The rotating-to-stationary data transmission system, attached to the rotor 
shaft ahead of the nose cone, is an integral part of the facility. It consists of a 150 
ring mechanical (brush/coin type) slip ring unit, and a specialized ten-channel low 
noise/signal ratio mercury slip ring unit. A 32 channel electronic pressure scanner 
unit is located in the rotating drum downstream of the turbine rotor. The electrical 
signals carrying the pressure information is carried to the stationary frame through 
the slip ring assembly. The rotor frequency is accurately determined by using an 
optical shaft encoder. 
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A completely automated data processing system is built around a computer 
with a clock rate of 25 MHz. The system consists of a 32 bit computer with 8 Mb 
random access memory, a disk operating system, 5 10 Mb hard disk storage space, 
and a HP laser printer. All the data from both stationary and rotating 
instrumentation can be processed on line. 

2.2 Design Features and Overall Performance Parameters 

The overall performance parameters are given in Table 2.1. The principal 
aerodynamic and geometric design features of the nozzle vane and the rotor blades 
designed by the General Electric Company (Aircraft Engines Group) personnel are 
presented in Table 2.2. Provisions exist for changing the vane-blade spacing from 
20 percent to 50 percent of nozzle axial chord. All the experiments presented in 
this thesis were performed at 20% vane-blade spacing. The actual operation 
conditions of the AFTRF are given in Table 2.3. Isometric views of the nozzle 
vane and the rotor are shown in Figures 2.2 and 2.3, respectively. A more detailed 
explanation can be found in Lakshminarayana et al. (1992). The aerodynamic 
design, while not representing any specific current or future GE product, does 
embody modem turbine design philosophy. Stage loading flow coefficient, 
reaction, aspect ratios, and blade turning angles are all within the ranges of current 
design practice. State-of-the-art quasi-3D design methods were used to design the 
airfoil shapes. It is felt that the design is fully capable of meeting the intended 
research applications. 




Table 2.1 

Design Performance Parameters 


Total temperature at inlet (°K); T 0 

DESIGN 

289 

Total pressure at inlet (kpa); P 0 

101.36 

Mass flow rate (kg/s); W 

11.05 

Specific work 
output (kJ/kg); Ah/W 

5.49 

Flow function; W>/t/p 
( kgV°K m 3 / kN S) 

1.85 

Speed function (rpm/V° K ); N/Vt 

77.69 

Rotational speed; N 

1300 rpm 

Total pressure ratio; P 01 / P o3 

1.0778 

Total temperature ratio;T 03 /T ol 

0.981 

Pressure drop (mmHg); P o3 - P o] 

56.04 

Hub reaction; R 

0.181 

Pitchline reaction; R 

0.3820 

Pitchline loading coefficient; *¥ 

3.76 

Hub loading coefficient; ¥ 

5.27 

Stator Zweifel coefficient 

0.7247 

Rotor Zweifel coefficient 

0.9759 

Power (kw) 

60.6 

Stator efficiency; r\ s 

9.9421 

Rotor efficiency; r| r 

0.8815 

Total-to-total 
Isentropic efficiency; Tj-^ 

0.8930 
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Hub Tip Ratio 

Tip Radius 

Blade Height 

Rotor Blade (tip) 

rel. mach number 

number 

chord 

spacing 

turning angle 

maximum thickness 

tip clearance 

Nozzle Guide Vane (tip) 
number 
chord 
spacing 
turning angle 
maximum thickness 

Midspan Axial Chord 
nozzle 
rotor 

Auxiliary Fan 
pressure rise 
mass flow 
power 

Vane Reynolds Number 
based on inlet vel. 
based on exit vel. 

Blade Reynolds Number 
based on inlet vel. 
based on exit vel. 


Tri te 2 < 2 

The Design Features of AFTRF 

0.7269 
0.4582 m 
0.1229 m 


0.240 (max) 

29 

0.1287 m 
0.1028 m 

95.42 deg (tip), 125.69 deg (hub) 

22 mm 

1.27 mm (actual 0.97 mm ave., 1.04 mm max, 0.77 mm min) 


23 

0.1768 m 
0.1308 m 
70 deg. 
38.81 mm 


0.1123 m 
0.09294 m 


74.72 mm Hg 

10.39 per sec. (22,000.0 cfm) 
149.1 kw 


(3 -4)x 10 5 
(9- 10) x 10 5 


(2.5 - 4.5) x 10 5 
(5 - 7) x 10 5 



Table 2.3 


Experimental Operating Conditions 


Mass flow rate ( m ) 

10.53 kg/s 

Pressure Drop (p o1 /p o3 ) 

1.078 

Rotational Speed 

1300 rpm 

Nozzle/Rotor gap 
tip 

midspan 

0.1 8C 
0.226C 
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hub 


Figure 2.3. Isometric View of the Rotor 
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2.3 Modifications to the AFTRF 


The following is a summary of the modifications made to the AFTRF since 
its arrival in February 1990. 

(1) An intake wooden structure (10’ x 10' x 5’), covered with wire mesh and 
0.125" thick foam sheet , was designed and built. 

(2) An exhaust structure to reduce the noise and recirculate air in the room 
was designed and built. The exhaust structure uses very dense acoustic absorbing 
material in a wooden frame to reduce the auxiliary fan noise. 

(3) A control system for the eddy current brake has been designed and 
installed, which is used to monitor the rotor speed and shut down the facility if the 
rotor speed increases above a predetermined limit. 

(4) Slotted windows and a circumferential traverse have been designed and 
built for circumferential traversing. The aluminum window frame contains a 
removable Plexiglas window which has a circumferential slot. Three windows have 
been built with slots at five different axial locations. More windows can be built if 
needed. The slots contain two strips of rubber foam. As the probe moves between 
the two strips of rubber foam, the rubber forms a seal around the probe so that air 
will not leak through the slot. 

(5) A computerized data acquisition system has been added to the facility. It 
consists of a 96 channel scanivalve, an IBM compatible personal computer and a 
Metrabyte DAS-20 data acquisition board. Software has been written to acquire 
and reduce the static pressure and five hole data which uses this equipment setup. 

(6) The bearing that is located just downstream of the rotor was replaced by 
a bearing that can withstand higher axial loads. Also, a thermocouple was 
connected to the bearing and an Omega CN100 six channel monitor was obtained 
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so that the bearing temperature can be closely monitored. 

(7) The six cylindrical struts that hold the slip ring in place have been 
replaced with three aerodynamically shaped struts. This decreased the turbulence 
level at the three measuring stations located one chord upstream of the nozzle from 
4 percent to 1.5 percent, which is the turbulence level without the struts in place. 
These measuring stations are directly downstream of the struts. 

(8) An optical shaft encoder was connected to the turbine shaft. This 
encoder enables LDV and hot wire measurements to be made in the rotating frame. 

(9) Glass and Plexiglas windows have been designed and built for LDV 
measurements in both the nozzle and rotor. 

2.4 Static Pressure Instrumentation 

The Axial Flow Turbine Research Facility is equipped with a large number 
of static pressure holes (nearly 500) at carefully selected locations. A list of the 
locations is as follows: (a) Static holes for performance measurements on the outer 
casing and inner casing of the rig at 4 axial planes. A circumferential average at 
each plane will be used for performance calculations, (b) 43 static holes on the 
nozzle (hub) endwall in one passage, (c) 43 static holes on the annulus wall of the 
same nozzle passage, (d) 154 static pressure holes at several axial and chordwise 
locations on both surfaces of the nozzle blades (one passage), (e) 52 static holes on 
the hub surface of one rotor passage, (f) 154 static pressure holes at several radial 
and chordwise locations of rotor blade on both surfaces (one passage). 
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2.5 Five Hole Probe 


A five hole probe was used to measure the flow field inside of the nozzle. 
The five hole probe is a miniature one, with a probe head diameter of 1 .67 mm. 

The five hole probe was calibrated so that it could measure flows with a velocity 
direction oriented ±30 degrees to the probe tip in both the yaw and pitch directions. 
The estimated errors in the five hole probe measurements are given by Sitaram et 
al. (1981), and are as follows: wall and blade vicinity effects, Reynolds and Mach 
number effects, misalignment of probe, probe blockage effect and turbulence 
effects. A detailed description of the error analysis and the accuracy of the five hole 
probe is given in Appendix B. 

2.6 Hot Wire 


A single sensor hot wire was used with a Dantex 55M01 anemometer, to 
measure the nozzle inlet turbulence intensity. The length of the tungsten sensor was 

9.5 mm and its diameter was 5 microns. The probe calibrations were corrected for 
temperature variations using the method of Kristensen (1973). The hot wire 
measurement errors are discussed in detail in Appendix C. 

2.7 Laser Doppler Velocimeter System 

The LDV is a TSI two color, four beam, two-dimensional measuring 
system. It consists of a seven Watt Argon-Ion laser tuned to the 488 nm (blue) and 

514.5 (green) lines. These beams are then split and one beam of each color is 
passed through a Bragg cell where it is frequency shifted. The beams then pass 
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through a 3.75 X beam expander and a 190 mm diameter lens which focuses the 
beams to a probe volume with a length of 0.65 mm based on the 1/e^ intensity 
points. The half angle of the intersecting beams is 4.95 degrees. The green and blue 
beam pairs travel down the same optical path and are used to measure axial and 
tangential components of velocity, respectively. Atomized mineral oil or Rosco fog 
fluid is used to seed the flow. Scattered light from the seed particles passing 
through the probe volume is collected on-axis, i.e. the receiving optics for the blue- 
green light is located on the blue-green transmitting optics. On-axis collection was 
used instead of off-axis to improve access into the rotor blade. 

The LDV system is mounted on an optical table which is attached to a 
mechanical traverse. The traverse can move vertically, parallel and perpendicular 
to the AFTRF axis (the axial and radial directions, respectively) and also can be 
tilted. The three linear degrees of freedom plus tilt enables positioning the probe 
volume very accurately. 

The flow in the AFTRF in the nozzle is measured through a curved 
Plexiglas window, conforming to the inside of the casing and in the rotor through a 
flat glass window. Both the Plexiglas and the glass are 3.175 mm (1/8 inch) thick. 
The beams are aligned through the glass by placing a blade solid surface inside the 
AFTRF at the probe volume, disconnecting the photomultiplier and attaching the 
eyepiece to the receiving optics to see the probe volume. The photomultipliers are 
then attached and the Bragg cells turned on. If both probe volumes are aligned 
correctly a clear sine wave with the downmix frequency set by the frequency 
shifter will be visible on the oscilloscope. 

The nozzle, rotor and end walls were painted with a high temperature flat 
black paint to reduce reflections. Both regular flat black and fluorescent orange 
paint were tried, but the laser beams burnt off these paints. 
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The LDV system is free to acquire velocity measurements whenever a seed 
particle crosses the LDV probe volume. In the rotor flow field , this results in the 
random acquisition of many velocity measurements during every rotor revolution. 
Thus, an optical shaft encoder has been attached to the turbine shaft. This encoder 
divides one revolution of the rotor into 6000 counts and since there are 29 rotor 
blades, this comes to 206.9 counts per rotor passage. One rotor passage is then 
divided into 50 measurement windows, each which have 4.14 counts and is 1.72 
mm in length at midspan. Each velocity measurement is tagged with the angular 
position of the rotor by means of the optical encoder. 

Since a two dimensional LDV is used for measurements in the AFTRF, 
only the properties in the axial and tangential directions are measured. In order to 
reduce the reflections of the laser beams from the glass or Plexiglas window on the 
receiving optics, the LDV system is offset in the tangential directions by 7.6 
degrees. This causes the tangential velocity to be in error by a small amount (< 1%) 
which was considered acceptable. 

There are two methods that can be used to acquire the data inside the rotor. 
The first method is for data to be recorded only when a measurement window is 
enabled. The advantage of this method is that the data is for only one measurement 
location which makes post processing easier. The disadvantage is that since LDV 
data is acquired only at one measurement window at a time, the AFTRF will have 
to run a long time to complete the measurements in an entire rotor passage. 

The second method is to take the data at all the rotor locations at once, 
tagging each of them with the appropriate rotor location using the encoder. Then 
after the measurement is done, the data is sorted according to rotor shaft location. 
The advantage of this method is that the measurements at all the tangential 
locations in a blade passage will be acquired at once, thus shortening the time it 
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takes to complete the measurements in the rotor passage. The disadvantage is that 
more post processing of the data needs to be done (because the data has to sorted 
according to rotor shaft location). The second method will be used since speed in 
acquiring the data is more important than speed in post processing the data. 

To process the data, two counter-type signal processors (one for each 
channel) are interfaced with a IBM PC compatible computer. The signal from the 
photomultiplier tube is fed into the counter through the frequency shifter. In the 
counter, the signal is passed through a band pass filter to reduce noise (and the 
pedestal if not frequency shifted). The counter measures the time a particle takes to 
cross eight fringes with a resolution of one nanosecond. To validate the signal, the 
transit time for crossing five and eight fringes is compared and when the two agree 
within a selected percentage, the data is accepted as a valid signal. The signal is 
then digitized and sent to a buffer where it is held. If a valid burst from the other 
channel is received within a preselected time interval (typically 50 psec in the 
nozzle and 10 psec in the rotor) the measurement is then sent to the computer and 
stored. If not the measurement is rejected and the process starts again. 

Atomized mineral oil was used as seeding in the nozzle. TSI's six jet seeder 
is used to atomize the oil. The atomized oil has an aerodynamic mean diameter of 1 
pm (TSI, 1987). A 25.4 mm (one inch) diameter tube is connected to the atomizer 
and is positioned 0.9 m (three feet) upstream of the AFTRF bellmouth intake (1.5 
m (five feet) upstream of the nozzle). The seeding tube can be moved radially and 
circumferentially in order to optimize the data rate each time the LDV probe 
volume is moved. Atomized mineral oil was tried in the rotor also, but the window 
became covered with oil very quickly (within 1 minute) and it became impossible 
to take measurements. A Rosco fog machine was used to seed the flow in the rotor. 
The seed particles exiting the fog machine have an aerodynamic mean diameter 



36 


and standard deviation of 1.1 and 1.9 pm, respectively (Wiedner, 1988). The fog 
machine produces a high data rate and does not dirty the window so quickly. 

2.8 LDV Data Analysis 

The data analysis procedure is based on the methods of Strazisar (1985) and 
Suder et al. (1987). The LDV does not measure the instantaneous velocity V(t) 
continuously as a hot wire does; instead it measures discrete samples of V(t). For 
each measurement position, the data at all the rotor shaft locations will be acquired 
at once, tagging each of them with the appropriate rotor location using the encoder. 
Then after the measurement is completed, the data are sorted according to rotor 
shaft location. In the stationary frame of reference, the rotor shaft positions can 
represent values of time, thus the instantaneous velocity V(t) can be represented by 
the discretized instantaneous velocity V , which is the velocity at one particular 
rotor shaft position during a given rotor revolution. The instantaneous velocity is 
than decomposed in the following manner, 

V(t) = V; = V + V rev + V rcv + V' (2.1) 

where V is the time-average velocity, V rev is the revolution periodic velocity 
(temporally fluctuating component), V rev is the revolution aperiodic velocity 

(passage to passage average) and V' is the unresolved velocity. The revolution 
periodic velocity can further be subdivided as follows: 

V_ = v b + V 

rev D r 


( 2 . 2 ) 
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where V b is the blade periodic velocity and V b is the blade aperiodic velocity. The 

time averaged value describes a steady-state flow field which is the same in each 
blade of a blade row. The periodic component is caused by the relative motion 
between the blade rows. The aperiodic velocity results from different number of 
nozzle or rotor blades in each stage. In a single stage turbomachine, the aperiodic 
velocity should be zero (provided there are no differences in each blade passage 
due to manufacturing tolerances and blade installation). 

Since the flow field in turbomachinery is not stationary, ensemble averaging 
is needed to decouple the periodic unsteady velocity from the unresolved velocity. 
The ensemble-average of the velocity for each measurement window can be 
calculated according to the equation 


V = 



(2.3) 


where V is the ensemble-averaged velocity, V, is the instantaneous velocity 
measured at a particular measurement window during a given rotor revolution and 
n is the total number of measurements at that measurement window. The 
unresolved velocity for each measurement window can also be calculated as 

V'= V - V (2.4) 


and the corresponding variance as 


V' 2 = 


Y-v) 


(n-l) 


(2.5) 


The level of unresolved unsteadiness in each measurement window is 
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determined by the variance. 

The time-averaged velocity V is obtained by averaging all the ensemble- 
averaged velocities in each measurement window as follows, 

<2 - 6) 

where NW is the total number of measurement windows. The time-averaged 
velocity is a time average of all measurements acquired at a fixed point in space. 
To find the revolution periodic and aperiodic components, first a time-averaged 
velocity over one blade passage must be found which is 

v ‘ = nwb |r v (27) 

where NWB is the number of measurement windows in one blade passage. The 
revolution periodic velocity is then obtained from the equation 

V„=V-V„ (2.8) 


and the revolution aperiodic velocity is 

V rev = V b -V (2.9) 

The blade periodic velocity can then be found by averaging the revolution periodic 
component for each blade passage over all the blade passages as follows 

i NB 

v =— y v 
b NBt? rev 


( 2 . 10 ) 
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where NB is the total number of blades. The blade aperiodic component is then 
calculated as 

V b =V rcv -V b (2.11) 

This velocity decomposition is shown in Figure 2.4 (which is described in more 
detail below). 

A program was written to decompose the velocity into the various 
components. In order to verify that the program works correctly a sample input was 
tested. The sample input, S, is composed of 

S = 0. 5 + 0.5 sin 29x + 0. 2 sin x (2.12) 

where sin 29x has a period of one blade passage and sin x has a period of one 
revolution. S can be thought of as the ensemble-averaged velocity for each 
measurement window. The input and output to the program is presented in Figure 
2.4. The figure shows that the program successfully decomposed the test input, S, 
into its various components. 

Using this program to decompose the measured data in the rotor, the rms of 
the total aperiodic unsteadiness is found to be low. The total aperiodic unsteadiness 
is equal to 


V = V 

V T v rev 


(2.13) 


The rms of the total aperiodic unsteadiness is less than 1.0% of the pitchwise 
averaged relative velocity everywhere in the rotor flow Held. This demonstrates the 
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periodicity of the flow field between the individual rotor passages of this turbine. 
Thus all ensemble-averaged velocity components can be spatially phase-locked 
averaged as follows 


* RB 

IX, 


G, =J=L 
J N 


( 2 . 14 ) 


RB 


where G represents either the ensemble-averaged velocity or variance, the subscript 
i determines the particular rotor passage and j is the measurement window location 
relative to the i'th rotor passage and is the number of rotor blades. Since the 
passage to passage variation in velocity is very small, calculation of the phase- 
locked average velocity and ensemble averaged velocity are done simultaneously. 
After the data is acquired, they are placed into the appropriate measurement 
window. The window from individual blade passages are then overlaid so that an 
"average" rotor blade passage is obtained. The ensemble-averaged velocity and 
variance is then calculated for each bin in the "average" blade passage. This is done 
to increase the number of measurements in each window, since the bin count 
(number of measurements in each window) is low in certain bins in a few blade 
passages such as bins in the blade boundary layer and wake center just downstream 
of the rotor. This increases the number of measurements that the ensemble- 
averaged velocity and variance is derived from, thus reducing the error. Since all 
results from now on are phase-locked averaged, the superscript u will be dropped 
hereafter. 

The aperiodic velocity components are small and can be neglected and thus 
equation 2.1 reduces to 


V(t) = V + V + V' 


( 2 . 15 ) 
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where V is the time-averaged velocity (the subscript b is dropped since V b — V 
when V b = 0 ) and V is the periodic velocity ( V b = V rev when V b = 0 ). The 

decomposition for this "average" rotor blade passage is shown in Figure 2.5. From 
the above velocity decomposition, the axial, tangential and cross velocity 
correlations (both periodic and unresolved) can be computed, and are discussed in 
more detail in Chapter 6. 


2.9 LDV Measurement Errors 

The error in locating the probe volume is +/- 0.03 degrees in the 
circumferential direction, +/- 0.2 mm in the radial direction and +/- 0.25 mm in the 
axial direction. Errors in the measurement arise due to fluctuations in the flow field 
and random noise in the photomultiplier tube signal, statistical or velocity bias and 
angle bias, among others. The LDV measurement errors are discussed in detail in 
Appendix D. 
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Figure 2.5. Velocity Decomposition for Phase-Locked Averaged Rotor Blade 
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CHAPTER 3 

RADIAL DISTRIBUTION OF FLOW PROPERTIES AT DESIGN AND 

OFF DESIGN CONDITIONS 

The overall performance of the AFTRF is discussed in this chapter. The 
overall performance consists of determining the pressure drop across the turbine 
stage at various corrected mass flow rates. These quantities were determined from 
a radial survey of the flow properties at both upstream of the nozzle and 
downstream of the rotor. The performance was measured before other 
measurements in the AFTRF were taken so that the experimental operating point 
could be selected. The axisymmetry of the inlet flow field of the turbine is also 
discussed in this chapter. A five hole probe was used for both the performance 
measurements and to determine the inlet flow velocity profile axisymmetry while a 
single sensor hot wire was used to measure the inlet turbulence intensity. The 
performance measurements were obtained using five hole probe radial surveys 
taken one chord upstream of the nozzle and two chords downstream of the rotor. 
The inlet flow field axisymmetry was found using radial surveys by both the hot 
wire and five hole probe taken at three equally spaced circumferential locations one 
chord upstream of the nozzle. 

3.1 Nozzle Inlet Flow Field 

In order to check for the flow axisymmetry upstream of the nozzle and to 
determine the turbulence intensity profile upstream of the nozzle, a single sensor 
hot wire probe and a miniature five-hole probe were used to measure the radial 
distribution of axial turbulence intensity and mean velocities at three tangential 








47 


locations (120 degrees apart) at the inlet. The axial velocity, presented in Figure 
3.1, is found to be nearly identical at these locations, being within 0.68% of each 
other at midspan, thus confirming the axisymmetric nature of the flow. The axial 
turbulence intensity, shown in Figure 3.2, is nearly constant at around 0.75 to 1.5 
percent, except near the hub and casing. The static and stagnation pressure profiles 
are shown in Figure 3.3. The stagnation pressure was nearly uniform, with the 
exception of locations inside the wall boundary layer. The static pressure was also 
uniform radially. The static pressure measured by the five hole probe is compared 
to the static pressure measured by static pressure taps on the casing and hub 
endwalls. The agreement is good, which gives confidence in the accuracy of the 
static pressure measurements by the five hole probe. The axial, tangential and 
radial velocity profiles are shown in Figure 3.4. The freestream velocity was 29.0 
m/s with a Reynolds number of 3.5 x 10^ based on a basic nozzle vane chord at 
midspan. The wall boundary layers were turbulent, with a thickness of 
approximately 5% blade span at the hub and 10% blade span at the tip. (The 
displacement and momentum thicknesses for the hub wall boundary layer were 
0.80 mm and 0.66 mm, respectively, while for the casing boundary layer, the 
displacement and momentum thicknesses were 1.15 mm and 1.02 mm, 
respectively.) These figures also show that the radial and tangential velocities are 
almost negligible upstream of the nozzle. 

3.2 Overall Performance 

In order to determine the overall performance of the turbine, radial surveys 
of flow properties were measured one chord upstream of the nozzle and two chords 
downstream of the rotor. The performance measurements were earned out at five 
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Figure 3.3. Total and Static Pressure Profiles 1 Chord Upstream of Nozzle 
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Figure 3.4. Velocity Profiles at 1 Chord Upstream of Nozzle 
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different mass flows at the corrected design speed of (1300 RPM), using a 
miniature (1.67 mm) five hole probe. The choice was limited, due to series 
operation of two downstream axial flow fans. The five-hole probe data included 
large number of radial stations, with close intervals inside the hub and annulus wall 
boundary layers. Hence the mass flow calculated includes the blockage due to the 
annulus and the hub wall boundary layers. The pressure drop coefficient (loading 
coefficient) is mass-averaged based on the local axial velocity at exit given by 

t 

f W 27crdr 

¥ = ^ (3.1) 

j\27trdr 

h 

The mass-averaged pressure drop coefficient is shown plotted against 
corrected mass flow in Figure 3.5. As expected, the loading coefficient varies 
linearly with an increase in the mass flow and the measured loading coefficients 
near the design mass flow are very close to the design value. The facility will be 
operated at m = 10.5 kg/s during the research phase of the program. The mass- 
averaged pressure drop coefficient is closely matched at the design and the 
operating conditions. 

3.3 Rotor Exit Measurements 

The radial distribution of total and static pressures, total, axial, radial and 
tangential velocities, and rotor exit flow angle at design conditions were measured 
two chords downstream of the rotor blade row using a five hole probe. Figures 3.6 
through 3.8 show the experimental data measured close to the design point 
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Figure 3.6. Total and Static Pressure Profiles 2 Chords Downstream of Rotor 
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Figure 3.7. Velocity Profiles 2 Chords Downstream of Rotor 






Figure 3.8. Radial and Tangential Flow Angles 2 Chords Downstream of Rotor 
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m = 10.5 kg/s). These results are compared with the design distribution at m = 

1 1.0 kg/s, which is close to the operating condition. 

The radial distribution of stagnation pressure drop coefficient (loading 
coefficient, 'P) and the static pressure drop coefficient ( V F S ) are shown compared 
with values at the design point in Figure 3.6. The experimental values are higher 
design as also is the mass averaged loading coefficient presented in Figure 3.5. The 
cumulative effect of nozzle and rotor secondary flow near the hub wall results in 
higher loading coefficient in this region. The leakage flow and the resulting 
undertuming is responsible for lower loading coefficient in the tip region (H > 
0.95). The undertuming is prevalent in 5 % of the span from the tip. The secondary 
flow, which overturns the fluid, has a major effect in increasing the pressure drop, 
reaching a maximum value at H = 0.90. The minimum pressure drop occurs near 
the midspan, while the maximum pressure drop occurs near the hub and the 
endwall region. Hence, this flow is dominated by the hub and the annulus wall 
secondary flows and the tip leakage flow. The presence of loss core in the hub 
region, which moves radially outward results in lower loading coefficient (due to 
higher losses) away from the hub wall. This is evident from the loading coefficient 
distribution away from the hub and the annulus walls. The profile boundary layer 
losses are also substantial in these regions. 

The radial variation of axial, absolute tangential and radial velocities are 
plotted in Figure 3.7. The radial velocities are negligibly small at this far 
downstream location. The presence of a loss core results in lower axial velocity in 
the hub region, and flow acceleration in other radial location (e.g. H=0.3 to 0.7). 
The presence of secondary flow and wall boundary layer result in low velocities in 
the hub and tip regions. These loss core regions tend to increase flow blockage 
giving rise to acceleration of the flow near the midspan region. Similar 
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distributions have been reported by Joslyn and Dring (1990) and Boletis and 
Sieverding (1991). The tangential velocity distribution shows similar trend. Near 
the hub wall (from H = 0 to 0.05) the effect of the wall boundary layer is evident. 
Higher tangential velocities are observed in the region 0.04 < H < 0.4. The 
tangential velocities are close to design at most other locations, with the exception 
of the tip regions. The effect leakage flows are clearly in the region H > 0.95, 
where the tangential velocities are lower. 

These distributions are consistent with angle distributions (absolute flow), 
shown in Figure 3.8. The pitch angles are very small. The overturning of the flow 
near the hub wall and undertuming in regions away from the hub wall is evident 
from Figures 3.7 and 3.8. The overturning region is confined to 0.04 < H < 0.3 in 
the hub region. Substantial deviation from the design distribution of yaw angle is 
evident from Figure 3.8. The yaw angles show substantial overturning in regions 
0.04 < H < 0.3, undertuming in the middle third of the blade and undertuming in 
regions H > 0.8. Hence, four distinct regions can be recognized. 

(1) 0 < H < 0.04: This is the inner region of hub wall boundary layer, where the 
flow is nearly axial, with boundary layer type of distribution. 

(2) 0.04 < H < 0.3: This is the secondary flow and the loss core region, with 
large overturning (reaching peak value near H = 0.1). The axial velocities 
are low, tangential velocities are high and the yaw angles are large in this 
region. The radial velocities are small. 

(3) 0.3 > H < 0.7: This is the inviscid core region. The axial velocities are 
higher, the tangential velocities are closer to the design and the outlet angles 
are smaller in this region. 

(4) 0.7 < H < 0.95: This region seems to be influenced by secondary flow. The 
changes in outlet angle, axial velocity and tangential velocities are not as 
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large as those near the hub. But the fact that the flow is overturned with 
respect to the design indicates the presence of secondary flow in this region. 
(5) H > 0.95: This region is dominated by tip leakage effect, with appreciable 
undertuming, low tangential and axial velocities. 

These exit flow distributions are consistent with the measurements by 
Joslyn and Dring (1990) and Boletis and Sieverding (1991) measured at the exit of 
a turbine stage. 

3.4 Off-Design Measurements 

The AFTRF was run at four off-design conditions; one at a mass flow rate 
of 5.21 kg/s, the second at a mass flow rate of 5.38 kg/s, the third at a mass flow 
rate of 9.72 kg/s and the last at a mass flow rate of 10.35 kg/s. Figures 3.9 and 3.10 
present the radial distribution of the velocities and flow angles two chords 
downstream of the rotor at the lowest mass flow rate ( m = 5.21 kg/s). The 
spanwise distribution of the properties upstream of the nozzle are not presented 
since they are similar in shape to the design condition except that their magnitude 
is less. Figures 3.9 and 3.10 show that the flow properties at the off-design 
condition have a much smoother spanwise variation than those at the design 
conditions. This results from the decrease in strength of the secondary flow 
vortices which occurs due to two factors. The first reason is that at this off-design 
condition the incidence angle entering the rotor is negative, which can be seen in 
the velocity triangles presented in Figure 3.1 1 for both the off-design and design 
conditions at midspan. Since the flow enters the rotor at negative incidence, it 
experiences less turning as it passes through the rotor, which results in a decreased 
secondary flow strength. The second factor which causes a decreased secondary 
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flow vortex strength is the decrease in Reynolds number compared to that at 
design. Hodson and Dominy (1987) state that a decrease in Reynolds number 
increases the rate of decay of the secondary flow regions, thus causing the 
secondary flow vortices to be smaller at a given distance downstream of the rotor. 
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CHAPTER 4 

NOZZLE PASSAGE FLOW FIELD 

The nozzle passage flow field measurements were carried out in order to 
obtaining a better physical understanding of the evolution of secondary flow inside 
the blade passage. The flow field in the nozzle was measured using a variety of 
measuring techniques. The static pressure on the nozzle pressure and suction 
surface along with the hub and casing endwalls was measured using static pressure 
taps. The flow field near midchord (X/C = 0.56) was measured using a two 
component LDV in order to capture both the velocity and the turbulence intensity. 
A miniature five hole probe was used to measure the flow field near the nozzle exit 
(X/C = 0.935), so that the both the losses and secondary flow velocities and 
streamwise vorticity could be characterized. These measurement locations are 
shown in Figure 4.1. Measurements at X/C = -1.0, 1.025 and 1.09 were also taken 
and are described in Chapters 3 and 5, respectively. 

4.1 Vane Surface and Endwall Static Pressure Distribution 

Figure 4.2 shows the vane static pressure distribution at nine spanwise 
locations on the nozzle pressure and suction surfaces where X L is the local axial 
distance from the nozzle leading edge and C L is the local axial chord. The pressure 
is normalized by the inlet axial velocity at midspan one chord upstream of the 
nozzle leading edge. The variation in static pressure from the tip to midspan is 
much smaller than the variation from midspan to the hub. This is reflected in the 
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Figure 4.1. Nozzle Flow Field Measurement Locations 
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nozzle vane shape variation presented in Figure 4.3. The variation of the vane 
shape is small from the tip to midspan of the vane, while it is larger from the 
midspan to the hub. The pressure distribution away from the endwalls (H=0.10 to 
H=0.90) is rather well balanced and agrees with the design values (not shown). The 
pressure distribution near the tip (H = 0.90 to H = 1.0) shows the maximum 
departure from the design (not shown). The suction peak occurs earlier than design 
and this is probably caused by the presence of secondary flow and the associated 
vortex. The distribution on the pressure side is well behaved with very little 
departure from design. The flow near the pressure surface is mostly inviscid, due to 
insignificant boundary layer growth. Further interpretation of this data will be 
given later. 

Figure 4.4 compares the experimental static pressure distribution with the 
static pressure distributions computed using MacFarland's (1981) panel code and 
the quasi-three-dimensional inviscid code of Katsanis (1977). The panel code 
matches the experimental data well at H = 0.30 and H = 0.50. At these locations the 
panel code distribution matches the experimental data exactly on the pressure side 
but on the suction side the measured flow accelerates less rapidly than the 
computed flow up to about X L /C L = 0.70 and H = 0.5. Good correspondence 
between the computed and experimental distributions at H = 0.3 and 0.5 are to be 
expected since the turbine nozzle secondary flow does not have a major influence 
at midspan. At other spanwise locations where the secondary flow affects the flow, 
the comparison between the computed and the experimental pressure distribution is 
not so good. The agreement between the data and Katsanis' code is excellent at 
midspan and at H = 0.10 and 0.30 , while the agreement deteriorates near the outer 
wall. In general, the three-dimensional inviscid code due to Katsanis shows better 
agreement with the data than the panel code. Hence, some of the departure between 



nee 

























70 


the design and the measured data is caused by the inviscid effects. 

As explained later, the static pressure increases linearly from hub to tip. 
This increase in static pressure is balanced by a decrease in total velocity (with 
increasing radius). Later observations indicate that is indeed true. Hence, the static 
pressure was also normalized by the total velocity at the exit of the nozzle (Figure 
4.5). The pressure side distributions collapse to nearly constant values. The suction 
side shows a spanwise variation from X L /C L = 0.05 to 0.8. This variation is due to 

the spanwise differences in loading, the three-dimensional inviscid effects and the 

secondary flow effects that exist near the endwalls. 

Contour plots of the casing and hub wall static pressure coefficients based 
on inlet dynamic head are shown in Figures 4.6a and 4.6b, respectively. The cross- 
channel pressure gradients exist over most the endwall for both the hub and casing. 
Sjolander (1975) found that the shear stress trajectories in the endwall regions were 
essentially parallel to the pressure gradient, indicating that the energy of the fluid is 
so small that the inertial effects are unimportant. Thus, the endwall pressure 
gradient is indicative of the direction of the endwall shear stress and consequently 

the direction of the flow. Accordingly, the cross passage flow extends over most of 
the hub and casing walls. The minimum static pressure region on the hub (C p = -13) 

occurs at X/C = 0.8 (downstream of the throat), while the minimum pressure 
region on the casing (C p = -9.5) occurs at X/C = 0.5 (upstream of the throat). (The 

throat is located at X/C = 0.71 along the suction surface for the hub and at X/C = 
0.68 on the suction surface for the casing.) Graziani et al. (1980) show that this low 
pressure region is the location where the passage vortex lifts off the endwall 
surface and begins to grow rapidly when approaching the nozzle suction surface. 

Contour plots of static pressure on the nozzle pressure and suction surfaces 
are presented in Figures 4.7a and 4.7b, respectively. The flow field is two- 
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dimensional over the whole pressure surface which is demonstrated by the 
relatively constant pressure coefficient values in the spanwise direction. The static 
pressure varies slightly from the leading edge to midchord. Downstream of 
midchord, the pressure decreases rapidly and the pressure gradient is perpendicular 
to the tr ailin g edge over most of the span, with a slight deviation from this 
direction near the tip. Since the trailing edge is not radial but has a radial lean (see 
Figure 4.1), the flow moves in a slightly radial direction, but is still two- 
dimensional. This demonstrates that the flow field is two-dimensional over most of 
the pressure surface with radial inward flow over most of the span and an increase 
in radial inward flow at the nozzle trailing edge/casing comer. 

The flow field on the suction surface of the nozzle can be divided into three 
regions. The first region occurs from the leading edge to X/C = 0.40. The flow 
field is two-dimensional over the whole span in this region, with the flow 
accelerating rapidly from the leading edge to X/C = 0.40. The second region is a 
three-dimensional region and starts at the hub wall at X/C = 0.80. As shown in 
Figure 4.6b, the cross passage pressure gradient drives the flow and thus the 
passage vortex toward the suction surface. The passage vortex meets the suction 
surface at X/C = 0.80 where the low pressure region occurs (C p = -13) and then 

travels up the suction surface and toward the trailing edge. The path of the hub wall 
passage vortex (HPV) along the suction surface is shown in Figure 4.7b. The 
position of the passage vortex at the trailing edge was determined from the flow 
field measurements in the trailing edge region discussed later. 

The third region on the suction surface is also three-dimensional and 
begins at X/C = 0.5 at the casing wall. The flow behavior is similar to that of the 
hub wall region with the flow and the casing passage vortex traveling toward the 
suction surface as shown in Figure 4.6a. The passage vortex meets the suction 
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surface at the low pressure point (C p = -9.5), which occurs at X/C = 0.50 and then 
the casing passage vortex (CPV) sweeps down the suction surface toward the 
trailing edge, following the path shown in Figure 4.7b. The position of the passage 
vortex at the trailing edge was determined from the flow field measurements in the 
trailing edge region. One interesting feature of this flow field, is that contrary to 
what occurs in a linear cascade, the casing passage vortex intersects the suction 
surface farther upstream (X/C = 0.50 ) than the hub passage vortex intersects the 
suction surface (X/C = 0.80). Since the passage vortex grows much more rapidly 
after it intersects with the suction surface, leaves the endwall and starts moving up 
the suction surface, this indicates that the casing passage vortex would cover a 
larger area than the hub passage vortex. This is demonstrated later in this chapter. 

4.2 Flow Field Near Midchord (X/C = 0.56) 

The flow field near midchord (X/C = 0.56) was measured using a two 
component LDV in order to capture both the velocity and the turbulence intensity. 
These measurement were taken from hub to tip and blade to blade. Since the nozzle 
axial chord decreases from hub to tip and the LDV measurement plane is at a 

constant axial distance downstream of the leading edge of the nozzle, the 
measuring location as a percent chord decreases from hub to tip (X L /C L = 0.61 at H 
= 0.05 and X L /C L = 0.52 at H = 0.95). X/C = 0.56, which is the percentage chord at 
midspan, will be used to define this axial plane. The third component could not be 
measured at this location because of limited access. 

The LDV measurements at X/C = 0.56 are shown in Figures 4.8 through 
4. 1 1 . The blank areas in the contour plots represent regions where data could not 
be taken due to limited access in the nozzle. Figure 4.8 shows a comparison of total 
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velocity derived from LDV, five hole probe and the vane static pressure at 
midspan. The tangential distance is normalized by the distance between the 
pressure and suction sides where S = 0 is the pressure surface and S = 1 is the 
suction surface. Both the LDV and five hole probe data matches the velocity 
profile quite well. This data shows a weakly nonlinear variation from pressure to 
suction surface and is similar to the data presented by Goldman and Seasholtz 
(1992). This excellent comparison between various types of measurements 
indicates that error in five hole probe and LDV data is negligibly small. 

Total velocity contours are presented in Figure 4.9. The total velocity 
variation follows the usual inviscid trend and is in agreement with Yamamoto and 
Yanagi's (1986) data. In addition to the inviscid flow field, one can observe the 
development of the endwall boundary layers, which are rather thin. This is because, 
according to Sharma and Butler (1987), all the inlet boundary layer fluid particles 
have either become part of the passage vortex or been swept toward the suction 
side, which results in a new endwall boundary layer being formed in the nozzle 
passage downstream of the separation line defined by the passage vortex (see 
Figure 1 .5) . 

The yaw angle contours, presented in Figure 4.10, show a rather well 
behaved flow with a ±2 degree variation at most locations. (The design mean flow 
angle is 40 degrees at this axial location.) As expected in inviscid flow, the turning 
angle near the suction surface is higher than that near the pressure surface. There is 
a region of high undertuming in the comer formed by the casing and suction side. 
The angle in that location is 28.5 degrees. This is due to the passage vortex. The 
overturning region associated with this vortex could not be observed due to 
constraints imposed by the LDV measurement limitations. The secondary flow 
seems to be rather weak and it is in the early stages of development. This is also 
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confirmed by the turbulence intensity contours (Figure 4.1 1) which show that the 
maximum turbulence intensity in this region is about three-fold compared to one 
percent at the inlet. No evidence of secondary flow is seen at the hub. This also 
corresponds to the nozzle endwall and surface static pressure contours presented 
earlier (Figures 4.6a, 4.6b and 4.7b). These static pressure contours showed that 
while the casing passage vortex has traveled across the passage and met the suction 
surface at midchord (where it begins to grow rapidly), the hub passage vortex is 
still being convected across the passage and does not intersect with the suction 
surface until farther downstream. Thus the hub passage vortex should still be rather 
weak and hard to detect. In addition to the secondary flow phenomenon, one can 
observe the increased turbulence intensity both in the hub and in the tip endwall 
boundary layer regions in Figure 4.1 1. 

4.3 Flow Field Near Trailing Edge (X/C = 0.935) 

Near the trailing edge, the flow field was surveyed with a five hole probe. A 
five hole probe was employed to facilitate stagnation pressure loss evaluation as 
well as to measure the blade endwall flows more thoroughly. The axial plane 
position is parallel to the trailing edge and ranges from X L /C L = 0.93 at the hub to 

0.94 at the tip. The axial plane position will be defined by the axial distance at 
midspan, X/C=0.935. Data were taken at 23 spanwise locations, clustered near the 
endwalls, and 40 to 55 tangential locations clustered near the vane surfaces. 
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4.3.1 Stagnation and Static Pressures 

Figure 4.12 shows the stagnation pressure loss coefficient 'F uoss . The two 
high loss regions located in the comer formed by the suction surface and the 
annulus wall as well as the suction surface and the hub wall show that the 
secondary flow is fairly well developed here. These high loss regions are 
associated with the passage vortex caused by the secondary flow. High losses occur 
in the core region and the maximum loss coefficient observed is 1.33 near the tip 
and 1.86 near the hub. The tip loss core is located farther from the casing than the 
hub loss core is from the hub, which results from the radial inward velocities that 
exist in both the freestream and the vane boundary layer. This is also the cause for 
the loss core near the hub to have a higher loss coefficient than the loss core near 
the tip. This is consistent with the measurements made slightly downstream of the 
trailing edge (Chapter 5) at X/C = 1.025. A comparison of the losses in the two 
axial planes show that the loss contours are very similar in the tip and the hub 
regions, indicating no major redistribution of the losses as the flow progresses from 
X/C = 0.935 to X/C = 1.025. 

Another noticeable feature of this flow field is that the casing passage 
vortex covers more area than the hub passage vortex. This results from three 
reasons. The first is the radial inward velocities cause the hub passage vortex to be 
pushed against the hub, which causes it to be compressed, while the casing passage 
vortex moves away from the casing. The second reason can be found in the 
discussion of the nozzle end wall and surface static pressure contours presented 
earlier (Figures 4.6a, 4.6b and 4.7b). These static pressure contours showed that the 
casing passage vortex travels across the passage and meets the suction surface 
farther upstream than the hub passage vortex does. Since the passage vortex begins 
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to grow rapidly after it intersects with the suction surface, the casing passage 
vortex has a larger area due to the longer distance it has traveled between the 
endwall/suction surface intersection point and this measuring location. The third 
reason is that the casing inlet boundary layer is larger than the hub wall boundary 
layer, which also results in the larger extent of the casing passage vortex. 

The static pressure contours at this location are given in Figure 4.13. The 
most striking feature is the strong radial pressure gradient that exists over the 
whole span. In the comer formed by the annulus wall and the suction surface there 
is a low pressure region (static pressure coefficient is 2.2) indicating the presence 
of the secondary flow vortex. A similar patch can be seen near the hub surface as 
well. 

4.3.2 Velocity and Flow Angles 

Figure 4.14 shows the total velocity contours. The total velocity increases 
from the tip to the hub due to the existence of the radial pressure gradient discussed 
earlier. A low velocity region is observed near the comer formed by the suction 
surface and hub and casing walls. This is the region of intense mixing of secondary 
flow and wall and blade boundary layers. 

One of the more interesting features of this flow field is the yaw angle 
shown in Figure 4.15. Most of the total flow turning has taken place by this 
location and the maximum overturning of about eight degrees occurs at H = 0.85 
close to the suction side. The maximum undertuming region is about four degrees. 
This is observed slightly below (H = 0.80) and a little farther from the suction side 
than the overturning region. This over and undertuming of the flow is characteristic 










of a vortex. Otherwise, the flow is well behaved with an angle variation of four 
degrees across the entire passage at midspan. A similar secondary flow region is 
observed near the hub. The maximum overturning of about three degrees and a 
similar magnitude of undertuming close to the suction side can by seen. This seems 
to indicate a conventional secondary flow pattern both in the hub and tip region. In 
addition, the flow very close to the walls, across the entire passage is overturned 
slightly. 

The magnitude of radial velocity can be discerned from the radial flow 
angle distribution shown in Figure 4.16. Negative values of radial flow angle 
indicate that flow is toward the hub over most of the passage. There are two 
regions with much higher and much lower radial flow angles than the rest of the 
flow. They are located in the suction surface hub and casing wall comers. Near the 
casing the maximum radial flow angle is -68 degrees and the minimum is one 
degree, while near the hub the maximum radial flow angle is -37 degrees and the 
minimum is six degrees. These high and low radial flow angles are another 
indication of the presence of secondary flow vortices. 

4.3.3 Secondary Flow Vectors 

The secondary flow velocity vectors in the r-0 plane, shown in Figure 4.17, 
are derived from the measured data. The primary flow angle at X/C = 0.935 was 
taken from the measured flow angle at the midspan at the respective tangential 
position, since the flow at midspan was not affected by the secondary flow, while 
downstream of the nozzle the primary flow angle was set equal to the design exit 
angle (70.0 degrees). Using this procedure, the secondary velocities were derived 
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Figure 4.17. Secondary Flow Velocity Vectors at X/C = 0.935 
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from 


V = rV.+nV 


= (V) measurcd -(V) 


design 


(4.1) 


where V is the radial velocity and V n is the velocity normal to the streamwise 

direction (estimated as indicated above). Strong radial inward flow can be seen to 

occur over most of the passage at X/C = 0.935 (Figure 4. 17). Since there is radial 

inward flow over most of the passage, the secondary flow phenomena is obscured. 

In order to perceive the secondary flow more clearly, the radial velocity at midspan 
and at midpitch ((V r ) ) is subtracted from the velocity vectors. Figures 4.18a and 

b show the secondary flow vectors with (V r ) subtracted from them. At 

X/C=0.935 (Figure 4.18a), the passage vortex (clockwise) can be clearly seen near 
the suction surface/casing comer. Since the pressure side leg of the horseshoe 
vortex has the same direction of rotation as the passage vortex, this vortex could be 
a combination of the passage vortex and the pressure side leg of the horseshoe 
vortex. Other researchers such as Langston (1977) and Sieverding (1985b) believe 
that the pressure side leg of the horseshoe vortex merges with the passage vortex 
and the results here do not show otherwise. There is no evidence of the suction side 
leg of the horseshoe vortex near the casing (which would rotate opposite to the 
passage vortex). 

Radial inward flow exists in the vane boundary layers caused by the 
imbalance of the centrifugal force and the pressure gradient inside the boundary 
layer. The tangential velocity decreases as the surface of the blade is approached, 
and thus the centrifugal force decreases, while the radial pressure gradient remains 
constant. This imbalance between the centrifugal and the pressure forces sets up 
radial inward flow, which continues even in the wake. This radial inward flow is in 
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the same direction as those induced by the casing passage vortex close to the 
suction side and thus it augments the passage vortex. At the hub, on the other hand, 
there is evidence of a passage vortex, but it is much weaker than the vortex near the 
casing. The radial inward flow in the suction surface boundary layer near the hub is 
in the opposite direction to those caused by the passage vortex, and thus the radial 
inward flow seems to overshadow those induced by the passage vortex. This is 
shown more clearly in Figure 4.18b which is an enlargement of the secondary 
velocity vectors near the hub/suction side comer. The passage vortex in this figure 
is labeled P. Above the passage vortex, there seems to be a weak vortex rotating in 
the opposite direction. This could be the suction side leg of the horseshoe vortex 
and it is labeled H ss in this figure. 

4,3.4 Streamwise Vorticitv 

The components of vorticity in the axial, tangential and radial directions are 
given by: 


_dv r 

rd0 

*rV e ) 

rdr 

(4.2) 
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dx 
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The streamwise vorticity is then given by: 

co s =co x cosoc p +a) e sina p + a> r siny (4.5) 
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where a p is the primary flow angle and y is the radial flow angle. The streamwise 

vorticity is normalized by C/V x j . The primary flow angle is determined as 
indicated above. The streamwise vorticity derived from the data at X/C = 0.935 is 
shown in Figure 4.19. The tangential and radial components of the vorticity 
involve axial gradients and since at X/C = 0.935, measurements were not made at 
an axial plane inside the vane close enough to this location to determine the axial 
gradients, another method had to be used to estimate the axial gradients. Gregory- 
Smith et al. (1988) have developed a method, which was extended for 
compressible flow by Niehuis et al. (1990), to estimate the axial gradients. Using 
their method, the axial gradients can be estimated as follows. Consider the Navier- 
Stokes equation for a compressible fluid with constant viscosity, which is 


p ^ = -Vp + V (^p VV + (VV) T ) + V[X(V • V)] (4.6) 


Nondimensionalizing this equation gives 

P - = _V„- + J-{ V . [vv- + (VV- ) T ]) + V[V(V V- )]} (4.7) 


where the superscript * represents the nondimensional variables. Assuming high 
Reynolds number flow, the viscous terms become very small and can be neglected. 
Converting back to dimensional form yields 



F 67.80 

E 60.83 

D 53.87 

C 46.91 

B 39.94 

A 32.98 

9 26.02 

8 19.06 

7 12.09 

6 5.13 

5 -1.83 

4 -8.80 

3 -15.76 

2 -22.72 

1 -29.69 
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Substituting the definition of vorticity which is 


to = VxV 


(4.10) 


into equation (4.9) and assuming steady flow yields 


VV 2 


-Vxffl 
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Solving equation (4.1 1) for the components of vorticity in the radial and tangential 
directions yields 


to, = 


i dp i av 2 

— p.+. 


v, [prae 2r ae 


Vffl, 


<o e =~: 


‘idp iav 2 _ _ 
p¥ + i~aT +v ^ 


(4.12) 

(4-13) 


Thus, with the radial and tangential vorticity now known, the streamwise vorticity 
can now be solved. The streamwise vorticity is normalized with respect to inlet 
axial velocity and vane chord at midspan. Since to derive equations (4.12) and 
(4.13), the viscous terms of the Navier-Stokes equations were dropped, the 
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streamwise vorticity presented in Figure 4.19 should be looked at qualitatively, not 
quantitatively. 

The passage vortex near the casing at XJC = 0.935 is shown clearly in 
Figure 4.19. The large positive vorticity (magnitude of 30.0) in the suction 
surface/casing comer corresponds to the passage vortex. Near the hub wall/suction 
surface comer, a large negative vorticity region (magnitude of -14.1) corresponds 
to the passage vortex also. The large positive vorticity region above the passage 
vortex could possibly correspond to the suction side leg of the horseshoe vortex. 
The vorticity is zero in the freestream giving confidence in this method to calculate 
vorticity correctly. 

Six criteria are used to determine the existence of a vortex in a flow field. 
These are high total pressure loss, low static pressure, high vorticity, over and 
undertuming of the yaw angle, positive and negative radial flow angles and vortical 
motion in the secondary flow vectors. A vortex exists when these six flow 
phenomena occur at the same position in the flow field. The existence of the 
passage vortices at both the hub and the casing of X/C = 0.935 is confirmed since 
they both meet all six of the criteria. On the other hand, the existence of the suction 
side leg of the horseshoe vortex is more doubtful. There is no evidence of its 
existence in the casing/suction side comer. While in the hub/suction side comer 
there is a region of high positive vorticity which could correspond to the suction 
side leg of the horseshoe vortex, the other criteria are not met. Thus, its existence 
cannot be proven conclusively. 
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CHAPTERS 

NOZZLE EXIT FLOW FIELD 


The nozzle exit flow field measurements were carried out with the 
following objectives: 

(1) To obtain a better physical understanding of the evolution of secondary 
flow downstream of the blade passage. 

(2) To determine the characteristics of the nozzle wake, including its decay 
rate and the growth of its wake width. 

(3) Measurement of the flow field upstream of the rotor was needed for the 
analysis of the rotor flow field measurements presented in Chapters 6 and 7. 

The flow field downstream of the nozzle was measured using a miniature 
five-hole probe, with a probe head diameter of 1.67 mm. The exit flow is measured 
at 50 to 80 tangential locations across the one passage and at 21 radial locations at 
X/C = 1.025 and 1.09, and at midspan at X/C = 1.007, 1.01, 1.12 and 1.16. Figure 
4.1 shows the measurement planes at X/C = 1.025 and 1.09. For the following 
plots, the tangential distance is normalized by the nozzle pitch, where S = 0 is the 
wake center at midspan at X/C = 1.025. The positive values of S are on the 
pressure side of the wake and the negative values of S are on the suction side. The 
pitchwise extent of the plots is a little less than one blade pitch for both X/C = 

1.025 and X/C = 1.09. 
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5.1 Total Pressure Loss and Static Pressure Drop 

The contours of total pressure loss coefficient (T'loss) are shown for X/C = 
1.025 and X/C = 1.09 in Figures 5.1 and 5.2. There are two high loss regions, one 
near H = 0.10 and the other near H = 0.80 for both axial locations. These two high 
loss regions are caused by secondary flow vortices and their dissipation. The 
maximum loss regions occur close to the suction surface. Furthermore, the tip loss 
core is located further from the casing than the loss core near the hub wall, which 
results from the radial inward velocities that exist in the wake and in the free 
stream region. Comparing the loss at both the tip and hub between these two axial 
locations and X/C = 0 935, one can see that the peak loss coefficient increases 
slightly between X/C = 0.935 and X/C = 1.025 and then decreases dramatically at 
X/C = 1.09. This can be seen more clearly in Figure 5.3, which presents the peak 
total pressure loss coefficient (('Floss ) ) in both the hub and tip secondary flow 

regions as a function of axial distance. The dramatic decrease between X/C = 1 .025 
and X/C = 1.09 results from the interaction of the passage vortex and the wake, the 
interaction of the passage vortex and the freestream flow and also the decay of the 
passage vortex as it travels downstream of the nozzle. 

On the other hand, the size of the passage vortex remains basically the same 
between X/C = 0.935 and 1.09. Figure 5.4, which presents the maximum width of 
the hub and casing passage vortices as a function of axial distance, shows this more 
clearly. The maximum width of the passage vortex, L s , is defined as the maximum 
pitchwise width of the passage vortex from the wake center (for the data 
downstream of the nozzle) or the nozzle surface (for X/C = 0.935) to the edge of 
the passage vortex (where V F L0SS = 0.2). In contrast, while the width (pitchwise 

extent) of the secondary flow loss regions has not changed between the two axial 
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locations, the width of the loss region outside of the secondary flow loss region has 
increased. The circumferentially mass-averaged total pressure loss coefficient for 
X/C = 1.025 is given in Figure 5.5. The two loss peaks can be seen clearly in this 
figure. The peak near the hub is larger in magnitude but smaller in radial extent 
than the one near the casing. The reasons for the larger extent of the loss region at 
the casing was discussed in section 4.3.1 of this thesis. The larger magnitude of the 
hub loss region is a result of the radial inward flow inside the wake which causes 
the low momentum fluid in the wake to move toward the hub. Both Yamamoto and 
Yanagi (1986) and Hunter (1982) report similar distributions in their annular 
nozzle cascades. 

At midspan, the secondary flow effects are absent and the mass averaged 
loss coefficient is the profile loss coefficient for the vane, which is 0.05. This is 
close to the value of 0.06 predicted from the correlation by Kacker and Okapuu 
(1982) which is shown in Figure 5.5. The area averaged total pressure loss 
coefficient ('F) is calculated to be 0.077 at X/C= 1.025 and 0.0818 at X/C=1.09. 
The losses increase as flow travels downstream of the rotor due to mixing of the 
wake, the decay of the vortices and the growth of the endwall boundary layer. The 
secondary flow loss is the loss calculated by the subtracting profile loss from the 
total loss. The mass averaged secondary flow loss coefficient measured in this 
nozzle is equal to 0.0262. This is less than the value of secondary flow loss 
coefficient predicted from the correlation by Dunham (1970) which is 0.0503. 
There are two main reasons for this difference. The first is that Dunham's formula 
overestimates the secondary loss parameter for low-aspect ratio blades. Kacker and 
Okapuu (1982) have derived a correlation for low aspect ratio blades to include the 
dependence on aspect ratio. Substituting Kacker and Okapuu's aspect ratio 
dependence formula into Dunham’s secondary loss prediction yields a value of 
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Figure 5.5. Circumferentially Mass- Averaged Pressures at X/C = 1.025 
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0.0361 for the secondary loss coefficient which is much closer, but still higher than 
the measured value. The second reason for this difference is that the secondary 
flow loss correlations were developed from cascade loss data taken between 0.7 
and 1.0 chord downstream of the blade, and thus include the additional losses due 
to mixing, decay of the secondary vortices and boundary layer growth. Moore and 
Adhye (1985) measured the losses at three axial planes (from X/C=l.l to X/C=1.4) 
downstream of a turbine cascade. Extrapolating their data, the nozzle secondary 
flow losses are found to increase by 30% from trailing edge to 0.70 chord 
downstream. Applying this correction to Dunham's correlation results in a value of 
0.0253 which is within 2% of the measured secondary flow loss coefficient 
(0.0262). 

Contours of static pressure loss coefficient for X/C = 1 .025 and 1 .09 are 
given in Figures 5.6 and 5.7 and the mass average distribution of static pressure 
loss coefficient is given in Figure 5.5. The mass averaged experimental values are 
compared to the design values that were calculated using a two-dimensional 
throughflow analysis that solved the circumferentially-averaged equations of 
motion in the meridional plane using a streamline curvature technique (see 
Lakshminarayana et al., 1992). The presence of strong radial pressure gradient is 
evident in these figures. The static pressures are generally lower in secondary 
vortex regions and higher in the wake regions. (Higher static pressure coefficient 
indicates lower static pressure.) While the radial variation in static pressure is 
similar for both X/C = 1.025 and X/C = 1.09, the pitchwise variation is not. The 
pitchwise variation of static pressure is much smaller and smoother at X/C = 1 .09 
than at X/C = 1.025, due to the rapid decay of the wake static pressure. Detailed 
interpretation of the static pressure across the wake is discussed later in this 
chapter. 
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5.2 Velocity and Flow Angles 

The circumferentially mass-averaged velocities at X/C = 1 .025 are shown 
plotted in Figure 5.8. The total and tangential velocity plots show a decrease in 
velocity near H=0.1 and 0.8, which is due to the hub and tip wall secondary flows. 
The axial velocity profiles show a decrease near the endwalls, due to the presence 
of the wall boundary layers. The total, axial and tangential velocities are in 
agreement with the design values near midspan region, while deviating from the 
design values in the secondary flow regions near the endwalls. 

Contour plots of the yaw angles (a) at X/C = 1.025 and 1.09 are shown in 
Figure 5.9 and 5.10, respectively. The flow is undertumed on the pressure side of 
the wake and overturned on the suction side of the wake away from the endwalls. 
This is because the suction surface blade angle is higher than the design exit flow 
angle (70 degrees) and the pressure side blade angle is lower than the design exit 
flow angle at the trailing edge. This over and undertuming of the flow decreases as 
the flow travels downstream of the trailing edge from a value of 6 degrees of both 
over and undertuming at X/C = 1.025 (midspan) to 2 degrees of both over and 
undertuming at X/C = 1 .09 (midspan), which is a result of mixing of the wake. The 
maximum undertuming occurs near H = 0.1 and H = 0.8. This local undertuming 
of the flow is caused by secondary flow vortices. In the core region of these 
vortices at X/C = 1.025, undertuming as high as 25 degrees is observed at H = 0.80 
while near the hub the maximum undertuming is 15 degrees. There is an 
overturning of 1 1 degrees near the suction side of the wake at H=0.80 and of 8 
degrees near the hub wall also at this axial location. This phenomena is as expected 
from secondary flow theories and this overturning and undertuming of the flow 
located next to each other is another indication of the presence of vortices at these 
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locations. Between X/C = 1.025 and 1.09, there is a dramatic decrease in maximum 
over and undertuming in the secondary flow regions. The maximum over and 
undertuming in the casing secondary flow region decreases by 8 and 20 degrees, 
respectively, between X/C = 1.025 and 1.09, while at the hub the peak values of 
over and undertuming decrease by 3 and 8 degrees, respectively. This results from 
the secondary flow vortex decay which is caused by the interaction of the 
secondary flow vortices with the nozzle wake and freestream flow. The mass- 
averaged yaw angles, presented in Figure 5.11, shows the overturning near the 
casing and near the hub endwall due to secondary flow and vortices. The closest 
measurement at the casing endwall is 4.5% of the span from the casing, while the 
closest measurement to the hub is 3.0% from the hub. As one moves away from the 
endwalls the flow becomes undertumed achieving design flow near 20% of span 
from the hub and 30% of span from the casing. The larger undertumed region near 
the casing is due to the larger secondary flow region near the casing. This is 
consistent with predictions from secondary flow theories and measurement by 
others. The experimental values of yaw angle match the design values well near the 
midspan, while deviating from design in the secondary flow regions, which is 
expected, since the design code is two-dimensional. 

The contours of radial flow angles (7) at both downstream locations are 
shown in Figures 5.12 and 5.13, respectively. The flow is directed toward the hub 
over most of the passage in the inviscid region. This is due to the radially inward 
lean of the nozzle trailing edge (see Figure 2.2). Binder and Romey (1982), who 
have a radially inward lean of their nozzle trailing edge, also see radially inward 
flow over most of the span just downstream of the trailing edge of their annular 
nozzle cascade. At X/C = 1.025, very high negative radial flow angles or pitch 
angles (-66 degrees) occur at H=0.83 near the suction side of the wake, while 
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positive angles occur at S= -0.25, at the same radial location. Also at the same axial 
location, high negative radial flow angles (-25 degrees) occur at H=0. 13 and S=0 
degrees while positive angles (21 degrees) occur closer to the hub. This is another 
indication of the presence of vortices, since local high and low radial flow angles 
occur across a vortex. Just as in the yaw angle, the difference between the 
maximum and minimum radial flow angle decreases between X/C = 1 .025 and 
1.09 by 21 degrees in both the hub and casing secondary flow regions. This results 
from the secondary flow vortex decay. The mass-averaged radial flow angles at 
X/C = 1.025 are shown in Figure 5.1 1. The radial flow angles are negative over 
most of the span, reaching zero at the casing. The negative radial flow angles at the 
hub are due to the slot which separates the rotating hub from the stationary hub, 
which is located just downstream of this location. 

5.3 Secondary Flow Velocity 

The secondary flow velocity vectors are shown in Figures 5.14 and 5.15, for 
X/C= 1.025 and 1.09, respectively. The secondary flow vectors in the r-0 plane are 
derived from the measured data and the design flow angle (70 degrees) using 
equation (4.1), as described previously. The presence of strong radial inward flow 
in the wake can be seen. Strong radial inward flow at the wake center near the 
casing is caused by the passage vortex (clockwise) which augments the radial 
inward flow of the wake. The passage vortex can be seen clearly, near H=0.8 at 
X/C= 1.09, while at X/C= 1.025 the radial inward flow is so strong that the vortical 
motion cannot be seen. Since the pressure side leg of the horseshoe vortex has the 
same direction of rotation as the passage vortex, this vortex could be a combination 
of the passage vortex and the pressure side leg of the horseshoe vortex as discussed 



121 


1.0 


0.8 


0.6 


H 


0.4 


SS 


PS 


N . N \ \ » * t t t t f f ////////zVVX 


"" ^ 


////////////// 1 
/////////////// j\ 



****sssssssssssssssssssss //////// / / / / 
**^S***SSSSSSSSSSSS/SSS///S/////////ftJ//ifft 
'S*'SSSs'//SS//////////S /////////////// /////////// 


///////i/i/l nil \'^**$sy 0 y//////////////////////////////////f/////Miii/ 
mutnxww n\N^ ^2^/////////////////////////////////////// / / \ 1 1 1 1 


u'Mii 


‘^//////////////// / / ///////;/// j j / n 1 1 m mini 


h "nm\uu\\\\N^ 


" m J^y///wwnmiiuuitiiiuiimnun\\\\\\ 


* * ■ - *yfy////////////////////////ss"""""‘ 


'^///////////////minimi i immiuunu 


<////////////////// //S/////J///W*" " " 
£ 



0.0 s 0.2 


Figure 5.14. Secondary Flow Velocity Vectors at X/C = 1.025 


0.2511 

m 



122 



Figure 5.15. Secondary Flow Velocity Vectors at X/C =1.09 
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earlier. There is no evidence of the suction side leg of the horseshoe vortex near the 
casing (which would rotate opposite to the passage vortex). Very close to the 
passage vortex, there is a vortical motion in the counterclockwise direction on the 
pressure side of the wake (This will be called the counter-rotating vortex). This is 
formed by the interaction of the passage vortex with the wake. 

The radial flow is directed inward over almost the whole inviscid flow field. 
This is because of the radial inward lean of the nozzle trailing edge, discussed 
earlier, which causes this flow movement toward the hub. The overturning on the 
suction side and the undertuming on the pressure side can be clearly seen at the 
wake center. The radial inward flow is stronger in the wake than in the freestream, 
resulting from an imbalance of the centrifugal force and the pressure gradient close 
to the surface of the blade, which was described earlier. This radial inward flow is 
in the same direction as those induced by the casing passage vortex and thus it 
augments the casing passage vortex. At the hub there is also evidence of the 
passage vortex at X/C= 1.025. But the radial inward flow in the wake near the hub 
is in the opposite direction to those caused by the passage vortex, and thus the 
radial inward flow seems to counteract those induced by the hub passage vortex. 
This is similar to interaction of the radial inward flow of the suction surface 
boundary layer at X/C = 0.935 and the passage vortices. The strong radial inward 
flow of the suction surface boundary layer, observed at X/C = 0.935 (Figure 4.18) 
seems to have disappeared at X/C = 1.025. In fact, the flow on the suction side of 
the wake near the hub has reversed direction and is moving radially outward. This 
is a result of the rotating hub that is located just downstream of this location. The 
rotation of the hub is toward increasing S in this figure. The hub wall boundary 
layer is highly skewed and undergoes sudden perturbation as shown in Figure 4.1. 
At X/C= 1.09 the hub passage vortex seems to have dissipated while the casing 
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passage vortex remains strong. This is also a result of the interaction of the hub 
passage vortex with the rotating hub which is located just downstream of this 
location (see Figure 4.1). This interaction causes the hub passage vortex to 
dissipate. The radial inward flow at the hub is due to the gap between the stationary 
and rotating hub which is at this location. 

Comparing the secondary flow vectors at X/C = 0.935 (Figure 4.18) with 
those at X/C = 1.025, no dramatic change in secondary flow is observed at the tip. 
The pitchwise and span wise position of the passage vortex relative to the wake 
center remains the same as one travels from X/C = 0.935 to X/C = 1.09. This is in 
contrast to the results of Moore and Adhye (1985) where the passage vortex in 
their turbine cascade migrates toward midspan and toward the middle of the 
passage as it progresses downstream. The secondary flow at the hub, however, 
undergoes a dramatic change going from X/C = 0.935 to X/C = 1.09. This is a 
result of the rotating hub which interacts with the secondary flow at the hub. 

The diameter of the passage vortices in the AFTRF nozzle are much smaller 
than the diameter of the passage vortices in other turbine cascades. The maximum 
pitchwise distance of the AFTRF passage vortices downstream of the nozzle are 
25% of the nozzle pitch (0.25S). This is in constrast to the rotor cascades of others 
such as Langson et al. (1977), Moore and Adhye (1985), and Gregory-Smith et al. 
(1988) whose maximum pitchwise width of the passage vortices are 100% of the 
blade pitch. This is also in constrast to the nozzles of Hunter (1982) and Sieverding 
et al. (1984) whose maximum pitchwise width of the passage vortices are 100% of 
the blade width, also. The rotor cascade's passage vortices are larger because the 
rotors have a much larger turning (120 degrees) compared to the AFTRF nozzle 
(70 degrees). The other researcher’s nozzles have larger passage vortices due to the 
larger nozzle inlet boundary layers. 
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5.4 Streamwise Vorticitv 

Even though the isocline angle plots and the secondary flow vector plots are 
useful in identifying vortices, they have their limitations. The isoclinal angle plots 
cannot determine the sign of the vortex, while the secondary flow vectors can 
change depending on what the primary flow is defined to be. A more logical 
approach in identifying both the vortices and their sign is to examine the strength 
of the streamwise vorticity . But this method has its problems, too. The scatter in 
experimental data can cause unrealistic gradients and thus false levels of vorticity. 
But, usually the biggest problem is that the axial distance between the data points is 
much larger than the radial and tangential distances between the data points, thus, 
the streamwise vorticity cannot be found explicitly since it would involve gradients 
of velocity in the axial direction. But in this experiment the axial distance between 
the measurement planes is small (7 mm) and comparable to the tangential distance 
(1 mm) and average radial distance (5 mm) between the data points. 

The components of vorticity in the axial, tangential and radial directions are 
given by equations 4.2 through 4.4 and the streamwise vorticity is calculated from 
equation 4.6 The primary flow angle is set equal to the design exit flow angle (70.0 
degrees). 

The contour plots of streamwise vorticity (< o s C/V x ) at X/C = 1 .025 and X/C 

= 1.09 are given in Figures 5.16 and 5.17, respectively. The existence of a strong 

casing passage vortex is clearly seen in these figures. The large positive vorticity 
(peak magnitude of (D S C/V X] = 17 at X/C = 1.025 and co s C/V = 18 at X/C = 1.09) 

near the casing corresponds to the passage vortex. At the hub, the passage vortex is 
revealed by the presence of large negative vorticity. The magnitude of peak 
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vorticity at the center of the vortex at X/C = 1.025 (cd s C/V Xi — -15) is similar in 
magnitude to the peak vorticity at the casing, which demonstrates that the hub 
passage vortex is similar in strength to the casing passage vortex at X/C = 1.025. 
The magnitude of peak negative vorticity at the hub decreases by 33% from 
X/C= 1.025 to X/C=1.09 which agrees with the secondary flow vectors that the hub 
passage vortex is decreasing in strength as it moves downstream of the nozzle 
trailing edge. 

At X/C = 1.025, there are two negative vorticity regions near the casing. 

The larger (in area) of the two (peak magnitude of <d s C/V x = -7) is due to the 

counter-rotating vortex (which is caused by the interaction of the passage vortex 
and the wake). The smaller negative vorticity region (peak magnitude of co s C/V x = 

-11) located just below the casing passage vortex, possibly could be the suction 
side leg of the horseshoe vortex. At X/C = 1.09, there is only one negative vorticity 
region near the casing. This vortex is also due to the interaction of the passage 
vortex with the wake. The small negative vorticity region seen at X/C = 1.025, 
seems to have disappeared. It is possible that this vortex has either merged with the 
large negative vorticity region between the two axial planes or has been dissipated 
by its contact with the casing passage vortex. 

Near the wake center from H=0.7 to H=0.1, there is a positive vorticity 
region which corresponds to the over and undermining regions in the wake. There 
is a positive vorticity peak at H = 0. 1 , which is caused by the interaction of the 

passage vortex and the wake. The peak magnitude of this positive vorticity 
decreases from 0) C/V = 25 at X/C = 1 .025 to CD C/V = 1 1 at X/C = 1.09. This is a 

S X] 5 Aj 

also a result of the interaction of the flow near the hub with the rotating hub. 

Comparing the streamwise vorticity upstream of the trailing edge (X/C = 
0.935) which is presented in Figure 4.20 with the vorticity at X/C = 1.025, the 
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magnitude of the casing passage vortex has decreased by half (from oo s C/V x = 30 
to 17). This is due to the interaction of the passage vortex and the wake, the 
interaction of the passage vortex and the freestream flow and also the decay of the 
passage vortex as it travels downstream of the nozzle. On the other hand, the 
negative vorticity region associated with the hub passage vortex has not changed 
very much in strength between X/C = 0.935 (( 0 s C/V Xi = -14) and X/C = 1.025 
(o) s C/V x - -15). Whereas at X/C = 0.935, the negative vorticity region corresponds 

to only the hub passage vortex, at X/C = 1 .025 the negative vorticity region results 
from both the hub passage vortex and the interaction of the endwall flow with the 
rotating hub. 

The streamwise vorticity was also calculated using the method of Gregory- 
Smith (1988), discussed in Chapter 4, which approximates the axial velocity 
gradients with an inviscid form of the Navier-Stokes equations. 

The pitch-averaged vorticity (m s C/V I ) plots are shown in Figure 5.18 for 

X/C= 1.025 and 1.09. Outside of the secondary flow regions the pitch-averaged 
vorticity is positive due to the over- and under-turning regions in the wake. A large 
positive vorticity region occurs at H=0.1, which is a result of the vortex formed by 
the interaction of the passage vortex and the wake. Moving closer to the hub, the 
vorticity becomes negative due to the existence of the passage vortex. The 
maximum negative vorticity occurs next to the hub and is a result of the interaction 
of the flow with the rotating hub. The variation in vorticity is not as large near the 
casing as it is near the hub, because the passage vortex and the counter-rotating 
vortex are side by side (at the same radial location) near the casing instead of at 
different radial locations as they are near the hub. The high positive vorticity at 
HMD. 85 is due to the passage vortex. 

The experimentally derived pitch-averaged secondary vorticity is compared 
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to the theoretical secondary vorticity of Squire and Winter (1951) (derived using 
classical secondary flow theory) in Figure 5.19. The experimental pitch-averaged 
secondary vorticity co^is defined as 


«*» = <**-«*. (5.1) 

where co s is the pitch-averaged streamwise vorticity at each radii and <5. m is the 

pitch-averaged streamwise vorticity at midspan. The secondary vorticity of Squire 
and Winter (Gl> ) is calculated as follows: 

x secsw ' 

=2to Di (a 2 -a,) (5.2) 

where co n is the vorticity at the nozzle inlet, a, is the inlet flow angle and a 2 is the 

mass-averaged exit flow angle. The agreement between the theoretical vorticity 
and the experimentally derived vorticity is good at the hub, while at the casing 
Squire and Winter's formula over predicts the secondary vorticity. This is 
somewhat unexpected, since the experimental secondary vorticity co^at the casing 
only contains the vorticity due to the secondary flow, while at the hub co^ contains 

both the vorticity due to secondary flow and due to the rotating hub. Therefore it 
must be just fortuitous that the prediction of the secondary vorticity matches the 
experimental secondary vorticity at the hub. Thus one can conclude that the 
theoretical secondary vorticity does not correctly predict the actual secondary flow 
vorticity in a turbine cascade. This is expected, since one of the assumptions in 
deriving Squire and Winter's secondary vorticity is that the turning angle (a 2 - a,) 

is small and this is violated in turbine cascades. The experimentally derived 
secondary vorticity was also compared to the secondary vorticity derived by Came 
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and Marsh (1974), but this vorticity is similar to that of Squire and Winter's and 
thus it is not shown. 

Comparison of the peaks in streamwise vorticity with the total pressure loss 
contours (Figures 5.1 and 5.2) shows that the peak total pressure loss does not 
occur at the peak vorticity location, but occurs somewhere between the positive 
and negative vorticity peaks. This is contrary to the results of linear cascades where 
the loss peaks coincide with the vorticity peaks (vortex centers), (see Gregory- 
Smith et al. (1988), and Hodson and Dominy (1987)). The interaction between the 
secondary flow vortices and the wake causes the accumulation of the low 
momentum and energy endwall boundary layer fluid into these two locations. 

The six criteria that are used to determine the existence of a vortex in the 
flow field were described in Chapter 4. A vortex exists when these six flow 
phenomena occur at the same position in the flow field. The existence of the 
passage vortices at both the hub and the casing are confirmed since they meet all 
six criteria. On the other hand, the existence of the suction side leg of the horseshoe 
vortex is more doubtful. There is no evidence for its existence at the hub. While at 
the casing there is a region of high negative vorticity at X/C= 1.025 that could 
correspond to the suction side leg of the horseshoe vortex, the other criteria are not 
met. Thus its existence cannot be proven conclusively, and has probably decayed 
by the time it has reached the nozzle trailing edge. This is in agreement with the 
findings of other researchers, both experimental (Gregory-Smith et al., 1988) and 
computational (Ho and Lakshminarayana, 1994), who found that the suction side 
leg of the horseshoe vortex had decayed by the time it had reached the trailing 
edge. 
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5.5 Nozzle Wake Characteristics 


An understanding of turbine nozzle wakes are important for the efficient 
design of axial flow turbomachinery. A major cause of noise and vibration 
characteristics of turbomachinery is caused by wakes. Turbine wakes represent a 
source of loss in efficiency, since the mixing of the wakes with the freestream 
dissipates energy. The three-dimensional characteristics of the wake, the decay 
characteristics and the path that it follows is important in the design of the 
following blade rows. This information is essential for both the prediction of the 
aerodynamic and mechanical performance of a turbine and for building quieter 
turbomachines. An understanding of the wake development and its decay is also 
essential because of the role it plays in the rotor-stator interaction . 

5.5.1 Nozzle Wake Static Pressure 

The static pressure variation across the wake for different radii at X/C = 
1.025 is shown plotted in Figure 5.20. Starting from the pressure side, the static 
pressure coefficient decreases sharply until the wake center, after which it increases 
rapidly across the suction side of the wake. A decrease in static pressure 
coefficient indicates an increase in static pressure as this is an accelerating flow. 

A hump at the suction side wake edge (A) has been caused by the over-turning on 
the suction side, while the dip in static pressure coefficient on the pressure side 
wake edge (B) is caused by the under-turning on the pressure side of the wake. 

The strong radial pressure gradient is evident with the static pressure 
decreasing more than 35% as one goes from the tip to the hub. The static pressure 
variation across the wake is as high as 20-25% of the value in the freestream. This 




Figure 5.20. Static Pressure (4^) Variation at X/C = 1.025 
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is consistent with the measurements of Lakshminarayana and Davino (1980) who 
found similar static pressure gradients across the wake of a compressor inlet guide 
vane. 

The existence of pressure gradients across the wake has been explained for 
a compressor rotor blade by Ravindranath and Lakshminarayana (1980), and can 
be comprehended by examining the equation of motion in a streamwise (s) and 
normal (n) coordinate system. The equation can be approximately written as 


i dp v s 2 3 r , 2 ^ 

pdn R c + dnl V " J 


( 5 . 3 ) 


where n is the direction normal to the streamline, V s is the streamwise velocity, R c 
is the radius of curvature of the streamline and (v') is the turbulent fluctuation in 

the n direction. It is clear from the above equation that in addition to the centrifugal 
force, the gradient of turbulent intensity in the n direction has a major influence on 
the pressure gradient dp/dn. The static pressure gradient is caused by the flow 

curvature, velocity change, turbulence intensity and possibly separated flow in the 
trailing edge region. 

5.5.2 Total Velocity 


Figures 5.21 and 5.22 show the total velocity profiles for X/C=1.025. The 
plot indicates that the suction surface boundary layer is thicker than the pressure 
surface boundary layer. Several interesting observations can be made. The wake at 
H=0.98 is well behaved. The influence of interaction between the wake and the 
secondary flow can be seen in the wake data at H=0.870 to H=0.810. The two 
distinct troughs from H=0.870 and 0.840 indicate that the wake and secondary 
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Figure 5.21. Total Velocity (V/U m ) Profiles at X/C = 1.025 
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Figure 5.22. Total Velocity (V/U m ) Profiles at X/C = 1.025 



139 


flows are still distinct and they are likely to merge as the flow progresses 
downstream. The region near H=0.810 is the location of the strongest secondary 
flow/vorticity and here, there is clear evidence of interaction between the wake and 
vortex resulting in a thick shear layer on the suction side. The region from H=0.766 
to 0.150 has a well behaved wake and the secondary flow/vortex interaction with 
the wake is again evident below this region. In general, the wake defect is higher in 
the secondary flow regions. 

5.5.3 Nozzle Wake Decay Characteristics 

A knowledge of the rate of decay for the wake defect is necessary for an 
understanding of the rotor-stator interaction. The decay of the velocity defect is 
influenced by many phenomena, including the pressure gradient, turbulence 
intensity, curvature and viscous effects. The endwall secondary flow and passage 
vortices also have an effect on the wake decay. The velocity defect plotted in 
Figure 5.23, is based on the maximum and minimum velocity in the wake, plotted 
against Z/coscc 0 , where a o is the local value of the vane outlet angle. This 

corresponds approximately to the streamwise distance downstream of the nozzle. 

The authors' data at midspan is compared with data from other sources. 
Goldman and Seasholtz (1982), (1992) did their measurements in an annular 
cascade, while Dring et al. (1987) measured the wake behind the nozzle with a 
rotor-stator spacing of 50% of nozzle axial chord. Sitaram and Govardhan (1986) 
did their measurements in a linear turbine rotor cascade. Ho and Lakshminarayana 
(1994) computed the flow field in a linear turbine rotor cascade. Raj and 
Lakshminarayana (1973) measured the wake behind a linear cascade of compressor 
blades and developed a correlation. They found that the wake from a compressor 
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Streamwise Distance, Z/cosa 0 


Figure 5.23. Decay of Total Velocity Defect with Streamwise Distance (Midspan) 
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cascade decays slower than that of a flat plate, cylinder or a symmetrical airfoil. 
They attributed this to the fact that the wake edge velocity decelerates in a 
compressor cascade, thus causing an adverse pressure gradient, which causes the 
wake to decay slower. 

The nozzle wake decays rapidly close to the trailing edge and less rapidly 
farther downstream. This can be seen more clearly in Figure 5.24, which shows the 
nozzle wake at midspan for several axial locations downstream of the nozzle. 
Lakshminarayana and Davino (1980) state that the rapid decay rate close to the 
trailing edge is due to pressure gradients, high turbulence intensities and wake 
centerline curvature. The wake data from Dring's cascade and Goldman's cascade 
decay slightly faster than Raj's cascade. This is because the wake edge velocity 
decay is almost zero for the turbine nozzle and rotor cascades (Figure 5.25), hence, 
the streamwise pressure gradient in the edge of the wake is zero. And according to 
Hill et al. (1963), wakes in an adverse pressure gradient will decay slower than 
wakes in the presence of zero pressure gradient. In addition, the wake in the 
presence of favorable pressure gradient is found to decay faster. Also, the 
increased loading on the nozzle vanes compared to Raj and Lakshminarayana's 
compressor cascade could cause an increase in the wake decay rate. 

The data from the AFTRF nozzle, on the other hand, have a much more 
rapid decay rate than the compressor cascade wake or Goldman and Seasholtz's 
turbine nozzle cascade. This is due to the presence of the rotor downstream at a 
very small rotor-stator axial spacing (rotor-stator spacing is 20% of the nozzle axial 
chord). The relative motion between the rotor and the stator causes periodic 
variations in the potential flow field around the blades. This unsteadiness causes 
the wake to decay faster than the wake of a cascade with no rotor behind it. The 
nozzle wake is also affected by the favorable pressure gradient in the rotor, which 
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Figure 5.24. Total Velocity Profiles at Midspan 
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Figure 5.25. Variation of Wake Edge Velocity with Downstream Distance 
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causes it to decay faster. Even though the nozzle of Dring et al. (1987) has a rotor 
downstream, the rotor-stator spacing is much larger (50% nozzle axial chord) than 
that of the AFTRF nozzle. Thus the potential flow interaction will have less effect 
on the nozzle wake in Dring's et al. (1987) case. Ho and Lakshminaray ana's turbine 
rotor cascade (computed) also decays faster than the compressor wake or Dring's or 
Goldman's nozzles. This is due to the fact that the loading on the rotor was much 
higher than on the nozzles, since the rotor turning is 1 10 degrees while the turning 
is around 70 degrees for all the nozzle data shown in Figure 5.23. Sitaram and 
Govardhan's turbine rotor wake decays a little faster than Ho and 
Lakshminarayana's rotor cascade since their loading is higher (rotor turning is 1 20 
degrees for Sitaram and Govardhan's cascade). 

The decay of maximum (absolute) radial velocity (Figure 5.26) normalized 
by the maximum velocity in the wake is constant from the trailing edge until 
Z/cosa 0 = 0.30, and then decreases slightly. The magnitude is similar to the radial 

velocity decay of a compressor stator reported by Lakshminarayana and Davino 
(1980). Figure 5.27 shows the radial variation of total velocity defect at the 
centerline at X/C = 1.025 and X/C = 1.09. Caution should be used in interpreting 
this data. Since the wake belongs to different blade sections with varying in 
boundary layer growth, the variation of defect in the radial direction should not be 
interpreted as the decay rate. The variation of wake defect with spanwise distance 
is clearly seen from this plot. There are two regions where the wake defect is larger 
than the others, one centered around H=0. 1 and the other centered around H=0. 8 . 
This is the region where the secondary flow vortices are located, and the 
interaction of these with the wake results in deeper and wider wakes. Thus the 
wake defect will be larger at these locations. The wake decay between X/C= 1.025 
to 1 .09 is almost constant from the hub to H=0.6 and then increases until the 
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Figure 5.26. Decay of Radial Velocity 
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maximum decay rate is reached at H=0.8. This is due to effect of the two counter- 
rotating vortices at this location. 

The variation of the semi-wake width at midspan with streamwise distance 
is shown in Figure 5.28. The semi- wake width is defined as the width of the wake 
at half the defect of total velocity. The semi-wake width grows rapidly close to the 
trailing edge and more gradually farther downstream. This growth is due to the 
exchange of momentum, mass and energy on both sides of the wake. The variation 
of the maximum static pressure difference ( *P d ) across the wake at midspan is 
shown vs. streamwise distance in Figure 5.28. The value of T,, decays to 50 
percent its value within a streamwise distance of 10 percent chord. This rapid 
decay is caused by the intense mixing and high turbulence intensities that occur 
close to the trailing edge (see Ravindranath and Lakshminarayana (1980)). The 
decay rate levels off as the wake moves downstream of the trailing edge reaching a 
value of 28 percent of its value at the trailing edge within a streamwise distance of 
40 percent of chord. Most researchers assume that the static pressure is uniform at 
the trailing edge, both inside the wake and in the free stream. This assumption is 
not valid in view of the data shown in Figures 5.28 and 5.20. 

5.6 Mass- Averaged Properties 

The circumferentially mass-averaged total pressure loss profiles at various 
locations from upstream of the nozzle to downstream of the nozzle are shown in 
Figure 5.29 (The properties upstream of the nozzle inlet for all the following 
figures are a single radial traverse). The stagnation pressure profile at the nozzle 
inlet shows the boundary layer thickness to be approximately 10% at the hub and 
15% at the tip. At X/C=0.935, 1.025 and 1.09, the upstream boundary layer profile 
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o 



Figure 5.28. Variation of Static Pressure and Semi-Wake Width with Streamwise 
Distance 
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D X/C = -1.0 
® X/C = 0.935 
a X/C = 1.025 
+ X/C = 1.09 
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Figure 5.29. Mass-Averaged Total Pressure Loss ('Floss) 
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losses have been subtracted, hence, those represent the losses occurring within the 
nozzle passage. At all three axial locations, there are two loss peaks. The peak near 
the hub is larger in magnitude but smaller in radial extent them the one near the 
casing. The larger inlet casing boundary layer results in the larger extent of the 
secondary flow region near the casing. The larger magnitude of the hub loss region 
is a result of the radial inward flow of the boundary layer and wake which cause 
the low momentum fluid in the boundary layer and wake to move toward the hub. 

The mass-averaged static pressure coefficient is shown in Figure 5.30. The 
radial static pressure gradient at the locations X/C = 0.935 to X/C = 1 .09 is nearly 
constant from hub to tip and is consistent with the velocity distributions shown 
later. 

The mass-averaged total velocity is shown in Figures 5.31. The boundary 
layer profile is clearly visible upstream of the nozzle inlet. At X/C = 0.935 to X/C 
= 1.09, the total velocity profiles show a sharp decrease in velocity near H = 0.10 
and a more gradual decrease in velocity near H = 0.80, which is due to the hub and 
tip wall secondary flow. The total velocity predicted by Katsanis' code is compared 
with the total velocity measured at X/C = 0.935 and X/C = 1.025. The agreement is 
excellent. 

The mass-averaged yaw angle is given in Figure 5.32. The yaw is derived 
from the mass averaged tangential velocity and the mass averaged axial velocity. 
At the nozzle inlet (X/C = -1.0), the yaw angle is about zero everywhere. The yaw 
angles at X/C = 0.56 show the presence of weak secondary flow. From X/C = 

0.935 to X/C = 1.09, the yaw angle shows overturning at the casing and hub 
endwalls. As one moves away from the endwalls, the flow becomes undertumed, 
achieving design flow near 30% span. This is consistent with the predictions from 
secondary flow theories and measurements by others. The largest overturning is 
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Computation (Katsanis' code) 

13 X/C = -1.0 
o X/C = 0.56 
a X/C = 0.935 
+ X/C = 1.025 
x X/C = 1.09 



V/U m V/U m V/U m V/U m V/U m 


Figure 5.31. Predicted and Measured Mass-Averaged Total Velocity ( V/U m ) 
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Computation (Katsanis' code) 

□ X/C = -1.0 
© X/C = 0.56 
A X/C = 0.935 
+ X/C = 1.025 
x X/C = 1.09 


o 



a (°) a (°) a (°) a (°) a (°) 


Figure 5.32. Predicted and Measured Mass- Averaged Yaw Angles (a) 
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about five degrees, and this is much smaller than those reported by others for rotor 
cascades. The yaw angle at X/C = 0.935 and 1.025 is compared to Katsanis' 
in viscid code predictions. The agreement is very good everywhere except very 
close to the hub. The prediction even captures the overturning at the casing. This 
demonstrates the usefulness of using an inviscid code in turbomachinery design. 

Figure 5.33 shows the area-averaged stagnation pressure loss for the PSU 
AFTRF and two other turbine cascades. In order to compare the losses between the 
three turbine blades, the stagnation pressure loss for all three are normalized by the 
inlet dynamic head. The AFTRF nozzle follows the same trend as the other turbine 
blades, increasing as one moves downstream. The large increase in losses going 
from X/C = 0.935 to X/C = 1.025 is partly due to separation that occurs at the 
vane's trailing edge and partly to the fact that the vane surface boundary layer 
could not be completely measured close to the wall (thus giving a lower area 
averaged loss), while downstream of the nozzle the losses in the entire wake could 
be measured. The loss growth downstream of the nozzle trailing edge is caused by 
the mixing of the flow downstream of the vane and by turbulence mixing. 





156 


CHAPTER 6 

STEADY AND UNSTEADY FLOW FIELD AT ROTOR MIDSPAN 

The flow around turbine blades is highly unsteady. This unsteadiness is 
caused by the aerodynamic interaction of the nozzle and rotor flow fields and is 
called rotor-stator interaction. Rotor-stator interaction can affect the aerodynamic, 
structural, thermal and acoustic performance of a turbine and is important for the 
efficient design of axial flow turbines. According to Dring et al. (1982), rotor-stator 
interaction can be divided into two parts. These are potential flow and wake 
interactions. The potential flow gradients extend both upstream and downstream of 
the blade and they decay exponentially with a length scale of the order of the pitch 
or chord of the blade row. If the axial gap between the blade rows is less than a 
chord (which it is in a typical axial flow turbine), then the potential influence can 
cause unsteadiness in both upstream and downstream of the blades. 

However, the wake is convected downstream and has a far field rate of 
decay much lower than that of the potential flow. The wake will still be felt several 
chords downstream. But in most modem axial flow turbines, which have a rotor- 
stator spacing close to 20 percent of a blade chord, both the potential and wake 
effects will occur together. And in the future as gas turbine designers try to reduce 
weight, and thus the rotor-stator spacing gets smaller, these effects will become 
more prominent. 

Even though unsteady flow plays a major role in axial flow turbines, 
turbines are designed using three-dimensional steady flow calculation methods 
(Sharma et al., 1992). Empirical correlations are used to account for the effect of 
the unsteadiness. Because actual models of the loss generating mechanisms in 
unsteady flow turbomachinery do not exist, these correlations are based on results 
from stationary cascade data and do not represent the actual fluid mechanics in the 
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flow field (Sharma et al., 1985). Thus these correlations must be multiplied by 
some factor to obtain a good estimate of the actual losses that occur in turbines. 
Although these correlations have worked well in the design of existing turbines, 
they do not represent the true physics of the flow field and are only useful in the 
areas from which they were obtained, namely design point predictions and turbines 
which are similar to existing designs (Hathaway, 1986). Thus, a more thorough 
knowledge of unsteady flow interactions is needed in order to increase both the 
design and off-design performance of existing turbines, and to design turbines 
which are considerably different than existing turbines. To obtain this knowledge, 
good time accurate data from inside the rotor is needed. This knowledge, in turn, 
can be used to model the unsteady flow mechanisms that are not currently in 
existing design codes. Thus, a better understanding of unsteady flow interactions 
can lead to an improvement in the ability to predict the performance of turbines 
and to corresponding improvements in the actual performance of turbines. 


6.1 Measurement Procedure 

Figure 6.1 shows the LDV measurement locations in the rotor. 
Measurements have been acquired at 37 axial measurement locations from just 
upstream of the rotor (X r /C r = -0.088) to one chord downstream of the rotor. Table 

6.1 lists these locations. Each measurement location is at midspan. Since a two- 
dimensional LDV was used for measurements in the rotor, only the velocities in 
the axial and tangential directions are measured. To account for the nonuniformity 
of the rotor absolute inlet flow field, measurements were made at six tangential 
locations in the absolute frame equally spaced over one nozzle pitch. These six 
tangential locations represent six different relative positions between the nozzle 
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Table 6.1 

Axial Measurement Locations in the Rotor 


X ,/ C , 

X ,/ C , 

- 0.088 

0.992 

- 0.063 

1.004 

- 0.042 

1.014 

- 0.021 

1.024 

- 0.011 

1.035 

0.000 

1.066 

0.010 

1.116 

0.040 

1.170 

0.091 

1.220 

0.142 

1.270 

0.194 

1.320 

0.297 

1.370 

0.398 

1.420 

0.502 

1.470 

0.604 

1.527 

0.706 

1.640 

0.808 

1.800 

0.911 

1.987 


0.962 



and the rotor (labeled nozzle/rotor locations 1 through 6) or if viewed from the 
nozzle frame of reference, six different time resolved positions of the rotor in 
relation to the nozzle. 
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As discussed in the LDV data analysis section (section 2.6), each rotor 
passage is divided into 50 measurement windows, which means that there are 
measurements at 50 circumferential locations across the rotor pitch. The velocity is 
then ensemble averaged at each measurement window according to the following 
equation 

V-i£v, (6.1) 

n i= i 

where V is the ensemble averaged velocity, V, is the instantaneous velocity 
measured at a particular rotor measurement window and n is the total number of 
measurements in that measurement window. The unresolved velocity for each 
measurement window can also be calculated as 

V'= V - V (6.2) 

and the corresponding variance as 

[£(v,-v) ! ] 

^ -Vi) (63) 

The level of unresolved unsteadiness in each measurement window is determined 
by the variance. 

Since the flow field between the rotor passages was demonstrated to be 
periodic (see section 2.6 in this thesis), all the ensemble average velocity 
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components were spatially phase-locked averaged as follows 




RB 


(6.4) 


where G represents either the ensemble averaged velocity or variance, the subscript 
m determines the particular rotor passage and j is the measurement window 
location relative to the m'th rotor passage and N,^ is the number of rotor blades. 

The successive application of equation (6.4) results in a description of the flow 
field at 50 equally spaced shaft positions across a representative rotor passage 
(Figure 2.10). Since all results presented from now on are spatially phase-locked 
averaged, the superscript u will be dropped hereafter. 

For the LDV measurements in the rotor, the instantaneous velocity, V, , is 
decomposed as follows 

V. = V + V + V' (6.5) 

where V is the time averaged velocity, V is the periodic velocity and V' is the 
unresolved velocity component as calculated in equation (6.2). The decomposition 
for the phase-locked averaged rotor blade passage is shown in Figure 2.10. The 
time averaged velocity V is obtained by averaging all the ensemble averaged 
velocities in each measurement window as follows. 



where NW is the total number of measurement windows. The time averaged 
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velocity is a time average of all measurements acquired at a fixed point in space. 
The periodic velocity is then obtained from the equation, 

V=V-V (6.7) 

From the above velocity decomposition, the axial, tangential and cross velocity 
correlations (both periodic and unresolved) can be computed, also. 

Figure 6.2 shows the rotor inlet velocity triangles. The velocity defect in the 
nozzle wake produces a slip velocity in the relative frame towards the rotor suction 
surface. A more detailed discussion of this phenomenon is provided later in this 
chapter. 


6.2 Cycle Average Properties 


The cycle averaged values are obtained by averaging the ensemble averaged 
properties in each rotor measurement window for one nozzle/rotor location over 
the six nozzle/rotor locations as follows: 


1 ^NBP 

G=— 


N 


NRP k=l 


( 6 . 8 ) 


where G represents any flow parameter (such as velocity or unresolved 
unsteadiness), the superscript ... stands for cycle averaged property, N nrp is the 

number of nozzle/rotor positions (6), the subscript j represents the individual 
measurement window location in the rotor passage, and the subscript k = 1 to 6 
denotes a particular nozzle/rotor position. All the rotor flow field contour plots 
presented in this thesis are for one ensemble averaged rotor blade pitch which is 
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Figure 6.2. Rotor Inlet Velocity Triangles 
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doubled so that two rotor pitches are shown in the contour plots. Since the LDV 
was oriented at 7.6 degrees to the circumferential direction (in order to reduce 
reflections from the laser beam on the casing window from saturating the 
photomultipliers), there is a shadow region on the blade suction side where no 
measurements could be obtained. 

The cycle averaged LDV data at the farthest measurement location 
upstream of the rotor (X r /C r = -0.088) was mass averaged over one rotor pitch and 

compared the mass averaged five hole probe data closest to this location ( X r /C r = - 

0.080 or X/C = 1.16). The agreement is excellent with the absolute total velocity 
being 0.98% of each other ( v/U m =1 .61 19 for the five hole probe and 
v/lJ m =1.6278 for the LDV) and the absolute flow angles being within 0.30% of 
each other (a=70.42 degrees for the five hole probe and a =70.63 degrees for the 
LDV). In Chapter 6, the superscript n represents mass-averaged properties, in 
order to differentiate mass-averaged from ensemble averaged properties. (The 
comparison of LDV and five hole probe data is presented in section 6.4 of this 
thesis.) 

6.2.1 Relative Total Velocity and Total Unresolved Unsteadiness 

The cycle averaged relative total velocity normalized by mean rotor speed 
(U m ) is presented in Figure 6.3a. (The bar over the velocity and angle notation in 

the figures in Chapters 6 and 7 represents ensemble averaging.) This figure shows 
the flow acceleration from the pressure to the suction surface of the blade. The 
velocity accelerates gradually on the suction side from the leading edge to X r /C r = 

0.80, after which it levels off and becomes fairly uniform until the trailing edge. On 
the pressure side the velocity change is fairly gradual from the leading edge to the 
trailing edge. The effect of the leading edge on the flow field is clearly shown with 
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the flow decelerating as the leading edge is approached. The relative total velocity 
decelerates to a value of w/u m = 0.2 just upstream of the leading edge. This can 

be seen more clearly in Figure 6.3b which presents an enlargement of the flow field 
near the rotor leading edge. Downstream of the blade, the rotor wake decays to 
negligible values within half a chord length downstream of the blade. The rotor 
wake thickness is small. The rotor wake is discussed in more detail in Chapter 7 of 
this thesis. Figures 6.4a, b and c show the blade-to-blade cycle averaged velocity 
profiles at three axial locations inside the rotor passage compared to the design 
blade surface velocity. Extrapolating the measured velocity to the blade surfaces, 
the design velocity matches the measured velocity well near the rotor leading edge. 
At midchord and near the trailing edge the match is not as good due to the 
boundary layer on the blade surfaces, which causes the measured velocity to be 
lower near the blade surface than the design velocity (The design velocity was 
calculated using an inviscid code.) 

The cycle averaged relative total velocity vectors are presented in Figure 
6.5. The measured velocity vectors follow the rotor blade contour almost 
everywhere, except near the leading and trailing edge. The potential effect can 
clearly be seen at the leading edge, while the over and undertuming of the flow 
field on the suction and pressure side, respectively, at the trailing edge can also be 
seen clearly. 

The cycle averaged relative total unresolved unsteadiness is shown in 
Figure 6.6. The total relative unresolved unsteadiness is defined as follows: 


Tu. = 


u /2 + v ' 2 )/2 


W 


x 100% 


(6.9) 


where 
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Figure 6.5. Cycle- Averaged Relative Total Velocity Vectors 
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Level Tu r (%) 
D 100.00 
C 50.00 

B 30.00 

A 20.00 

9 18.00 

8 16.00 

7 14.00 

6 12.00 

5 10.00 

4 8.00 

3 6.00 

2 4.00 

1 2.00 


Figure 6.6. Cycle-Averaged Relative Total Unresolved Unsteadiness Tu 
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and 



_ £(v.,-v,) ! ] 

#2 L i=I J 


v' 2 = 


(n — 1) 


( 6 . 10 ) 


( 6 . 11 ) 


where V is the ensemble averaged velocity, V, is the instantaneous velocity 
measured at a particular rotor measurement window and n is the total number of 
measurements at that measurement window. The relative total unresolved 
unsteadiness is then cycle averaged using equation (6.8). The relative unresolved 
unsteadiness is low inside the rotor passage and high in the rotor wake. The highest 
level is observed at the leading edge. This is demonstrated more clearly in Figure 
6.7 which is a blade-to-blade profile of cycle average total relative unresolved 
unsteadiness just upstream of the rotor leading edge. Figure 6.8 shows the total 
cycle-averaged unresolved unsteadiness at the leading edge with the total 
unresolved unsteadiness (Tu t ) normalized by the mean rotor speed instead of the 
local relative velocity. The total unresolved unsteadiness is defined as follows: 

J(u^ + v^)/2 

Tu, = x 100% (6.12) 


The total unresolved unsteadiness is then cycle averaged using equation (6.8).The 
high levels of total unresolved unsteadiness are still visible at the leading edge. 
This demonstrates that the increase in relative total unresolved unsteadiness at the 
rotor leading edge shown in Figures 6.6 and 6.7 is not only due to low relative 
velocities, but also due to an increase in the absolute magnitude of the velocity 
fluctuations. This is in contrast to measurements made at the leading edge of a 
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rotor cascade by Priddy and Bayley (1988), who showed that the absolute 
magnitude of the velocity fluctuations slightly decreased as the leading edge was 
approached. But it agrees with the results of Hobson and Shreeve (1991), who 
measured an increase in the absolute magnitude of the velocity fluctuations at the 
leading edge of a compressor cascade. This can be explained by examining the 
turbulent kinetic energy equation, which is, 


dk+U jdk^ IL v > ^ 

dt 1 dx; dXj ^G k dXj 


dU. 

1 UjUj dx 


— £ 


(6.13) 


I 


n 


m 


iv v 


where term I is the rate of change of kinetic energy (k), term D is the convection 
term , III is the diffusion term, IV is the production term and V is the dissipation. 
The kinetic energy is a measure of the intensity of the velocity fluctuations. As the 
leading edge is approached, the large increase in mean flow velocity gradients 
cause an increase in production term which overshadows the dissipation, diffusion 
and convection terms. Thus, the total unresolved unsteadiness should increase as 
the leading edge of the blade is approached. In addition, the total unresolved 
unsteadiness is observed to be higher near the pressure side. This is because the 
unresolved unsteadiness in thin shear layers is highly sensitive to streamline 
curvature in the plane of the mean shear. The unresolved shear stress and 
unsteadiness are increased due to curvature when the angular momentum of the 
flow decreases in the direction of the radius of curvature. 
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6.2.2 Axial Velocity and Relative Flow Angle 

Figure 6.9 and 6.10 show the cycle averaged axial velocity and relative flow 
angle contours, respectively. The potential effect of the blade is seen to have a 
large effect on the axial velocity upstream of the leading edge, which is shown 
more clearly in Figure 6.11, which is a blade-to-blade cycle averaged axial velocity 
profile just upstream of the leading edge. The deceleration of the flow field in the 
presence of the stagnation point at the rotor leading edge and the acceleration away 
from it are dramatically illustrated in this figure. The axial velocity in the 
stagnation region is 30% of the free stream value at this location. Inside the 
passage, the axial velocity increases almost linearly from the pressure to the 
suction side. The relative flow angle shows the characteristic large change in flow 
angle for a turbine rotor, with the flow being turned around 110 degrees in the rotor 
passage. The effect of the leading edge on the flow field is also evident by the large 
change of flow angle at the leading edge. The blade-to-blade profiles of the cycle 
averaged relative flow angle presented in Figures 6.12a and b show that the rotor 
leading edge has a significant effect on the flow field even 9% of the rotor axial 
chord upstream of the leading edge, with the change in flow angle being 18 degrees 
across the passage at this location. This effect increases as the rotor leading edge is 
approached with the change in angle increasing 35 degrees just upstream of the 
rotor leading edge. 

6.2.3 Unresolved and Periodic Unsteady Velocity Correlations 

Figures 6.13 and 6.14 show the cycle averaged unresolved unsteady 
velocity cross correlations for both the turbine coordinate system (uV) and 
streamwise-normal coordinate system (u'v' ). The unresolved velocity correlations 
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Figure 6.14. Cycle-Averaged Unresolved Velocity Correlation 
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are defined as follows: 


r-. K, - V,K - V.) 

r 

ui 


U V = 


x 100% 


(6.14) 


v'v'(cos 2 p-sin 2 Pj + ^u 72 - v ,2 )cos(3sinp 


U V — ■ 

s n 


u 


x 100% 


(6.15) 


where 

<^ = ( v ,-v < X v ,-v.) 

and the subscript i represents instantaneous velocity, the superscript — represents 
the ensemble averaged velocity in each rotor measurement window and P is the 
relative flow angle. The unresolved unsteady velocity cross correlations for both 
the turbine coordinate system (uV) and streamwise-normal coordinate system 
(uV ). are then cycle averaged using equation (6.8). The unresolved cross 

correlations are zero almost eveiywhere except at the leading edge and near the 
pressure surface of the rotor and in the rotor wake. This trend is expected since the 
production of the unresolved velocity cross correlations is brought about by 
velocity and turbulence intensity gradients. While the axial-tangential unresolved 
cross correlation is negative in the rotor wake, the streamwise-normal unresolved 
cross correlation is positive on the pressure side of the wake and negative on the 
suction side of the wake. The positive and negative streamwise-normal velocity 
cross correlation distribution about the wake center results from opposite gradients 
of streamwise mean velocity about the wake center. 

The cycle averaged periodic velocity correlations are presented in Figures 
6.15, 6.16 and 6.17. The periodic velocity correlations are defined as follows: 
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Level 

u u (%) 

F 

24.34 

E 

22.72 

D 

21.10 

C 

19.47 

B 

17.85 

A 

16.23 

9 

14.60 

8 

12.98 

7 

11.36 

6 

9.74 

5 

8.11 

4 

6.49 

3 

4.87 

2 

3.25 

1 

1.62 


/ — \ 

uu 

V 


Figure 6.15. Cycle-Averaged Periodic Velocity Correlation 
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Level vv(%) 
F 20.38 

E 19.02 

D 17.66 

C 16.30 

B 14.95 

A 13.59 

9 12.23 

8 10.87 

7 9.51 

6 8.15 

5 6.79 

4 5.43 

3 4.08 

2 2.72 

1 1.36 


Figure 6.16. Cycle- Averaged Unresolved Velocity Correlation 
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x r /c r 


Level uv(%) 

F 9.40 

E 7.40 

D 5.40 

C 3.39 

B 1.39 

A -0.62 

9 -2.62 

8 -4.62 

7 -6.63 

6 -8.63 

5 -10.63 

4 -12.64 

3 -14.64 

2 -16.64 

1 -18.65 


Figure 6.17. Cycle-Averaged Unresolved Velocity Correlation 
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( v - v ) p ( v - v ) 

V p V q =- ^ — x 100% (6.16) 

m 


where the subscripts p and q correspond to the axial and /or tangential velocity 
components. The periodic velocity correlations are then cycle averaged using 
equation (6.8). All three show high levels of periodic unsteadiness at the rotor 
leading edge, near the pressure and suction surfaces and in the rotor wake. The 
high level of periodic unsteadiness at the leading edge is a result of the large 
potential effect of the rotor blade on the flow field, while the high levels near the 
blade surfaces is a result of the large velocity change across the rotor pitch. The 
rotor wake periodic unsteadiness decays fairly rapidly, decaying to negligible 
values within one half chord downstream of the rotor. 


6.3 Rotor Time Resolved Flow Field Including Nozzle Wake Propagation 
through the Rotor 

To understand the propagation of the nozzle wake through the rotor, one 
needs to look at the flow field in the rotor at different relative positions of the rotor 
with respect to the nozzle. This will show of the chopping the nozzle wake and its 
transport through the rotor, and is accomplished by looking at the rotor flow field 
at the six different nozzle/rotor blade locations measured (as shown in Figure 6.1). 

There are several criteria that can be used to identify the nozzle wake in the 
rotor passage. Compared to the free stream flow field outside of the wake, the wake 
has a velocity defect, higher unresolved unsteadiness, a variation in flow angle 
across the nozzle wake, and higher shear stress. These criteria will be used in this 
section to determine the presence and propagation of the nozzle wake through the 
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rotor passage. 

The level of the interaction between the nozzle and rotor flow field can be 
determined by examining the ratio of the time it takes the rotor to traverse one 
nozzle pitch versus the time it takes for fluid particles to travel through the rotor, 
blade passage. This is called the reduced frequency and it is given by 


n = 


C,/V. 

s/u m 


(6. 17) 


where Q, is the reduced frequency, C r is the rotor axial chord, V x is the axial 

velocity at theinlet to the rotor and S is the nozzle pitch. This ratio determines the 
number of nozzle wakes in each rotor passage at any instant in time. For the 
turbine in this investigation, the reduced frequency is 1.5, which means that there 
should be one and a half nozzle wakes in each rotor passage for each nozzle/rotor 
blade location, which is demonstrated in the following paragraphs. 


6.3.1 Relative Total Unresolved Unsteadiness 


Figures 6.18a through 6.18f show the relative total unresolved unsteadiness 
at the six different nozzle/rotor locations. They represent six different "snapshots" 
of the rotor flow field, and since these six different locations are equally spaced 
over one nozzle pitch, they can be viewed sequentially from nozzle/rotor location 
one to location six and then back to one again. Examining position 2 first,a region 
of increased unresolved unsteadiness upstream of the rotor leading edge can be 
seen as compared to the cycle averaged relative total unsteadiness upstream of the 
rotor presented in Figure 6.5. This is the nozzle wake. Moving to position 3, the 
nozzle wake enters the rotor passage, and is subsequently chopped into individual 
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segments by the rotor blades. The nozzle wake is bowed because the convection 
velocity at midpitch is higher than at the rotor leading edge. These individual 
segments of the nozzle wake can now move independently of each other resulting 
in a mismatch between segments that were originally part of the same nozzle wake 
by the time they reach the rotor exit. There is also increased unresolved 
unsteadiness near the pressure surface of the rotor blade, just downstream of the 
point where the nozzle wake interacts with the rotor pressure surface. This increase 
in unresolved unsteadiness is a result of the interaction of the nozzle wake with the 
rotor pressure surface boundary layer which, along with the concave curvature 
effects discussed earlier, destabilizes the flow. 

At position 4, the nozzle wake is becoming distorted as it travels through 
the rotor passage, with the region of the nozzle wake near the rotor suction surface 
moving faster than the region near the pressure side. This distortion of the nozzle 
wake is due to the large differential in the convection velocity between the pressure 
and suction surface, especially at the leading edge. Moving to positions 5 and 6, the 
distortion of the nozzle wake is continuing with the nozzle wake turning clockwise 
in the rotor passage. Continuing onto positions 1 and 2, the nozzle wake has turned 
more than 30 degrees from its orientation at the rotor leading edge and it is now 
parallel to the rotor pressure surface. At position 3, the nozzle wake is stretched 
along the rotor pressure surface. Thus, it can not be easily identified in the rotor 
passage since it is close to the rotor pressure surface where is cannot be measured 
by the LDV. This is in contrast to other measurements of the nozzle wake in the 
rotor passage (Hodson, 1984 and Binder et al., 1985) and computations (Hodson, 
1985 and Korakianitis, 1992) who show that at the rotor trailing edge, the nozzle 
wake stills spans the rotor passage from the pressure to suction surfaces. The 
reason for this can be found by looking at the blade surface velocity distributions. 
In both Hodson's and Binder’s turbine rotors the difference between the pressure 
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and suction surface velocities is much less than the difference between the pressure 
and suction surface velocities for the AFTRF rotor. (Korakianitis (1992) calculates 
the flow in Binder's turbine rotor, while Hodson (1985) calculates the flow in his 
own turbine rotor.) At midchord, the suction surface to pressure surface velocity 
ratio of the AFTRF turbine rotor is twice that of Hodson's and Binder's turbine 
rotors suction to pressure surface velocity ratios. This large difference between the 
pressure and suction surface velocities in the AFTRF turbine rotor causes the 
region of the nozzle wake near the rotor suction surface to travel much more 
rapidly than the region of the nozzle wake near the pressure surface. Thus by the 
time the region nozzle wake near the rotor suction surface has reached the rotor 
trailing edge, it has rotated so that it hugs the pressure surface. In Hodson and 
Binder's turbine rotors, the difference between the pressure surface and suction 
surface velocities is not as great, thus the nozzle wake does not turn as much and 
by the time it reaches the rotor trailing edge it still spans the rotor passage from the 
pressure to suction surface. 

Using an averaged convection velocity along the AFTRF rotor pressure 
surface between X/C=0.50 and X/C=1.0 (where the nozzle wake is located at 
position 2) to calculate the distance the nozzle wake moves between each position, 
the nozzle wake should be completely inside the rotor wake after six nozzle/rotor 
locations. Thus six nozzle/rotor locations from position 2 is position 2. 

This is confirmed by examining the flow field downstream of the rotor. It 
can be seen that the flow field is not the same at every nozzle/rotor position, but 
changes from one position to another. Sharma et al. (1985) also noticed this 
phenomena downstream of their rotor. Measuring the flow field at 10 % axial chord 
downstream of a rotor with a three sensor hot wire, they showed that there are two 
distinct flow fields downstream of the rotor which they called the minimum and 
maximum interaction between the nozzle and rotor wakes. The maximum 
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interaction occurs when the nozzle wake is directly inside of the rotor wake while 
the region outside the wake shows low total unsteadiness. The minimum 
interaction occurs when the upstream nozzle wakes are in between the rotor wakes, 
which is shown by the high total unsteadiness in this region. A similar feature can 
be seen in the AFTRF rotor flow field just downstream of the rotor (from the 
trailing edge to one half chord downstream). The maximum interaction occurs at 
position 2 with high relative total unsteadiness in the rotor wake (maximum 
relative total unsteadiness of 16% at X r /C r = 1.12) and low total unsteadiness in 
the region between the wakes (relative total unsteadiness of 4% at X r /C r = 1.12). 
This low relative total unsteadiness region occupies more than 60% of the rotor 
pitch at 10% axial chord downstream of the rotor trailing edge. 

The minimum interaction occurs at position 5 with the low relative total 
unsteadiness region between the wakes (Tu r = 4%) occupying only 15% of the 
rotor pitch at X r /C r =1 . 10. At this location the maximum relative total 

unsteadiness in the rotor wake is 12% as compared to 16% at the maximum 
interaction. This is because the nozzle wake is located in between the rotor wakes 
at position 5, causing higher relative total unsteadiness in between the rotor wakes 
while at position 2 the nozzle wake is located in the rotor wake, thus causing 
higher relative total unsteadiness in the rotor wake. 

Blade to blade profiles of total unresolved unsteadiness for nozzle/rotor 
positions 4 and 1 are presented in Figures 6.19a and 6.19b, respectively, at five 
typical axial locations in the rotor passage. The total unresolved unsteadiness is 
normalized by U m here, so that the magnitude of the unsteadiness at different axial 
locations can be compared. While near the leading edge (X r /C r =0.01 and 0.04) the 
unresolved unsteadiness increases from 6% outside the nozzle wake (S r =1.2) to 
10% inside the nozzle wake (S r =1.0), near the trailing edge (X r /C r = 0.706) the 
unresolved unsteadiness increases from 2.5% outside the nozzle wake (S r =0.7) to 
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Figure 6.19. Blade to Blade Profiles of Total Unresolved Unsteadiness (Tu,) 
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10% inside the nozzle wake (S r = 0.9). This demonstrates that the unresolved 
unsteadiness in the nozzle wake does not decay significantly in the rotor passage. 
This is discussed in more detail in Chapter 7. 

6.3.2 Relative and Absolute Total Velocity 

Figures 6.20a through 6.20f show the relative total velocity in the rotor at 
the six nozzle/rotor positions. The location of the nozzle wake in the rotor can be 
identified by the region that has a lower total velocity than the cycle averaged total 
velocity shown in Figure 6.3. Using this criteria the nozzle wake can be identified 
at nozzle/rotor position 2, upstream of the rotor leading edge at the same location 
as the increased relative total unsteadiness shown in Figure 6.18. Examining the 
successive nozzle/rotor locations, this velocity defect is located at the same 
location as the higher relative total unsteadiness. This gives added confirmation 
that the nozzle wake is located at these positions. Downstream of the rotor trailing 
edge, one can see the flow field is not the same at every nozzle/rotor position, but 
changes from one position to another, similar to the total unsteadiness plots. At the 
maximum interaction position (position 2), the free stream flow field outside of the 
rotor wakes has a uniform velocity over almost 70% of the rotor pitch at 
X r /C r =1.10. At this position the nozzle wake is located in the rotor wake (in the 

region from the rotor trailing edge to one half chord downstream). On the other 
hand, at the minimum interaction position (position 5), the free stream flow field at 
X r /C r =1.10 covers only 28% of the rotor pitch. At this position the nozzle wake 

is in between the rotor wakes. Also, the free stream velocity is lower at the 
minimum interaction position (from the trailing edge to X r /C r = 1.50) than at the 

maximum interaction position thus providing added evidence that the nozzle wake 
is located between the rotor wake at the minimum interaction position. 
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Figure 6.20. Relative Total Velocity (w/U m )at the Six Nozzle/Rotor Locations 
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Figure 6.20 (Cont.). Relative Total Velocity (w/U m )at the Six Nozzle/Rotor 
Positions 
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Blade-to-blade profiles of absolute total velocity for nozzle/rotor positions 4 
and 1 are presented in Figures 6.21a and 6.21b at five different axial locations in 
the rotor passage. The profiles are plotted over two rotor blade pitches. These axial 
positions represent the location of the nozzle wake in the rotor passage. The nozzle 
wake is just downstream of the rotor leading edge at nozzle/rotor position 4 (see 
Figure 6.18), while it spans the rotor passage from X r /C r = 0.2 to 0.8 at position 1. 

The velocity defect and wake width varies over the rotor pitch at each location. 
Both the wake width and the velocity defect grow larger as one moves from the 
rotor pressure to the suction side, which results from the nozzle wake acting like a 
negative jet and drawing the low momentum fluid from the rotor pressure side 
boundary layer to the suction side. The addition of this low momentum boundary 
layer fluid to the nozzle wake cause its defect and width to increase. 

6.3.3 Relative Flow Angle 

The relative flow angle is presented in Figures 6.22a through 6.22f. Starting 
with nozzle/rotor position 2, an overturning of the flow angle is associated with the 
regions that correspond to the nozzle wake as identified by the relative total 
unsteadiness and velocity. Flow overturning in the rotor is identified by a lower 
relative flow angle (or higher negative angle), since a lower angle indicate more 
turning of the flow in the rotor. This overturning in the nozzle wake results from 
the slip velocity toward the rotor suction surface shown in Figure 6.2. This causes 
the nozzle wake to act like a negative jet drawing fluid from the rotor pressure side 
to the suction side. 

Figure 6.23 presents the relative flow angle in the nozzle wake upstream of 
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Figure 6.21. Blade to Blade Profiles of Absolute Total Velocity (v/U m ) 
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the rotor, that was measured by the five hole probe. Close to the nozzle trailing 
edge, there is a large change in angle across the wake, which decreases as the wake 
travels downstream. At the closest measuring location to the rotor (X r /C r = -0.08), 

overturning in the nozzle wake still exists, with the flow being overturned by 4 
degrees at this location. When the nozzle wake is at the rotor leading edge, the 
overturning in the nozzle wake causes the flow to enter the rotor at a negative 
incidence. Both Hodson and Dominy (1987) and Tremblay, Sjolander and 
Moustapha (1990) observed that at negative incidence a large separation bubble 
occurs on the pressure side of the blade. This separation bubble causes the losses to 
increase in the blade passage thus causing a decrease in performance. Thus the 
nozzle wake passing through the rotor blade is detrimental to the turbine 
performance. 

The blade to blade profiles of relative flow angles in the nozzle wake at 
positions 4 and 1 are presented in Figures 6.24a and 6.24b, respectively. The 
overturning in the nozzle wake, which results from the nozzle wake acting as a 
negative jet, is seen clearly. The overturning in the nozzle wake at position 4 is 
greater than the overturning in the nozzle wake at position 1 . This results from the 
decrease in the slip velocity (the negative jet) in the nozzle wake as the nozzle 
wake travels through the rotor passage. 

6.3.4 Unresolved Velocity Cross Correlations 

The axial-tangential cross component of the unresolved velocity correlation 
for each nozzle/rotor position is shown in Figures 6.25a through 6.25f. The 
unresolved velocity cross correlation in the nozzle wake is higher than in the 
surrounding fluid. These higher velocity cross correlation regions are located in the 
same regions as where the nozzle wake was identified using the relative total 
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unsteadiness, relative total velocity and relative flow angle, thus providing added 
confidence for the existence the nozzle wake. The unresolved velocity cross 
correlation starts out being negative in the nozzle wake at position 2. (The nozzle 
wake is located upstream of the rotor blade at this position.) At position 3, the 
velocity cross correlation changes sign in the region of the nozzle wake near the 
pressure surface, with this region being positive while the region of the nozzle 
wake near the suction surface is negative. Continuing on to the next locations, the 
turning and the distortion of the nozzle wake can be seen. These figures show that 
as the nozzle wake travels downstream through the rotor passage it thins out near 
the rotor pressure surface and thickens near the rotor suction surface. This is due to 
the nozzle wake acting like a negative jet as discussed in the previous paragraph. 
Examining positions 1 and then 2, the nozzle wake is seen to have elongated and 
thinned out considerably. This is a result of two reasons. The first is an inviscid 
phenomena which was discussed by Smith (1966), and is due to the fact that 
vorticity must be conserved in the nozzle wake, thus as the wake length grows, the 
wake width must become smaller. The second is a result of the large variation in 
convective velocity across the rotor pitch. Since the velocity is faster along the 
rotor suction surface than near the pressure surface, the region of the wake near the 
suction surface moves faster than the region near the pressure surface. Thus the 
wake stretches, becoming narrower and longer. 

Downstream of the rotor trailing edge, one can see the flow field is not the 
same at every nozzle/rotor position, but changes from one position to another, 
similar to the total unsteadiness and relative velocity plots. For the region just 
downstream of the trailing edge at the maximum interaction (position 2), the 
unresolved velocity cross correlation is high in the near rotor wake (at X r /C r =1.04 

it is a maximum of -3.0) and is low (zero) in the free stream region outside of the 
rotor wake. This low unresolved velocity cross correlation region takes up over 
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60% of the rotor pitch at X r /C r =1.12. The minimum interaction position (position 

5) has a maximum unresolved velocity cross correlation that is lower in the rotor 
wake (maximum unresolved velocity cross correlation of -2.0 at X r /C r =1.04) and 

a higher unresolved velocity cross correlation in the region between the rotor wake 
(maximum unresolved velocity cross correlation of 0.40 at X r /C r =1.12) as 

compared to the maximum interaction position. The low unresolved velocity cross 
correlation (u'v'= 0) occupies only 25% of the blade pitch at X r /C r =1.1. This is 

because at the maximum interaction position, the nozzle wake is inside of the rotor 
wake, thus causing higher unresolved velocity cross correlations in the wake and 
outside of the rotor wake the unresolved velocity cross correlations are low, while 
at the minimum interaction position the nozzle wake is located in between the rotor 
wake, thus causing higher unresolved velocity cross correlations in the region 
between the rotor wake. 

Figures 6.26a and 6.26b show the blade to blade profiles of unresolved 
velocity cross correlation in the nozzle wake at positions 4 and 1, respectively. The 
increase in unresolved velocity cross correlation in the nozzle wake is easily 
identifiable. Outside the wake the cross correlation is zero, while inside the wake it 
has a value of 0.5% which is constant in the nozzle wake as it migrates from the 
rotor leading edge to the trailing edge. As with the total unsteadiness, this shows 
that the unresolved velocity cross correlation does not decay any significant 
amount in the rotor passage. This is discussed in more detail in Chapter 7. 

6.3.5 Axial and Tangential Components of Unresolved Velocity Correlations 

The axial and tangential components of the unresolved velocity correlations 
are given in Figures 6.27a through 6.27f and 6.28a through 6.28f, respectively. The 
axial and tangential components of the unresolved velocity correlations are defined 
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Position #2 




Figure 6.27. Unresolved Velocity Correlation (uV) at the Six Nozzle/Rotor 
Positions 







Figure 6.27 (Cont ). Unresolved Velocity Correlation (u'u') at the Six 
Nozzle/Rotor Positions 
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Figure 6.28 (Cont.). Unresolved Velocity Correlation (vv) at the Six 
Nozzle/Rotor Positions 
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as follows: 


(Y-vWv-v) 

VV = — x 1 00% 


p q 


Uf 


(6.16) 


where the subscripts p and q correspond to the axial and /or tangential velocity 
components. The nozzle wake and its propagation through the rotor passage can be 
clearly identified by the higher unresolved velocity that occurs in the nozzle wake 
as compared to the surrounding fluid. Starting from position 2, where the nozzle 
wake is upstream of the rotor, and continuing onto the successive positions, it can 
beseen that where the nozzle wake interacts with the rotor leading edge an increase 
in the magnitude of the velocity fluctuations occur for both the axial and tangential 
components of unresolved velocity correlations. 

Another interesting feature is the increase in the tangential component of 
the unresolved velocity correlation in the nozzle wake at rotor midpitch (shown in 
Figure 6.28) . At position 2 (where the nozzle wake is located upstream of the 
rotor), the maximum tangential component of the unresolved velocity correlation is 
1.75. Continuing onto positions 3 and 4 where the nozzle wake moves into the 
rotor passage, the maximum tangential component of the unresolved velocity is 2.0 
and 2.5, respectively. Binder et al. (1987) and Hathaway (1986) also notice this 
feature. While Hathaway does not provide an explanation for this feature in his 
compressor, Binder attributes this increase in random velocity fluctuations within 
the chopped nozzle wake segments at the turbine rotor leading edge midpitch to the 
cutting of the secondary flow vortices by the rotor resulting in the vortices breaking 
down. While this is a plausible explanation for the increase in random velocity 
fluctuations in Binder's turbine, it is not a correct explanation for this phenomena 
in the present turbine. This is because the nozzle secondary flow vortices were 
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located at midspan in Binder's axial flow turbine, which was his measurement 
location, while for the present turbine there are no nozzle secondary flow vortices 
at midspan (the measurement location). The increase in fluctuating velocity can be 
explained by looking at the nozzle wake as a negative jet drawing fluid from the 
rotor pressure to the rotor suction surface as discussed earlier. This movement of 
low momentum fluid toward the suction surface causes the high momentum fluid 
in the free stream to move in an opposite direction to replace the migrated wake 
fluid. This interaction induces two counter-rotating vortices, one on either side of 
the nozzle wake. The generation of these vortices causes an increase in the velocity 
fluctuations. Moving on from position 4 to position 5 the maximum tangential 
component of unresolved velocity correlation in the nozzle wake decreases and this 
decrease continues as the nozzle wake travels downstream in the rotor until 
position 2 where the maximum tangential component of the unresolved velocity in 
the nozzle wake is 0.50%. 

The flow field downstream of the rotor trailing edge shows different 
features at each nozzle/rotor location due to the nozzle wake interaction with the 
rotor wakes. As with the previous properties, the maximum interaction position 
occurs at nozzle/rotor position 2 (for the near wake region). The unresolved 
unsteadiness is low between the rotor wakes and high in the rotor wakes since the 
nozzle wake is located inside of the rotor wake at this position. The minimum 
interaction for the near wake region occurs at nozzle/rotor position 5. Compared to 
the maximum interaction position the unresolved unsteadiness is lower in the rotor 
wake and higher in between the rotor wakes, since the nozzle wake is located in the 
free stream region between the rotor wakes. 
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6.3.6 Periodic Velocity Correlations 

The periodic velocity correlations include the variation of velocity across 
the passage caused by potential and viscous effects. This is substantial even in the 
absence of the nozzle wake. This quantity is useful in evaluating the average 
passage equations used in industry for prelimary design (see Suder et al., 1987). 
The effect of the nozzle wake on rotor relative mean flow is also captured in this 
correlation. 

The axial and tangential components of the periodic velocity correlation are 

presented in Figures 6.29a through 6.29f and 6.30a through 6.30f, respectively. The 

periodic velocity correlations are about five times higher than the unresolved 

velocity correlations in the rotor. This is due to the large variation of the rotor flow 

field in the pitchwise direction which is a steady flow field phenomena in the 

relative frame and can also be seen in the cycle averaged periodic velocity 

correlation presented in Figures 6.15 and 6.16. One can also see that while the 

nozzle wake can be identified in the tangential component of periodic velocity 

correlation plots, it cannot in the axial component of periodic velocity correlation 

plots. The reason for this is that the axial velocity has a much larger pitchwise 

variation in the rotor blade passage than the tangential velocity, thus the nozzle 

wake is obscured in the plots. Comparing the tangential component of the periodic 

and unresolved velocity correlations in the nozzle wake, it is clear that their 

magnitudes are similar. This can be seen clearly at positions 4, 5 and 1 . Inside the 

nozzle wake at position 4, the maximum tangential component of unresolved 

velocity correlation ( v V ) is 2.5 while the maximum tangential component of 

periodic velocity correlation ( vv ) is 3.0. Inside the nozzle wake at position 5, 

(vV) = 2.25 and fiv) = 3.0. And at position 1, (vV) =1.0 and (Vt) = 

' /max V /max ' /max V /max 

1.0 inside the nozzle wake. Thus both the tangential component of the periodic and 
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unresolved velocity correlations have an equal impact on the nozzle wake 
propagation through the rotor. 

Figures 6.31a through 6.3 If show the axial tangential cross components of 
the periodic velocity correlations for all six nozzle/rotor positions. In contrast to the 
plots of unresolved velocity cross correlations shown in Figures 6.25a through 
6.25 f, the nozzle wake cannot be identified in all the periodic velocity plots. Since 
the periodic velocity cross correlations are an order of magnitude greater than the 
unresolved cross correlations (due to the large pitchwise variation in the flow 
field), the effect of the nozzle wake would be hard to discern in the regions where 
there are large values of periodic velocity cross correlations. When the nozzle wake 
is in the areas where the periodic cross correlations are small, the nozzle wake can 
be seen. At nozzle/rotor locations 4, 5, 6 and 1, an increase in periodic cross 
correlations occurs in the region where the previous plots identified the nozzle 
wake. The magnitude of the periodic cross correlations in the nozzle wake is 
similar to the magnitude of the unresolved cross correlations in the nozzle wake, 
also providing added proof that the periodic and unresolved velocities have an 
equal impact on the nozzle wake propagation through the rotor. 


6.3.7 Unsteady Velocity Vectors 


The unsteady velocity vector field at the six different nozzle/rotor locations 
are presented in Figures 6.32a through 6.32f. The axial and tangential components 
of the unsteady velocity vector is calculated first as follows: 


V* = v — v 


(6.18) 


W, 


= w e -w e 


(6.19) 
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Figure 6.32. Unsteady Velocity Vectors at the Six Nozzle/Rotor Positions 
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where k is the nozzie/rotor location, — stands for the ensemble averaged velocity 
at each nozzle/rotor location and ... represents the cycle averaged velocity. These 
two components are then combined to form to unsteady velocity vector using the 
equation: 


W . , = xV. +ew e (6.20) 

Since the cycle averaged velocities were only averaged over six positions, the 
unsteady velocity vectors should be looked at qualitatively, not quantitatively. 

At position 2, the nozzle wake was identified upstream of the rotor leading 
edge by the relative total unsteadiness and shear stress plots, among others. 
Examining this region upstream of the rotor leading edge at position 2, the 
unsteady velocity vectors are seen to be moving from the pressure to the suction 
surface, which is a result of the slip velocity in the nozzle wake, discussed earlier 
(see Figure 6.2). Proceeding onto positions 3, 4 and 5 the slip velocity in the nozzle 
wake continues, causing the nozzle wake to act as a negative jet drawing fluid from 
the pressure to the suction surface, as discussed earlier. The two counter-rotating 
vortices on either side of the nozzle wake can be seen also. As the nozzle wake 
propagates through the rotor passage the counter-rotating vortex that is upstream of 
the nozzle wakes decays faster than the counter-rotating vortex that is downstream 
of the nozzle wake. At nozzle/rotor position 6, the upstream vortex has decayed 
and is not visible anymore. The nozzle wake is parallel to the rotor pressure surface 
and the downstream vortex is still strong, drawing fluid upstream along the 
pressure surface into the nozzle wake. Moving onto position 1, the downstream 
vortex has decayed considerably and by position 2 it is not visible anymore. 
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6.4 Mass-Averaged Properties 

The mass-averaged properties are obtained by circumferentially mass- 
averaging each property over one rotor blade pitch as follows: 

s r 

J QVjde 

Q = \ ( 6 . 21 ) 

jv,rd9 

0 

where Q is any property, the superscript n represents mass-averaging, r is the 
radius (midspan), 0 is the tangential distance, V K is the local axial velocity and S r 
is the rotor pitch. The mass-averaged cycle averaged relative total velocity is 
presented in Figure 6.33a. The data measured with the LDV is also compared to 
the mass-averaged relative velocity derived from the five hole probe in the nozzle 
wake. The comparison between the data derived from the two measuring 
techniques is excellent, being within 1.0% of each other at X r /C r = -0.08. The 
velocity is fairly constant from the nozzle trailing edge (which occurs at X r /C r = - 
0.27) through the rotor passage to about X r /C r =0.40 after which it accelerates 

rapidly, with an almost linear velocity distribution, until trailing edge. Just 
downstream of the blade the velocity drops rapidly until X r /C r =1.05 after which is 

decreases gradually until 1 .0 chord downstream of the trailing edge. 

The mass-averaged relative velocity for each nozzle/rotor location is shown 
in Figure 6.33b. The position of the nozzle wake can be identified by a dip in 
velocity. Using this method, the propagation of the nozzle wake through the rotor 
passage is clearly shown in this figure. At position 2, the nozzle wake is located 
upstream of the rotor leading edge. Continuing on to position 3, the nozzle wake 
has traveled downstream into the rotor passage. At position 3, the nozzle wake has 
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Figure 6.33b. Mass-Averaged Relative Total Velocity 
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moved to X r /C r = 0.05, while at position 4, the nozzle wake has traveled to 
X r /C r = 0.10. The nozzle wake continues moving downstream, reaching 20% rotor 
axial chord at position 5 and 40% rotor axial chord at position 6. By the time it has 
reached position 1, the nozzle wake has been sheared and rotated so that it is 
parallel to the rotor pressure surface. Thus the velocity defect would be spread out 
over a large axial distance, and it cannot be detected at position 1 and 2. 

Comparing this figure to the contour plots previously shown, the location of the 
nozzle wake coincides for both the contour and mass-average plots. 

The mass-averaged cycle averaged absolute total velocity is shown in 
Figure 6.34. Upstream of the rotor, the comparison between the mass-averaged 
absolute velocity derived from the five hole probe and the mass-averaged cycle- 
averaged absolute velocity derived from the LDV is also excellent, being within 
1.0% of each other. The absolute velocity is constant from the nozzle trailing edge 
(at X r /C r = -0.27) to the rotor leading edge. In contrast to the relative total 

velocity, the absolute velocity drops through the rotor. This is because a turbine 
extracts kinetic energy from the flow field to produce work. 

Figure 6.35a presents the mass-averaged cycle average relative flow angles. 
The mass-averaged relative flow angles are calculated from the mass-averaged 
relative tangential velocity and axial velocity as follows: 
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= AT AN 


f ^ } 

1 W, ' 


V v„, 


( 6 . 22 ) 


The large turning of the flow (110 degrees) that is a characteristic of turbine rotors 
is seen in this figure. (Higher flow turning produces larger pressure and 
temperature drops, and thereby higher work output.) Upstream of the rotor, the 
comparison between the mass-averaged relative flow angle derived from the five 
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hole probe and the mass-averaged cycle-averaged relative flow angle obtained with 
the LDV is excellent being within 0.5% of each other. 

The mass-averaged relative flow angle for each nozzle/rotor position are 
shown in Figure 6.35b. The location of the nozzle wake can be identified by the 
lower relative flow angle occurring at that axial location. This results from the 
overturning of the flow in the wake discussed earlier. Using this criteria the nozzle 
is seen upstream of the rotor leading edge at position 2. Continuing onto the next 
position, the nozzle wake is seen move through the passage until position 1, where 
it cannot be identified any longer. This is because the nozzle wake has rotated until 
it is parallel to the rotor pressure surface as discussed earlier. 

Figure 6.36a presents the mass-averaged cycle-averaged axial velocity in 
the rotor. Just as with the previously discussed mass-averaged values the 
comparison between the mass-averaged axial velocity derived from the five hole 
probe and the LDV is excellent upstream of the rotor. Just upstream of the rotor 
leading edge there is an increase in axial velocity. This is due to the potential effect 
of the rotor blade. Traveling downstream of the leading edge, the axial velocity 
increases steadily until 60% rotor axial chord, which is caused by the decreasing 
area of the rotor passage. The maximum velocity is at X r /C r = 0.60 which is just 
upstream of the blade throat (the blade throat is at X r /C r = 0.65). Downstream of 
65% rotor axial chord until the trailing edge, the axial velocity decreases sharply 
due to the increase in rotor passage area. The axial velocity then level off to a fairly 
constant velocity downstream of the rotor trailing edge. 

Examining the mass-averaged axial velocities for each nozzle/rotor position 
(presented in Figure 6.36b), the nozzle wake can be distinguished clearly by the 
velocity defect associated with it (as compared to the mass-averaged cycle- 
averaged axial velocity shown in Figure 6.36a). Using this method, the location of 
the nozzle wake in the rotor passage is the same as in the previous mass-averaged 
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Figure 6.35b. Mass-Averaged Relative Flow Angle 
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plots and the trajectory of the nozzle wake can be identified. While inside the rotor 
passage, the trajectory of the nozzle wake can be clearly seen, this is not the case 
downstream of the rotor. The nozzle wake trajectory downstream of the rotor is 
much more difficult to interpret for three reasons. The first is that the absolute flow 
angle changes direction (see Figure 6.37), which would cause the nozzle wake to 
move in the opposite direction downstream of the rotor as compared to inside the 
rotor passage (since the nozzle wake follows the flow direction), and the second is 
the absolute flow velocity decreases. The third reason is that in some nozzle/rotor 
positions, the nozzle wake is in the rotor wake freestream and decays slower 
(position 5), and in other locations (position 2), the nozzle wake is inside the rotor 
wake and decays faster. Despite this, the nozzle wake can be identified in the flow 
field downstream of the rotor at several locations. At nozzle/rotor position 5, there 
is an axial velocity defect at X r /C r = 1.2, which probably corresponds to the 
nozzle wake. At position 4, there is a defect at X r /C r = 1.6 and at position 3, there 
is a defect at X r /C r = 2.0, both of which could correspond to the nozzle wake. 

The mass-averaged cycle-averaged absolute flow angle in the rotor is shown 
in Figure 6.37.The flow angle continuously decreases from the rotor leading edge 
to the trailing edge, changing sign just upstream of the rotor trailing edge (at 
X r /C r = 0.8). Downstream of the trailing edge the absolute flow angle remains 

constant. Figure 6.37 also shows the comparison between the mass-averaged 
absolute flow angle derived from the five hole probe and the mass-averaged cycle- 
averaged absolute flow angle obtained with the LDV upstream of the rotor. As 
with the previously discussed mass-averaged properties, the comparison is 
excellent being with 0.3% of each other at X r /C r = -0.08. 

The mass-averaged cycle-averaged relative total unresolved unsteadiness is 





247 


shown in Figure 6.38a. Starting from upstream of the rotor and moving 
downstream, the unresolved unsteadiness decreases as the rotor leading edge is 
approached. This is caused by the nozzle wake decay. At the leading edge there is a 
sharp increase in unresolved unsteadiness as discussed earlier. The total unresolved 
unsteadiness decreases as one continues to travel downstream, since the relative 
velocity increases. Downstream of the trailing edge there is a sharp increase 
between the tr aili ng edge and 10% axial chord downstream. While the unsteadiness 
downstream of the rotor includes both nozzle wake and the rotor wake 
unsteadiness, this increase is only caused by the unsteadiness in the rotor wake. 
Continuing on downstream, the rate of increase in unresolved unsteadiness 
decreases until close to one chord downstream where the magnitude of relative 
total unsteadiness levels off. 

Figure 6.38b presents the mass-averaged relative total unresolved 
unsteadiness for all six nozzle/rotor locations. Using the increase in total 
unresolved unsteadiness in the rotor blade passage to identify the nozzle wake, the 
position of the wake is at the same locations where the previous mass-average plots 
found it to be at. Downstream of the blade the nozzle wakes are harder to identify, 
since the overall level of unresolved unsteadiness is higher than in the blade 
passage and thus the nozzle wakes are not as noticeable. But looking closely at the 
plot, there are regions of higher unsteadiness that could be associated with the 
nozzle wake. Starting with position 5, there is a high unresolved unsteadiness 
region just downstream of the trailing edge (at X r /C r = 1.10). At position 4, there 
is a region at X r /C r = 1.30, while at position 3, the peak unsteadiness occurs at 
X r /C r = 1.50. At both nozzle/rotor positions 2 and 1, there are high unresolved 
unsteadiness regions (at X r /C r = 1.8 and 2.0, respectively), that possibly could 
correspond to the nozzle wake. While the regions of high unresolved unsteadiness 
matched the axial position of the other mass-averaged properties in the rotor blade, 





249 





A O 






x x x 


*«&< 


XXXx xxxx)^ XXxxxx xxx 

^ + 

+ + + 

+++ ^ . .Jr 


++ 

++ + 


*++ ++ ++ +++ 


A + 

A 

£ A 

/ A 

A A 

A A A A A a 

% 


A A A A AAA A 


0 


q 00©°000 


°°°°°oo 0 'r‘° 


00 


oo 


o© 00 


a 


+ 


A 


© 


❖ 


X 


+ 


A 


© 


A 


X 


+ 


A 


O 


Pr - 


iS 

cl 


'%&m aa an a * aaaaa 'n 


Dnnirn^DBDGjB 0130 


□ □ ° 


JL 


-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

Xr/Cr 


D □ 


i I i i lJ 

1.8 2.0 


^ Position 
x Position 
+ Position 
A Position #3 
© Position # 2 
E Position #1 
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the region of high unsteadiness do not occur at the same location as the low axial 
velocity regions downstream of the rotor blade. This is because, although in the 
near wake the maximum unsteadiness occurs at the wake center (position where 
velocity defect is largest), far downstream of the nozzle the maximum turbulence 
intensity occurs away from the center due to the spread of the wake. Hah and 
Lakshminarayana (1982) also noticed this feature in the wake of an isolated airfoil. 

The mass-averaged cycle averaged unresolved velocity cross-correlation 

• • • 

(u'v') is presented in Figure 6.39a. The unresolved velocity cross correlation is 
zero upstream of the rotor leading edge. There is a sharp increase in unresolved 
velocity cross correlation at the leading edge due to the high unresolved velocity 
cross correlation in the stagnation region. The high levels of mass-averaged 
unresolved velocity correlation at midchord results from the high unresolved 
velocity cross correlation levels in the pressure surface boundary layer. The sharp 
increase in negative mass-averaged unresolved velocity cross correlations just 
downstream of the trailing edge is because high levels of negative unresolved 
velocity cross correlations in the rotor wake. Moving downstream of the trailing 

edge the mass-averaged unresolved velocity cross correlation decays very rapidly 
to zero at X r /C r = 1.25. This is mainly caused by the rapid decay of the rotor 

wake. The mass-averaged unresolved velocity cross correlation for each 
nozzle/rotor location is shown in Figure 6.39b. The regions of increased unresolved 
velocity cross correlation correspond to the nozzle wake location. These regions 
corresponds with the regions of increase total relative unsteadiness, thus giving 
added confirmation of the nozzle wake trajectory, especially downstream of the 


rotor. 
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Figure 6.39b. Mass-Averaged Unresolved Velocity Correlation 



253 


6.5 Nozzle Relative Flow Field (Time-Averaged Properties) 


The time-averaged properties at each stationary measuring location are used 
to discern the wake pattern in the absolute frame without any reference to the rotor 
blade positions. The time-averaged values Q are defined as 


Q = 


1 IN W 

— Sq 

NW" 


(6.23) 


where Q is any ensemble averaged property and NW is the total number of 
measurement windows. The time-averaged velocity is a time average of all 
measurements acquired at a fixed point in space. This is equivalent to equation 6.6. 

Figure 6.40 are the time-averaged total unresolved unsteadiness (normalized 
by U m to see the nozzle wake more clearly). The data represents two nozzle pitches 
(the data from one pitch was doubled). The rotor blade is not visible due to the 
integration method. The nozzle wake, identified by increased unresolved 
unsteadiness, moves on a curved path through the rotor blade, changing direction 
from the inlet of the rotor to the exit. It seems to follows the path of the absolute 
flow. The absolute flow path can be seen more clearly by looking at the time- 
averaged absolute velocity vectors shown in Figure 6.41. By comparing the path of 
the nozzle wake with that of the absolute flow, the nozzle wake is seen to follow 
the absolute flow. 

The time-averaged axial velocity is presented in Figure 6.42. The regions of 
low velocity correspond to the nozzle wake. By comparing the path of the low 
velocity regions with that of the high unresolved unsteadiness, the paths coincide in 
the rotor blade, while downstream of the rotor blade they deviated from one 
another. This is due to wake spreading which was discussed earlier. The time- 
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Figure 6.40. Time-Averaged Total Unresolved Unsteadiness lTu t 
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Figure 6.41. Time- Averaged Absolute Total Velocity Vectors 
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averaged unresolved velocity cross correlation is shown in Figure 6.43. The 
increased unresolved velocity cross correlation resulting from the nozzle wake can 
be seen. The path of the nozzle wake is clearly marked out downstream of the 
rotor, while in the rotor some of the high unresolved velocity cross correlation 
region result from the increased unresolved velocity cross correlation at the 
pressure surface, thus making the determination of the nozzle wake path more 
difficult. 
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CHAPTER 7 

ROTOR FLOW FIELD WAKE CHARACTERISTICS 

An understanding of the steady and unsteady characteristics of turbine 
nozzle and rotor wakes are important for the efficient design of axial flow 
turbomachinery. A major cause of noise and vibration characteristics of 
turbomachinery is caused by wakes. Turbine wakes represent a source of loss in 
efficiency, since the mixing of the wakes with the ffeestream dissipates energy. 

The three-dimensional characteristics of the wake, the decay characteristics and the 
path that it follows is important in the design of the following blade rows. This 
information is essential for both the prediction of the aerodynamic and mechanical 
performance of a turbine and for building quieter turbomachines. 

7.1 Nozzle Wake Characteristics in Rotor Passage 

The nozzle wake absolute velocity profiles upstream of the rotor blade, in 
the rotor blade passage and downstream of the rotor blade are presented in Figures 
7.1 and 7.2a, and b respectively. (The bar over the velocity and angle notation in 
the figures in Chapter 7 represents ensemble averaging.) The measurement 
upstream of the rotor were acquired with the five hole probe which was discussed 
in Chapter 5. For the measurements in and downstream of the rotor, each axial 
location represents the axial location for each nozzle wake segment where the 
maximum velocity defect occurs. (A nozzle wake segment is defined as follows. 

As the nozzle wake enters the rotor passage, it is divided in to individual segments 
the rotor leading edge and these segments subsequently propagate through the rotor 
passage independent of each other.) The data is displayed over two rotor blade 
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Figure 7.1. Nozzle Wake Total Velocity (v/U m ) Profiles Upstream of Rotor (Five 
Hole Probe Measurements) 
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Figure 7.2 Nozzle Wake Total Velocity (v/U m ) Profiles in Rotor Passage 
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pitches (the data from one rotor pitch is doubled) to aid interpretation. 

Upstream of the rotor, the nozzle wake decays very rapidly as the rotor 
leading edge is approached. As the nozzle wake moves through the rotor passage, 
the absolute magnitude of its defect grows and then decreases. Downstream of the 
rotor trailing edge, the nozzle wake defect remains constant until X r /C r = 1.64. 

(The nozzle wake is the smaller of the two velocity defects in the flow field 
downstream of the rotor while the rotor wake is the larger defect.) The variation in 
nozzle wake defect is demonstrated more clearly in Figure 7.3, which shows the 
nozzle wake defect normalized by the local maximum velocity. (The local 
maximum velocity is the velocity that would exist at the location where the nozzle 
wake occurs if the nozzle wake did not exist.) This is consistent with the 
normalization of the nozzle wake velocity defect upstream of the rotor measured by 
the five hole probe, which is also presented in this figure. Examining this figure, 
one can see that the velocity defect decreases sharply from the nozzle trailing edge 
(located at X r /C r = -0.27) to the rotor leading edge. This rapid decay is due to 

pressure gradients, high turbulence intensities and wake centerline curvature along 
with the effect of the rotor. The relative motion between the rotor and the nozzle 
causes periodic variations in the potential flow field around the blades which 
causes the wake to decay faster. A more detailed discussion of the nozzle wake 
decay between the nozzle and rotor is given in Chapter 5 of this thesis. 

The nozzle wake then enters the rotor passage where the velocity defect 
increases from the rotor leading edge up to X r /C r = 0.30 and decreases after that 

until the trailing edge. Downstream of the trailing edge the velocity defect is 
constant. Matsuuchi and Adachi (1983) have observed a similar feature in their 
axial flow fan. The maximum velocity defect for each upstream stator wake 
segment in their axial flow fan rotor is presented also in Figure 7.3. As the stator 
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Figure 7.3. Variation of Peak Nozzle Wake Velocity Defect in Rotor Passage 
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wake travels through the rotor its velocity defect increases. This is because, 
according to Hill et al. (1963), a positive (adverse) pressure gradient slows the 
wake decay and if the pressure gradient is large enough the wake decay will be 
stopped completely and the wake will grow in size. This is what is happening here, 
since there is a positive pressure gradient in a compressor rotor. 

The reasons for the velocity defect increase and decrease in the turbine rotor 
are more complex. A possible explanation can be given by looking at the mass 
averaged cycle-averaged relative velocity, presented in Figure 6.33a. The relative 
velocity slightly decreases from the leading edge until X r /C r = 0.30, increases 
sharply from X r /C r = 0.40 until the trailing edge and then is constant downstream 

of the trailing edge. The change in wake defect seems to correspond to the relative 
velocity change, with the defect increasing until X r /C r = 0.30 and then decreasing 
downstream of X r /C r = 0.40. Downstream of the rotor trailing edge, the wake 
defect is constant. As discussed earlier, Raj and Lakshminarayana (1973) found 
that the wake decay rate corresponds to the variation in wake edge velocity. In the 
region where the relative velocity decreases, an adverse pressure gradient occurs, 
thus causing the nozzle wake velocity defect to increase. Between X r /C r =0.40 and 

the trailing edge, the relative velocity increases, resulting in a favorable pressure 
gradient which causes the nozzle wake to decay. 

Blade to blade profiles of total unresolved unsteadiness (defined by eqn. 
6.12) at the axial location where the maximum total unresolved unsteadiness 
occurs in each nozzle wake segment are shown in Figure 7.4. The total unresolved 
unsteadiness is normalized by U m here in order to compare the fluctuating velocity 
at different axial locations, since the velocity outside of the nozzle wake changes 
significantly in the rotor. One interesting feature is that the peak magnitude of total 
unresolved unsteadiness in the nozzle wake does not continuously decay as it 
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Figure 7.4. Nozzle Wake Total Unresolved Unsteadiness (Tu,) in Rotor Passage 
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travels through and downstream of the rotor passage. This is also established in 
Figure 7.5, which presents the peak total unresolved unsteadiness for each nozzle 
wake segment. The peak total unsteadiness increases slightly from the rotor leading 
edge to X r /C r = 0.30 after which it decreases dramatically. 

This increase and decrease in total unresolved unsteadiness is related to the 
nozzle wake velocity defect increase and decrease in the rotor passage shown in 
Figure 7.3. Where the nozzle wake velocity defect increases in the rotor passage 
(from the rotor leading edge to X r /C r = 0.30) , the total unresolved unsteadiness 
increases and where the nozzle wake decays (from X r /C r =0.40 to the rotor trailing 

edge) the nozzle wake total unresolved unsteadiness decays. An additional effect 
that could cause the decrease in total unresolved unsteadiness from X r /C r =0.40 to 

the rotor trailing edge can be seen by examining Figure 6.33a, which is the mass- 
averaged cycle-averaged relative velocity in the rotor. The relative velocity 
increases sharply downstream of X r /C r = 0.30 and since the total unsteadiness 

decays in an accelerating flow field, the maximum total unsteadiness in each 
nozzle wake segment decreases also. Downstream of the rotor, the nozzle wake 
total unsteadiness increases in magnitude until at one half chord downstream of the 
rotor it is at the same value as it is at the rotor inlet. Sharma et al. (1985) also 
notice in their axial flow turbine that the nozzle wake total unsteadiness 
downstream of the rotor is the same order of magnitude as it is upstream of the 
rotor. Another interesting feature of this flow field can be seen downstream of the 
rotor in Figure 7.4. Close to the trailing edge, the total unsteadiness of the rotor 
wake is much higher than the nozzle wake. As the wakes travel downstream, the 
rotor wake total unsteadiness decreases while the nozzle wake total unsteadiness 
increases until at one half chord downstream of the rotor the total unsteadiness for 
both the nozzle and rotor wakes are equal in magnitude. 
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Blade to blade profiles of unresolved velocity cross correlation (defined by 
eqn. 6.14) at the axial location where the maximum unresolved velocity cross 
correlation occurs in each nozzle wake segment are presented in Figures 7.6a and 
b. As the nozzle wake travels through the rotor passage, the velocity cross 
correlation grows until midchord, after which is decreases sharply. This is clear 
from Figure 7.7, which is the variation of peak unresolved velocity cross 
correlation in each nozzle wake segment. The unresolved velocity cross correlation 
increases as the nozzle wake travels through the rotor passage until midchord, after 
which is decreases. This decrease occurs at the same location as the decrease in 
total unresolved unsteadiness and they are probably related. 

7.2 Rotor Wake Characteristics 

The rotor wake can be classified into three different categories, the trailing 
edge region, the near wake and the far wake. The trailing edge region is confined to 
the area just downstream of the trailing edge. The velocity defect is very large in 
the region. In the near wake region, the physical characteristics of the blade and the 
aerodynamic loading on the blade have a major impact on the development of the 
wake, causing the wake to be asymmetric. The wake defect is of the same order of 
magnitude as the mean velocity in this region. In the far wake the wake structure is 
almost symmetric and the physical characteristics and the aerodynamic loading 
have minor effects on the development of the wake. The velocity defect is also 
small in the far wake region. For the turbine rotor wake data presented in this 
thesis, the individual regions are defined as follows: 
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Figure 7.6. Nozzle Wake Unresolved Velocity Correlation (uv) in Rotor Passage 
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Trailing Edge Region: X r /C r < 1.014 

(Z/cosp 0 < 0.035) 

Near Wake Region: 1.014 < X r /C r <1.22 

(0.035 < Z/cosP 0 < 0.555) 

Far Wake Region: X r /C r > 1.22 

(Z/cosp 0 >0.555) 

where Z/cosP 0 is the streamwise distance downstream of the rotor (Z is the axial 
distance downstream of the rotor with Z = 0.0 at the trailing edge and P 0 is the 
rotor blade outlet angle). 

7.2.1 Cycle-Averaged Properties 

Figures 7.8a and 7.8b show the cycle-averaged relative velocity profiles at 
several axial locations in the rotor wake. (As in chapter 6, all velocities presented 
in chapter 7 are ensemble averaged.) The abscissa in the figure represents the 
tangential distance normalized by the rotor blade pitch, with Y equal to 0.0 at the 
wake center. The velocity gradients in the tangential direction are very large just 
downstream of the rotor trailing edge. This feature results from the development of 
the flow as it transitions from the blade boundary layer to the wake. Farther 
downstream the gradient becomes much smaller due to wake spreading and mixing 
with the freestream as well as interchange of momentum and energy on either side 
of the wake. Also, the wake profiles are asymmetrical about the wake center due to 
the differential growth of the boundary layer on the two surfaces of the blade. The 
suction side of the wake has a larger width than the pressure side wake width since 
the suction surface boundary layer at the blade trailing edge is larger than the 
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pressure surface boundary layer. As the wake moves downstream, the differences 
between the pressure side and the suction side of the wake diminish until X r /C r = 

1 .22, where the wake profile becomes symmetrical about the wake center. This 
also results from wake spreading and mixing as discussed earlier. 

The cycle-averaged total unresolved unsteadiness (normalized by U m to 
show the absolute level of the fluctuations and defined by equation 6.12) is 
presented in Figures 7.9a and 7.9b. The total unresolved unsteadiness is highest 
near the trailing edge and decreases further downstream. The high values near the 
trailing edge result from the vortex street shed from the trailing edge and the high 
production of turbulence in this region. The unsteadiness profiles are asymmetric 
about the wake center in the near wake becoming symmetric farther downstream. 
This results from the differential growth of the boundary layer on the pressure and 
suction surfaces of the blade. The unsteadiness profiles have a dip in unsteadiness 
at the wake center with higher values of unsteadiness on either side of the dip. This 
is to be expected since the total unsteadiness is zero on the blade surface with the 
maximum value occurring slightly away from the surface. This dip in total 
unsteadiness disappears due to wake spreading and mixing, as the wake travels 
downstream. The highest total unsteadiness is observed on the pressure side of the 
wake which demonstrates that the pressure surface boundary layer has a higher 
total unsteadiness than the suction surface boundary layer. This is due to the 
concave curvature effect on the pressure surface and the interaction of the nozzle 
wake with the pressure surface boundary layer, discussed earlier. This increases the 
unresolved unsteadiness which persists even in the wake. The maximum total 
unsteadiness occurs close to the wake center in the near wake. Farther downstream, 
the maximum total unsteadiness occurs away from the wake center due to the 
spread of the wake. 
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The cycle-averaged unresolved velocity cross correlation in the streamwise- 
normal coordinate system (defined by equation 6.15) is shown in Figures 7.10a and 
7.10b. The streamwise-normal coordinate system is used to aid physical 
interpretation. The unresolved velocity cross correlation is very small in the free 
stream and reaches a maximum value in the wake center. This is to be expected 
since the unresolved velocity cross correlation results from velocity and turbulence 
intensity gradients. The velocity cross correlation changes sign near the wake 
center due to the opposite velocity gradients on either side of the wake center, but 
close to the trailing edge the zero unresolved velocity cross correlation does not 
occur at the point of minimum velocity. Other researchers, such as 
Lakshminarayana and Reynolds (1979) also see this feature in their rotor wakes. 
Beyond X r /C r = 1.035, the maximum value of unresolved velocity cross 

correlation occurs at the point of minimum velocity. In the near wake, the 
unresolved velocity cross correlation profiles are asymmetric about the wake 
center, with the magnitude of the cross correlation on the pressure side being larger 
than the magnitude on the suction side. The difference in the unresolved velocity 
cross correlation between the pressure and suction side of the wake decreases as the 
wake travels downstream until X r /C r = 1.22, where the cross correlation profile 

becomes symmetric about the wake center. The unresolved velocity cross 
correlation decays rather slowly with significant values of this correlation still 
occurring one half chord downstream of the rotor. 

The relative flow angle in the cycle-averaged rotor wake is presented in 
Figure 7.1 1. Just downstream of the rotor trailing edge, the flow is overturned on 
the suction side of the wake and undertumed on the pressure side of the wake. This 
is due to the difference in flow direction on the pressure and suction surface of the 
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Figure 7.1 1. Rotor Wake Cycle-Averaged Relative Flow Angle ( (3 ) Profiles 



278 


rotor. This over and undertuming of the flow angle decreases rapidly from a 
maximum of 17 degree difference at the trailing edge to a three degree difference 
in flow angle at 10 % chord downstream. 

7.2.2 Cvcle-Averaged Wake Decay Properties 

A knowledge of the rate of decay for the rotor wake defect is necessary for 
an understanding of the rotor-stator interaction as discussed earlier. The decay of 
the velocity defect is influenced by the pressure gradient, turbulence intensity, 
curvature and viscous effects. The velocity defect decay for the cycle-averaged 
rotor wake is shown in Figure 7.12. The very rapid decay in the trailing edge 
region results from the high unresolved unsteadiness and possibly the three 
dimensional effects in this region. In the near and far wake regions, the velocity 
defect decay is less rapid. The rotor wake velocity defect is also compared to the 
decay of the AFTRF nozzle wake, and the decay of two linear rotor cascades, 
Sitaram and Govardhan's (1986) and Ho and Lakshminarayana’s (1994) who 
computed the flow field in Gregory-Smith’s rotor cascade, which has similar flow 
turning as the present turbine rotor. The AFTRF rotor velocity defect decays much 
slower than the velocity defects of the other blades. This is in contrast to a 
compressor rotor wake which decays faster than the wake of a compressor cascade. 
This reason for this is not understood as of yet. It is possible that the low 
momentum fluid from other spanwise locations is being convected into the 
midspan wake and thus causing the midspan wake to be deeper and to decay 
slower. Future measurements at midspan and other spanwise locations that include 
the radial velocity are needed to verify this. 
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Schlichting's (1979) analysis indicates that a plane wake defect decays far 
downstream according to 



(7.1) 


where V c is the velocity defect, is the maximum velocity in the free stream 
outside of the wake, and S is streamwise distance. Using this relationship to 
correlate the rotor wake velocity defect results in the following equation 



= 0.14 


( Z v 


COS p 


O J 


(7.2) 


where W c is the cycle-averaged velocity defect, is the cycle-averaged 
maximum relative velocity in the free stream outside of the wake, Z/cos p c is the 

streamwise distance downstream of the rotor (Z is the axial distance downstream of 
the rotor with Z = 0.0 at the trailing edge and P n is the rotor blade outlet angle). 
This correlation can be seen to match the data quite well over most the region 
downstream of the rotor in Figure 7.12. It is not very good just downstream of the 
rotor trailing edge. In this region, the flow is very complex and is dominated by 
flow separation and trailing edge vortices along with the high levels of total 
unsteadiness, thus the wake decays much faster than farther downstream. And since 
the correlation was developed for a far wake, it does not hold in the trailing edge 
region. The constant in eqn. 7.2, 0.14, probably depends on the aerodynamic 
properties of the blade such as blade loading and turbulence intensity. These 
properties need to be varied for a more general wake decay correlation to be 
derived. 

The variation of the cycle-averaged rotor semi-wake width with streamwise 
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distance is shown in Figure 7.13. Semi- wake width is defined as the width of the 
wake where the total relative velocity defect is half. The data from the Sitaram and 
Govardhan's (1986) rotor cascade is also presented there for comparison. Both 
wake width's increase rapidly just downstream of the trailing edge (in the trailing 
edge region) and then grow more gradually in the near and far wake regions. Very 
far downstream (beyond Z/cosP 0 = 1.5), the AFTRF rotor wake width grows at a 
much slower rate. Since the interchange of mass, energy and momentum is 
continuous on both side of the wake as the wake travels downstream, the growth 
should be continuous and this is seen in Figure 7.12. 

Schlichting's (1979) analysis indicates that the increase of rotor semi-wake 
width is proportional to streamwise distance as follows, 

LocS* (7.3) 


where S is streamwise distance and L is semi-wake width (the wake width at half 


the maximum velocity defect). Using equation (7.3), the variation of semi-wake 
width over most of the streamwise distance can be represented quite accurately by 
the equation 


L/S = 0.389(Z/cosf) o )5 
for 0.088 < (Z/cos p 0 ) < 2.5 


(7.4) 


The semi-wake width variation of Sitaram and Govardhan (1986) turbine rotor also 
matches this equation. Their turbine rotor has similar loading as the AFTRF 
turbine rotor (AFTRF rotor flow turning is 1 10 degrees and Sitaram and 
Govardhan's rotor flow turning is 120 degrees). Since wake width depends on the 
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Figure 7.13. Variation of Cycle- Averaged Semi-Wake Width with Streamwise 
Distance 
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aerodynamic properties of the blade and this should be a function of the blade 
loading, equation (7.4) probably would need to be modified to match the wake 
width data of a turbine rotor with different loading. 

An interesting feature in the variation of semi-wake width is observed 
between the trailing edge and near wake regions. The semi-wake width grows 
rapidly in the trailing edge region (from the trailing edge to Z/cosP Q = 0.061) and 

then decreases sharply between the trailing edge region and the near wake region 
(between Z/cosP 0 = 0.061 and 0.088). It then increases more gradually beyond 
Z/cosp 0 = 0.088. Reynolds et al. (1979) also see this feature in their compressor 
rotor wake width. This phenomenon is a characteristic of rotor wakes, since the 
two turbine cascade wake widths do not have this decrease in wake width. It is 
caused by the effect of three-dimensional flow (radial flows) on the wake which 
are not present in cascade wakes. The radial transport of mass, momentum and 
energy could be responsible for the small decrease in wake width between the 
trailing edge and near wake regions, although the wake defect is decaying steadily 
there. 

The decay of the maximum cycle-averaged relative total, axial and 
tangential unsteadiness is presented in Figure 7.14. The relative total unresolved 
unsteadiness is defined by equation (6.9) while the relative unresolved axial and 
tangential unsteadiness is defined as follows: 
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Figure 7.14. Decay of Rotor Wake Cycle- Averaged Maximum Relative Total 
Unresolved Unsteadinesses 
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The relative toted, axial and tangential unresolved unsteadiness is then cycle- 
averaged using equation (6.8). While the unsteadiness is very high in the trailing 
edge region, it decays rapidly as the wake travels downstream. In the far wake 
region, beyond Z/cosf3 0 = 1.50, the relative unsteadiness has decayed to its 

magnitude upstream of the rotor blade. Each of the relative unsteadinesses do not 
decay at the same rate but can be shown to decay at a rate given by 


(Tu) = 

V /max 


f Z 
^COsP 


V 


O ) 


(7.5) 


according to Raj and Lakshminarayana (1976), where K and n assume different 
values for each unsteadiness. For the AFTRF turbine rotor, the rate of decay of the 
maximum cycle-averaged relative total, axial and tangential unsteadinesses are 
given by the following expressions, based on the data: 
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The above expressions match the data quite well, as seen in Figure 7.14. The good 
match of the equation (7.5) with both Raj and Lakshminarayana's (1976) 
compressor rotor maximum unsteadiness and the AFTRF turbine rotor data shown 
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here gives confidence that it can be used as a general equation for the decay of 
Reynolds' stresses downstream of turbine and compressor rotors. Since the 
exponent n varies from -0.18 to -0.22 for the turbine rotor and -0.19 to -0.24 for the 
streamwise and normal components of relative unsteadiness in the compressor 
rotor (see Raj and Lakshminarayana, 1976), assuming the exponent is equal to -0.2 
will give a fairly good fit for the expression. The constant K probably depends on 
factors such as the blade drag coefficient, the inlet unsteadiness and the mean 
velocity. Including the effect of these factors into the equation could collapse the 
data into a universal curve and lead to good correlation. For this to be done, 
though, the unsteadiness in the rotor wake at other conditions, such as different 
blade loading, inlet unsteadiness and mean velocity needs to be measured since the 
turbine rotor data here represents the only turbine rotor unsteadiness data the 
author has knowledge of. 

The decay of the maximum cycle-averaged streamwise-normal unresolved 
velocity cross correlation is presented in Figure 7.15. The unresolved velocity 
cross-correlation is very high in the trailing edge region and it decays very rapidly 
as the wake travels downstream. The rate of decay seems to follow the same trend 
as the decay of the relative unsteadinesses and using an equation of the form of 
equation (9.5), the maximum cycle-averaged unresolved velocity cross correlation 
decay rate can be given by 
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uX 

V /max 


= 0.335 


cosf) 


\- 0.644 


o 7 


(7.9) 


This correlation agrees with the data quite well as can be seen in Figure 7.15. 
The decay of the maximum cycle-averaged periodic and unresolved 
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Figure 7.15. Decay of Rotor Wake Cycle- Averaged Ma ximum Unresolved 
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unsteady velocity axial tangential and cross correlations are shown in Figures 7.16, 
7.17 and 7.18. For all three correlations, the maximum periodic unsteady velocity 
correlations are much larger than the maximum unresolved correlations in the 
tr ailin g edge region. This is because the periodic velocity correlations are due to 
the periodic variation in velocity over the rotor blade pitch, which is large in the 

rotor wake. But the periodic correlations decay much faster than the unresolved 

• •• 

correlations so that the maximum periodic axial velocity correlation (uu) is less 

• • • 

than the maximum unresolved axial velocity correlation (u'u') downstream of 

Z/cos(3 0 =0.25, while the magnitude of the maximum periodic tangential velocity 

• • • 

correlation ( vv) becomes lower than the maximum unresolved tangential velocity 

• « • 

correlation (vV) downstream of Z/cos(3 o =0.75. While the maximum periodic 

• •• 

unsteady velocity cross correlation (uv) decays more rapidly than the maximum 

• •• 

unresolved velocity cross correlation (uV), the magnitude of the maximum 
periodic cross correlation never is lower than the magnitude of the maximum 
unresolved cross correlation. Both decay to negligible values in the far wake. The 
decay of the maximum cycle-averaged unresolved velocity correlation reflect the 
same trend as the relative unsteadiness and can be modeled using equation (7.5) as 
follows: 
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Figure 7.18. Decay of Rotor Wake Cycle- Averaged Maximum Unresolved 


and Periodic Velocity Correlation 
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The decay of the maximum periodic velocity correlations, on the other hand, only 
follow this trend in the near and far wake regions; in the trailing edge region the 
decay is linear. Thus in the trailing edge region the maximum cycle-averaged 
periodic velocity correlation decay rates conform to the following relationships: 
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(7.15) 


while in the near and far wake regions the maximum periodic velocity decay rate 
correspond to the equations: 
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where 0.17 < Z/cos(3 0 <2.5. The above correlations correspond the experimental 
data quite closely. They also give added confidence that equation (7.6) can be used 
to correctly model both the decay of each component of the Reynolds stress tensor 
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and the periodic velocity correlations, too. Since the periodic velocity correlations 
also correspond to the velocity defect, the above correlations could be used to 
model the total velocity defect decay, too. 

The cycle-averaged wake momentum thickness and shape factor variations 
downstream of the rotor are shown in Figures 7.19 and 7.20. The momentum 
thickness for the rotor wake was determined using the equation, 

... ' — N 

c J 

^ 0 W c 

where the integration was performed in the tangential direction over one blade 
spacing. The momentum thickness implies the loss of momentum in the wake as 
compared with potential flow. The wake shape factor in the wake was found using 
the following equation. 
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where 5 is the displacement thickness. The displacement thickness is defined by 
the expression. 
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where the integration was also performed over one blade spacing in the tangential 


direction. 
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Figure 7.20. Variation of Rotor Wake Cycle-Averaged Shape Factor with 
Streamwise Distance 
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Figure 7.19 shows that the momentum thickness decreases in the trailing 
edge region, then increases as the wake travels downstream and then decreases in 
the far wake region. Both Raj and Lakshminarayana (1973) and Ravindranath 
(1979) exp lain ed the variation in momentum thickness in their compressor cascade 
and rotor wakes, respectively, on the basis of the von Karman momentum integral 
equation with zero wall shear stress as follows. 


— + (2 + H)— ^ = 0 
ds W ds 


(7.22) 


Equation (7.22) shows that the increase of decrease of 0 depends on the variation 
of W 0 . If W o increases then 0 decreases and if W 0 decreases then 0 increases. 

While this holds for both the compressor cascade and rotors, it does not hold for 
the turbine rotor wake presented here. The wake edge velocity continuously 
decreases downstream of this turbine rotor while the momentum thickness 
variation does not correspond to the wake edge velocity decrease. This is due to the 
radial transport of properties in the rotor wake which cause the variation in 
momentum thickness. 

The variation of shape factor with streamwise distance is presented in 
Figure 7.20. The shape factor decreases sharply in the trailing edge region and then 
decreases at a slower rate as the wake travels downstream. The high value of shape 
factor just downstream of the trailing edge (H = 1.95) shows that the flow has a 
tendency to separate there. Since a turbine blade has a thick trailing edge, the flow 
does have a tendency to separate there. Other researchers such as Hobson and 
Lakshminarayana (1990) have also predicted flow separation at the trailing edge of 
a turbine blade. 

The variation of the shape factor with streamwise distance downstream of 
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the trailing edge of an isolated airfoil was given by Spence (1953), 



where H te is the shape factor at the trailing edge. Raj and Lakshminarayana (1973) 
demonstrated that eqn. (7.23) can be used to accurately used to predict the variation 
of shape factor downstream of a compressor cascade. Figure 7.20 shows the 
comparison of this equation with the experimentally measured shape factor. In the 
trailing edge region the comparison is not good, but farther downstream the 
agreement between eqn. (7.23) and the experimental data is better. In the rotor far 
wake region, eqn. (7.23) accurately predicts the magnitude of shape factor. The 
poor agreement in the trailing edge region is due to the large three-dimensional 
nature of the flow in the trailing edge region, while the good agreement in the far 
wake region results from the reduction of the three-dimensional effects far 
downstream of the trailing edge. 



Figures 7.21a and b show the relative velocity profiles in the rotor wake at 
different axial locations close to the rotor trailing edge. The data is presented at two 
selected nozzle/rotor locations that have large difference in flow properties 
between them (locations 3 and 5). The wake width at position 3 is wider than that 
at position 5. This is because the nozzle wake is located on the suction side of the 
rotor wake at position 3 and is located in the freestream region outside of the rotor 
wake at position 5. While the magnitude of the minimum velocity in the wake 
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Figure 7.21. Rotor Wake Relative Velocity (w/U m ) Profiles 
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center is similar close to the trailing edge at both positions, the freestream velocity 
is quite different. The nozzle wake causes the freestream relative velocity to be 
lower at position 5 than the freestream velocity at position 3. 

The total unresolved unsteadiness in the rotor wake at nozzle/rotor positions 
3 and 5 are shown in Figure 7.22a and b. At both positions the unsteadiness profile 
is asymmetrical about the wake centerline which results from the different 
unresolved unsteadiness profiles in the blade surface boundary layers upstream of 
the trailing edge. The unsteadiness is much higher in the rotor wake and the region 
of increased unsteadiness occurs over a larger area at position 3 than at position 5 
due to the presence of the nozzle wake inside of the rotor wake at position 3. On 
the other hand, the free stream region of position 5 contains twice the magnitude of 
total unresolved unsteadiness as the free stream region of position 3. This also 
results from the nozzle wake which is located in the free stream region of position 
5. 

Figures 7.23a and b and 7.24a and b present the streamwise-normal 
unresolved velocity cross correlation and the relative flow angle, respectively, at 
nozzle/rotor positions 3 and 5. The magnitude of unresolved velocity cross 
correlation in the pressure side of the rotor wake at position 3 is twice that in the 
pressure side of the rotor wake at position 5. On the suction side of the rotor wake 
at position 5 the unresolved velocity cross correlation is negligible while the 
suction side of the wake at position 3 contains a large area of negative unresolved 
velocity cross correlation. The higher unresolved velocity cross correlation in the 
position 3 rotor wake results from the presence of the nozzle wake inside the rotor 
wake at this position. 

A noticeable difference between the relative flow angle in the rotor wake at 
these two nozzle/rotor positions is also seen in Figures 7.24a and b. At position 3, 
the overturning in the suction side rotor wake is larger than the overturning in the 
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Figure 7.23. Rotor Wake Unresolved Velocity Correlation (u'v' 
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Figure 7.24. Rotor Wake Relative Flow Angle (p) 
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suction side rotor wake of position 5, while the overturning in the freestream 
region outside of the rotor wake at position 5 is higher than that in the freestream 
region of position 3. The higher overturning regions result from the presence of the 
nozzle wake in these areas. 

Figure 7.21 through 7.24 show that rotor wake varies between individual 
nozzle/rotor positions. Thus one can conclude that the rotor wake is not steady in 
the rotor time frame. This can be seen more clearly by looking at Figures 7.25 and 
7.26 which show the rotor wake relative total velocity and total unresolved 
unsteadiness just downstream (X r /C r = 1.066) of the rotor trailing edge at all six 

nozzle/rotor positions. There is a variation in both the velocity inside the wake and 
in the freestream for all six positions (there is an 8% difference between the highest 
and lowest freestream velocity and an 18% difference in velocity inside the wake). 
The lowest velocity at the wake center occurs at position 3 due to the presence of 
the nozzle wake inside the rotor wake at this position. At position 4, the nozzle 
wake has moved toward the rotor wake freestream region. This causes the suction 
side wake width to increase and the freestream velocity to decrease, since the 
nozzle wake is located partly in the suction side of the rotor wake and partly in the 
freestream region. The nozzle wake has moved out of the rotor wake at position 5, 
since the wake center velocity has increased and into the freestream region where 
the velocity has decreased. 

The total unresolved unsteadiness, presented in Figure 7.26, shows the 
increased levels of unsteadiness in suction side of the rotor wake at position 3 due 
to the nozzle wake. At position 4, the level of the unsteadiness has decreased in the 
wake but has increased in the freestream region since the nozzle wake now is in 
both the rotor wake and freestream region. And at position 5, the nozzle wake has 
moved entirely into the freestream region as shown by the higher level of 
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7.25. Rotor Wake Relative Velocity (w/U m ) Profiles at X t /C, - 1.066 




Figure 7.26. Rotor Wake Total Unresolved Unsteadiness (Tu,)at X r /C r = 1.066 
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unsteadiness in the freestream region and lower level in the rotor wake. 

The decay of the rotor wake velocity defect with streamwise distance is 
shown in Figure 7.27 for each of the different nozzle/rotor positions. The decay 
rate for each of the positions is similar but the magnitude of the velocity defects are 
different at each streamwise location. In the trailing edge and near wake region the 
defects for positions 3 and 4 are larger than the defects at 5 and 6 due to the 
presence of the nozzle wake inside the rotor wake at positions 3 and 4. Thus while 
the nozzle wake does not seem to have an impact on the velocity defect decay rate 
it does influence the magnitude of the velocity defect. 

The nozzle wake interaction with the rotor wake also affects the rotor wake 
semi-wake width as shown in Figure 7.28. The semi-wake width is defined as the 
width of the wake at half the defect of total velocity. The rate of increase in wake 
width for each of the nozzle/rotor positions is similar but the magnitude of the 
semi-wake width at each streamwise location is not the same for each nozzle/rotor 
location. In the trailing edge and near wake regions, the rotor semi-wake width for 
positions 3 and 4 are larger in magnitude than the semi-wake width for positions 5 
and 6 due to the presence of the nozzle wake inside the rotor wake at positions 3 
and 4. Thus, just as for the velocity defect, the nozzle wake influences the 
magnitude of the semi-wake width but does not affect the semi-wake width rate of 
increase. 

Figures 7.29 and 7.30 present the decay of the maximum total unresolved 
unsteadiness and streamwise-normal shear stress, respectively. Just as with the 
velocity defect and the semi-wake width, the rate of decay is similar for all six 
nozzle/rotor positions; just the magnitude of the unsteadiness and shear stress is 
different for each streamwise location. Therefore, the nozzle wake interaction with 
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Figure 7.27. Decay of Rotor Wake Velocity Defect with Streamwise Distance for 
Each Nozzle/Rotor Location 



308 



Stream wise Distance, Z/cos/S 0 


Figure 7.28. Variation of Rotor Wake Semi-Wake Width with Streamwise Distance 
for Each Nozzle/Rotor Location 




Figure 7.29. Decay of Maximum Relative Total Unresolved Unsteadiness 

(Tu ) with Streamwise Distance for Each Nozzle/Rotor Location 
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Figure 7.30. Decay of Maximum Unresolved Velocity Correlation (u'v' n )with 
Streamwise Distance for Each Nozzle/Rotor Location 



311 


the rotor wake does not change the rate of decay or increase of the various 
properties, but it does change the magnitude of the properties for each nozzle/rotor 
position at each streamwise location. 
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CHAPTER 8 

SUMMARY AND CONCLUSIONS 

The three-dimensional steady flow field inside an axial flow turbine nozzle 
and the two-dimensional unsteady flow field inside a turbine rotor have been 
investigated experimentally. The nozzle flow field was studied in order to better 
understand the nozzle secondary flow and the nozzle wake properties. A complete 
flow field survey was carried out at two locations inside the nozzle passage at X/C 
= 0.56 and 0.935 and at two locations downstream of the nozzle at X/C = 1 .025 and 
1.09. The nozzle surface and endwall static pressures were also measured. The 
nozzle flow field measurements were carried out with a five hole probe, a single 

sensor hot wire and a two component LDV. 

The rotor flow field was measured at midspan with a two component LDV 
in order to better understand the steady and unsteady flow field in a turbine rotor. 
Measurements were acquired at 37 axial locations from just upstream of the rotor 
to one chord downstream of the rotor and at 50 tangential locations (relative to the 
rotor) over one rotor pitch. To account for the non-uniformity of the rotor inlet 
flow field (due to the nozzle wake), measurements were made at six tangential 
locations (relative to the nozzle) equally spaced over one nozzle pitch. The rotor 
wake was also studied in detail in order to understand the rotor wake properties and 
its decay characteristics. Some important conclusions that can be drawn based on 
the results of the present investigations are presented below. 
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8.1. Nozzle Flow Field 

The major conclusions of the nozzle flow field investigation are as follows: 

1. The nozzle surface static pressures were predicted better by the three 
dimensional inviscid code of Katsanis than by the panel code, thus demonstrating 
that some of the departure between the design and the measured data is caused by 
the three-dimensional inviscid effects. The agreement of the yaw angle and the 
total velocity prediction by Katsanis' code and the experimental data is excellent, 
demonstrating the usefulness of using an inviscid code in early design. 

2. The endwall static pressures show that the minimum static pressure on 
both endwalls occurs at the hub wall/suction surface comer. The minimum static 
pressure location on the hub occurs downstream (X/C = 0.80) of the minimum 
static pressure region on the casing, which occurs at midspan. This low pressure 
region is the location where the passage vortex lifts off the endwall and begins to 
grow rapidly as it moves up the suction surface. Thus the casing passage vortex 
covers a larger area than the hub passage vortex by the time it reaches the nozzle 
trailing edge, due to the longer distance it has traveled between the endwall/suction 
surface intersection point and the tr ailin g edge. 

3. Near midchord, the LDV measurements indicate the existence of the 
casing passage vortex by the high turbulence intensities and the undertuming of the 
flow in the casing wall/suction surface comer. On the other hand, at the hub there 
is no indication of the hub wall passage vortex. This is in agreement with the 
endwall static pressure distribution, discussed above. While the casing passage 
vortex meets the suction surface at midchord, where it begins to grow rapidly and 
starts to move up the suction surface (away from the endwall), the hub passage 
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vortex is still being convected across the passage at midchord and is small. Thus 
the hub passage vortex should be rather weak and hard to detect at midchord. 

4. Just upstream of the trailing edge (X/C = 0.935), the five hole probe data 
clearly show that the casing passage vortex is very strong. Radially inward flow in 
the suction surface boundary layer augments the casing passage vortex. The suction 
side leg of the horseshoe vortex is not visible near the casing. On the other hand, 
the passage vortex in the hub/suction surface comer is very weak. The radially 
inward flow in the suction surface boundary layer is in the opposite direction as 
those induced by the passage vortex and thus counteracts it. A vortex rotating in 
the opposite dir ection as the hub passage vortex is visible above the passage vortex. 
This might be the suction side leg of the horseshoe vortex. 

5. Downstream of the nozzle, weak radial inward flow was observed over 
the whole span, which was more pronounced in the wake. Casing and hub passage 
vortices have been identified on the suction side of the wake. There is also 
indication of other vortices located in the wake center, near the casing and the hub. 
They rotate opposite to the passage vortices and are caused by the interaction of the 
passage vortices and the wake. On the other hand, there is conflicting evidence on 
the existence of the suction side leg of the horseshoe vortex. While there is no 
evidence for its existence at the hub, there is a region of negative vorticity near the 
casing at X/C = 1.025 that could correspond to this vortex, but the secondary flow 
vectors do not show the existence of a vortex. It has probably decayed by the time 
it has reached this axial location (X/C = 1.025). 

6. Comparing the nozzle flow field just upstream of the trailing edge to that 
downstream of the nozzle, the casing passage vortex is seen to remain strong as it 
progresses and no dramatic change in secondary flow occurs at the tip. At the hub, 
however, the radially inward flow of the suction surface boundary layer has 
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reversed direction and is moving outward as a result of the rotating hub located 
downstream of the nozzle. The rotating hub also causes the hub passage vortex to 
dissipate, which can be seen from the secondary velocity vectors at X/C = 1.09. 

7. A nozzle wake survey was also done at several axial locations at midspan 
to determine the wake decay characteristics. The AFTRF nozzle wake decays much 
faster than a compressor cascade wake, an annular turbine nozzle cascade wake or 
a turbine nozzle wake with a large rotor-nozzle spacing. This is due to the presence 
of a rotor in close vicinity, as well as the influence of a favorable pressure gradient 
downstream. 

8. The radial variation of nozzle wake total velocity defect shows larger 
velocity defects in the hub and casing secondary flow regions. These larger defects 
result from the interaction of the passage vortices and the wake which cause deeper 
and wider wakes. 

8.2 Rotor Passage Steady and Unsteady Flow Field 

The major conclusions of the rotor passage flow field investigation are as 
follows: 

1. Detailed measurements were made near the rotor leading edge. These 
measurements showed that the rotor leading edge has a major influence on the flow 
field, with large velocity gradients and flow angle changes in the vicinity of the 
leading edge. The axial velocity in the stagnation region decelerates to a value of 
30% of the ffeestream axial velocity. The rotor leading edge has an influence on 
the flow field even 9% of the rotor axial chord upstream of the rotor leading edge, 
with the change in flow angle between the freestream and stagnation region being 
18 degrees at this location. This effect increases as the rotor leading edge is 
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approached with the change in flow angle increasing to 35 degrees just upstream of 
the rotor leading edge. 

2. High levels of relative total unresolved unsteadiness are observed at the 
leading edge. This is caused not only by a reduction in relative velocity, but also by 
an increase in the absolute magnitude of the velocity fluctuations at the leading 
edge. As the leading edge is approached, this large increase in the mean flow 
velocity gradients cause an increase in the production of unresolved unsteadiness 
which overshadows its dissipation. 

3. Higher levels of relative total unresolved unsteadiness are also observed 
near the rotor pressure surface. This increase in unresolved unsteadiness is due to 
the interaction of the nozzle wake with the pressure surface boundary layer, along 
with the concave curvature effects, which destabilize the flow. Concave curvature 
effects results from the fact that the unresolved unsteadiness in thin shear layers is 
hi ghl y sensitive to streamline curvature in the plane of the mean shear. The 
unresolved shear stress and unsteadiness are increased due to curvature when the 
angular momentum of the flow decreases in the direction of the radius of curvature. 

4. The nozzle wake can be clearly identified in the rotor passage. Compared 
to the flow field outside of the wake, the nozzle wake has a velocity defect, a 
higher unresolved unsteadiness, overturning of the flow angle and a higher 
unresolved velocity cross correlation. 

5. As the nozzle wake enters the rotor passage, it is chopped into individual 
segments by the rotor blades. These individual segments can now move 
independently of each other. At the rotor inlet, the nozzle wake becomes bowed 
because the convection velocity at midpitch is higher than at the rotor leading edge. 
As the nozzle wakes travels further through the rotor passage, it becomes distorted 
with the region of the nozzle wake near the rotor suction surface moving faster than 
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the region near the pressure surface. This is due to the large differential in the 
convection velocity between the pressure and suction surface. Inside the rotor 
passage, the nozzle wake acts like a negative jet, drawing fluid from the rotor 
pressure surface to rotor suction surface, which causes the nozzle wake to thin out 
near the rotor pressure surface and thicken near the suction surface. This movement 
of low momentum fluid toward the suction surface also causes the high momentum 
fluid in the free stream to move in an opposite direction to replace the migrated 
wake fluid. 

6. The nozzle wake propagates through the turbine rotor until it has turned 
more than 30 degrees from its orientation at the rotor leading edge and is spread 
out along the rotor pressure surface from midchord to the trailing edge. This is in 
contrast to other measurements of the nozzle wake in the rotor passage (Hodson, 
1984 and Binder et al., 1985) who show that at the rotor trailing edge, the nozzle 
wake still spans the rotor passage in the pitchwise direction, from the pressure 
surface to the suction surface. This is caused by the large difference in velocity 
between the rotor suction and pressure surfaces. In the present turbine, the ratio of 
suction surface velocity to pressure surface velocity is twice the suction to pressure 
surface velocity ratio of Hodson's and Binder's turbine rotors. This large difference 
between the velocities near the suction and pressure surfaces of the author's 
turbine, causes the region of the nozzle wake near the rotor suction surface to travel 
much more rapidly than the region of the nozzle wake near the rotor pressure 
surface. Thus the nozzle wake rotates about the region of the nozzle wake near the 
rotor pressure surface, ending up spread out along the pressure surface. In Hodson’s 
and Binder's turbine rotors, the difference between the suction and pressure surface 
velocities is not as great, thus the nozzle wake does not turn as much in the rotor 
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passage and by the time it reaches the rotor trailing edge it still spans the rotor 
passage in the pitchwise direction, from the pressure to suction surfaces. 

7. By the time the nozzle wake has propagated through the rotor passage 
and is spread out along the rotor pressure surface from midchord to the trailing 
edge, it has elongated and thinned out considerably (compared to its shape at the 
rotor leading edge). This results from two reasons. The first is the fact that vorticity 
must be conserved in the nozzle wake, thus as the wake length grows, the wake 
width must become smaller. The second is due to the large variation in convective 
velocity across the rotor pitch. Since the velocity is faster along the rotor suction 
surface than near the pressure surface, the region of the wake near the suction 
surface moves faster than the region near the pressure surface. Thus the wake 
stretches, becoming narrower and longer. 

8. As the nozzle wake interacts with the rotor leading edge, an increase in 
unresolved unsteadiness at the leading edge is observed. The tangential component 
of the unresolved unsteadiness also increases in the region of the nozzle wake at 
rotor midpitch. This increase in unresolved unsteadiness results from the two 
counter-rotating vortices, one on each side of the nozzle wake. The generation of 
these vortices cause an increase in the velocity fluctuations. 

9. The periodic velocity correlations are similar in magnitude to the 
unresolved velocity correlations in the nozzle wake, which indicates that both have 
an equal impact on the through-flow mixing in this turbine rotor. 

10. As the nozzle wake moves through the rotor passage the magnitude of 
the velocity defect increases from the leading edge until X r /C r = 0.30, after which 
it decreases. This can be explained by examining the mass-averaged cycle- 
averaged relative velocity. The relative velocity decreases from the rotor leading 
edge until X r /C r = 0.30, causing an adverse pressure gradient, which causes the 
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nozzle wake velocity defect to increase. Downstream of X r /C, = 0.30, the relative 
velocity increases, resulting in a favorable pressure gradient, which causes the 
nozzle wake to decay. Downstream of the rotor trailing edge the nozzle wake 
velocity defect is constant. 

1 1 . The peak magnitude of the total unresolved unsteadiness in the nozzle 
wake increases initially and then decreases as the nozzle wake travels through the 
rotor passage, until at one half chord downstream of the rotor trailing edge it is at 
the same value as it was upstream of the rotor. 

12. Both the mass-averaged and time-averaged properties show that inside 
the rotor blade, the region of high unresolved unsteadiness that corresponds to the 
nozzle wake occurs at the same location as the low velocities due to the nozzle 
wake. On the other hand, downstream of the rotor, the region of high unresolved 
unsteadiness due to the nozzle wake does not occur at the same location as the low 
velocities due to the nozzle wake. This is because, although in the near wake the 
maximum unresolved unsteadiness occurs at the wake center (position where the 
velocity defect is largest), far downstream of the nozzle, the maximum unresolved 
unsteadiness occurs away from the center of the wake due to the spread of the 
wake. 

83. Rotor Wake Properties and Characteristics 

The major conclusions based on the rotor wake investigation are as follows: 

1 . The rotor wake mean velocity profiles in the trailing edge and near wake 
regions are asymmetrical and tend to become symmetrical in the far wake region. 

2. The rotor wake total unresolved unsteadiness is highest in the trailing 
edge region and decreases further downstream. This is caused by the vortex street 
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shed from the trailing edge and the production of turbulence in this region. Also in 
the tr ailin g edge region, the unresolved unsteadiness profiles are asymmetric with a 
dip in unsteadiness at the wake center. The asymmetry results from the differential 
boundary layer growth on the pressure and suction side of the rotor blade, while the 
dip in unste adin ess is due to the fact that the total unresolved unsteadiness is zero 
on the blade surface with a maximum value occurring slightly away from the blade 
surface. The asymmetry and dip in unresolved unsteadiness disappear as the rotor 
wake travels downstream. 

3 The streamwise-normal unresolved velocity cross correlation in the rotor 
wake is also asymmetric close the rotor trailing edge, becoming symmetric farther 
downstream. Close to the trailing edge, the zero unresolved velocity cross 
correlation point does not occur at the point of minimum velocity , but farther 
downstream it does. The unresolved velocity cross correlation decays rather slowly 
with significant values still occurring one half chord downstream of the rotor. 

4. The rotor wake velocity defect decays very rapidly in the trailing edge 
region, becoming less rapid in the near and far wake regions. The very rapid decay 
in the trailing edge region results from the high unresolved unsteadiness, separation 
and possibly the three dimensional effects in this region. The rotor wake velocity 
defect decays much slower than the velocity defect of turbine cascades. The decay 
of the rotor wake velocity defect decays according to 1 / -Jz / cos|3 0 where z / cos (3 0 

is the streamwise distance downstream of the rotor blade. Using this property of the 
rotor wake, a correlation has been presented which matches the decay of the 
velocity defect in the near and far wake regions. 

5. The rotor wake semi- wake width increases rapidly in the trailing edge 
region and then grows less rapidly in the near and far wake regions. The rate of 
semi- wake width increase is proportional to i/z / cos p 0 . Using this property of the 
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rotor wake, a correlation has been presented which matches the growth of the semi- 
wake width. 

6. Decay of the maximum unresolved unsteadiness and maximum 
unresolved velocity cross correlation is very rapid in the trailing edge region and 
this trend slows in the far wake region. Equations have been presented which 
match their decay rates. 

7. In the trailing edge region, the maximum periodic velocity correlations 
are much larger than the maximum unresolved velocity correlations. But the 
periodic velocity correlations decay much faster than the unresolved correlations. 
Equations were presented that matched the experimental decay of the maximum 
periodic and unresolved velocity correlations in the rotor wake. 

8. The rotor wake shape factor decreases sharply in the trailing edge region 
and then decreases at a slower rate farther downstream. The high value of shape 
factor just downstream of the rotor trailing edge indicates that the flow has a 
tendency to separate there. The variation of the rotor wake shape factor with 
streamwise distance downstream of the rotor was found to match the correlation 
developed by Spence (1953) quite well in the far wake region and not as good in 
the near wake region. 

9. While the interaction of the nozzle wake with the rotor wake does not 
influence the decay rate of the various wake properties, it does change the 
magnitude of the properties. 
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CHAPTER 9 

RECOMMENDATIONS FOR FUTURE RESEARCH 

The present investigation has provided a good understanding of the whole 
nozzle flow field. One unanswered question in the nozzle flow field is what 
happens to the suction side leg of the horseshoe vortex. Flow visualization on the 
nozzle endwall surfaces and hot wire or five hole probe surveys close to the 
endwalls could be used to trace the path of the suction side leg of the horseshoe 
vortex as it travels through the nozzle passage. 

The nozzle wake could be investigated in more detail. By measuring the 
properties of the nozzle wake with the rotor removed, the difference between the 
nozzle wake decay with and without the rotor could be quantified. 

This investigation has also provided an understanding of the rotor flow field 
at midspan and the interaction of the nozzle wake with the rotor flow field. To 
better understand the rotor flow field including the secondary flow and the tip 
leakage vortex, a knowledge of the flow field at other radii is critical. Measurement 
of the radial velocity component is also essential for a complete understanding of 
the flow field. The LDV can be used to measure the radial velocity and to measure 
the flow field at other spanwise locations. Measurement of the losses in and 
downstream of the rotor would also be useful. This could be accomplished with a 
rotating five hole probe. The AFTRF has a rotating traverse that the five hole probe 
could be installed in. With this traverse the losses in the entire flow field could be 
measured. 

There are numerous static pressure taps on the rotor blade surface and 
endwalls which should be used to measure the pressure distribution on the rotor 
blade. The unsteady pressures and the shear stress on the rotor blade surface can be 
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measured by the Kulite transducers and the shear stress gauges, respectively on the 
rotor blade surface. The rotor blade boundary layers could also be measured by 
using a two sensor hot wire mounted in the rotating traverse. By completing all 
these measurements, the complete flow field in the rotor could be measured. 

The influence of the inlet turbulence intensity on the flow field could be 
quantified by repeating selected measurements with the inlet turbulence grid 
installed and comparing them with the measurements completed without the 
turbulence grid. The influence of the Reynolds number and the incidence angle on 
the turbine flow field could also be quantified by measuring the flow field at 
several locations at off-design conditions. 
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APPENDIX A 
ERROR ANALYSIS 

All experimental results have a certain degree of error associated with them. 
This error needs to be quantified in order for the results to be correctly interpreted. 
Thus an error analysis based on the ASME measurement uncertainty methodology 
(Abernathy et al., 1985) is given below. 

An error is the difference between the measurement and the true value. 
Uncertainty is the maximum error that might be expected from a measured 
quantity. The total error is a combination of both bias (fixed) and precision 
(random) errors. The precision error, S, is the random part of the total error which 
for N measurements (X,,X 2 ,...,X N ) of the parameter X , is 


X(X,-X) 2 


s = 


-1 J=L 


N- 1 


(A.l) 


where X , the average value (mean) of X is 


i N 

x=-yx, 

N" 


(A.2) 


The precision index of the sample mean X can be also found from 


S v = 


Vn 


(A. 3) 
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The bias error is a fixed or systematic error. It remains constant during an 
experiment. In repeated measurements, each measurement has the same bias. There 
is no statistical equation, as for the precision error, to define the bias error, B. It 
must be estimated. This estimate can be based on calibrations or comparison with 
other experimental techniques, but in general the estimate of bias must be based on 
judgment. 

To obtain the precision error of a given parameter, the root sum square 
(RSS) method is used to combine the precision errors from K sources of error as 
follows, 


S = [Sf+S 2 +...+S 2 ]* (A.4) 

The bias of a given parameter can be similarly found as 

B = [B? + B 2 +...+B 2 ]* (A.5) 

The total uncertainty is obtained for a 95% confidence level) by combining the 
precision an bias errors in the following manner 

U = [B 2 + P 2 ]' (A. 6) 

where P = tS* . For a large number of samples (>30), t=2.0. U is then the 
uncertainty in the measured value of a single variable. If, on the other hand, the 
value of the experimental result is not directly measured, but the values of several 
variable are measured and then they are combined into a data reduction equation to 
obtain the value of the desired variable, the following is used. 
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Consider a general case in which an experimental result, r, is a function of J 
variables X i 

r = r(X 1 ,X 2 ,...,X J ) (A.7) 


The precision limit of the result is given by 


S = 


dr 

v ax, 


+ 


— s x 

dX, 



(A. 8) 


where the S are the precision limits in the measured variables X,. It is assumed 

that the relationship given by equation (A.7) is continuous and has continuous 
derivatives in the domain of interest and that the precision limits S x are 

independent of one another. 

This same method can be used to calculate the experimental bias error from 
the bias limits of each measured variable that the experimental result is a function 
of (shown in equation A.7) as is done as follows 
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The uncertainty of the result is again given by (for 95% confidence level) 


U = [B r 2 + P r 2 ]- 


(A. 10) 


where P r = tS r and t=2.0 when N>30. 
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If it is difficult to separately estimate the bias errors and the precision errors 
for each measured variable, an uncertainty which is a combination of the bias and 
precision errors can be estimated for each measured variable as follows. 


U = [U 2 +U 2 +...+U 2 ] 


(A. 11) 


where U ,U ,...,U are the sources of error for the measured variable X j5 and U 
is the uncertainty in X, . Then the total uncertainty of the experimental result, r, 
which is a function of J variables X; as shown in equation (A. 7), is then 
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(A. 12) 
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APPENDIX B 

FTVE-HOLE PROBE AND STATIC PRESSURE UNCERTAINTY 

ANALYSIS 


B.l Sources of Error 

A detailed discussion on the errors involved in the five hole probe and static 
pressure measurements in given by Sitaram, Lakshminarayana and Ravindranath 
(1981) and by Zierke, Straka and Taylor (1993). The present discussion is based on 
the above mentioned references using the ASME measurement uncertainty analysis 
described in Appendix A to calculate the uncertainty. 

The accuracy of the five hole probe measurements is dependent on several 
factors including transducer error, wall and blade vicinity effects, Reynolds number 
and compressibility effects, misalignment of problem probe blockage effects and 
turbulence effects. Each of these factors are discussed in detail below: 

(1) Transducer uncertainty: The error in the transducer is 0.15 inches of 
water. This error affects both the five hole probe and static pressure measurements. 
The remaining errors discussed below only affect the five hole probe 
measurements. 

(2) Wall Vicinity Effect: When a probe is located near a solid surface, the 
flow accelerates in that region introducing errors in the measurements. Sitaram, 
Lakshminarayana and Ravindranath (1981) found that for a five hole probe, the 
error is negligible when the distance between the probe and the solid surface is 
more than two probe diameters. Thus this error is negligible since all 
measurements were made at a distance greater than two probe diameters from the 


wall. 
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(3) Reynolds number effects: The calibration of the five hole probe was 
carried out at a Reynolds number approximately equal to that in the measured flow. 
Thus the Reynolds number effects are accounted for in the calibration. 

(4) Compressibility effects: Although the five hole probe was calibrated at 
the same Reynolds number as that in the measured flow, the calibration was done 
in an incompressible flow (water). To account for compressibility effects of the 
flow, the pitch and yaw calibration constants were modified using the Prandtl- 
Glauert rule as follows. 


C 


Ppudxc 



(B.l) 


C 


Pyiw.c 



(B.2) 


where the subscripts c and i represent compressible and incompressible calibration 
coefficients. The total velocity is then calculated using. 


V 2 


2c 2 



Y-l 


VP 



(B.3) 


Thus the error due to compressibility can be neglected since it is accounted for. 

(5) Misalignment of the probe: The probe was aligned using a protractor. 
The error in aligning the probe is 0.25 degrees. 

(6) Probe and stem blockage effects: The error in probe and stem blockage 
is negligible since the probe was calibrated with its stem and the ratio of the probe 
diameter to both the vane spacing and pitch is 1.4% and 1.5%, respectively. 
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(7) Turbulence Effects: The five hole probe was calibrated in a low 
turbulence calibration tunnel. Measurements were taken in the nozzle wake where 


the turbulence intensity is around 10%. Goldstein (1936) has theoretically 
investigated the effect of turbulent velocity components on pressure probe 
measurements. Sitaram, Lakshminarayana, and Ravindranath (1981) have modified 
Goldstein's analysis to estimate the error due to turbulence on the five hole probe 
measurements. Based on their analysis, the error in the calibration angle and 
pressure coefficients is 0.67%, assuming 10% turbulence intensity in the nozzle 
wake. 

(8) Pressure and Velocity Gradient Error: The five hole probe used in this 
investigation was calibrated in a uniform flow, but was used to measure flow where 
steep gradients in pressure and velocity exist (such as in the nozzle wake). Since 
each hole is located in a differing pressure field, an error is introduced. This is 
corrected by using a linear interpolation of the pressures at two adjacent measuring 
locations in the direction of the gradient, as devised by Prato (1992). Thus, the 
pressure and velocity gradient error can be neglected. 

B.2 Uncertainty Analysis 

The error in static pressure results mainly from the error in the pressure 
transducer which is 0. 15 inches of water. Thus for the static pressure measurements 

^ = 1.0% (B.4) 

P 

For the five hole probe measurements, the data reduction equations are based on 
the method of Treaster and Yocum (1979) which is modified to include 



compressibility and shear gradient effects. The equations for this analysis are 
presented and then the propagation of uncertainty from the measured variable is 


given. 


The first step in the five hole probe data analysis is to calculate an average 


pressure. 


P = ^(P2 + P 3 +P4+P5) 


(B.5) 


where p is the average pressure, p 2 and p 3 are the pressures measured in the yaw 
plane, and p 4 and p 5 are the pressure measured in the pitch plane. The uncertainty 


equation is then 


U. =-U„ 
p 2 Pj 


(B. 6 ) 


Next, the yaw and pitch calibration coefficients are calculated as 

f V 2 ZI 1 
r' _ V Pi _ P > 




1 -M 


(B.7) 


P 4 Ps 
c - VP'~P 
Pp “ 1 i-m! 


The general uncertainty equation for these coefficients is then 


U r = 


P2-P3J U2-P3J \pi-p) Ipi-pJ u-m 2 ; 


(B. 8 ) 


(B.10) 


u„ V ( -U V -U V ( U- V f-MU Y 
U c = — ^ + M + p —\ + ■ H 

IP4-P5J IP4-P5J vp, - py vPi-py U-M 2 ; 


(B .11) 


After the uncertainty in C p ^ and C are computed, then the uncertainty in 
the angles can be calculated. Since the values for the pitch angle, y, and yaw angle, 
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a, come from the calibration plots, the partial derivatives in the general uncertainty 


equation are replaced with finite differences to obtain 


u =! 
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(B.12) 
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(B. 13) 


The calibration pressure coefficients, which are used to compute the pressures from 
five hole probe calibration plots, are calculated next. The uncertainty in these 
coefficients are 


and 
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The total and static pressures are then calculated using the following equations 

p 0=p,-c;jp,-p) (B.16) 

p = p -c L.(pi-p) (BI7 ) 


Thus the general uncertainty equations for the total and static pressures are 
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u ^(( 1 - O u J Ju J 4 p - p, )U C J 2 ? (B.18) 


u p =f((- c LJ u J 2+ (( i+c ;„«) u p) J+ ((p-p.) u c-_) a T (B.19) 


The total velocity can then be calculated as follows. 


r-i Ip 


(B.20) 


and the uncertainty equation for V 2 is 


U, = V 2 


U T V 2RT 


t ) + [ y2( a 


^ U P 


(B-21) 


Finding the magnitude of the velocity vector from the squared value leads to 




(B.22) 


Next, the axial, radial and tangential components of velocity in the probe 


coordinate system are computed as 


V = Vcosy cosa 


V r = Vsiny 


(B.24) 


V e = Vcosy sin a 


(B.25) 
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The uncertainty equations for these components are then 

U Vx = [(u v cosy cosa) 2 +(-U r V sin ycosa) + (-U T VcosYsina) j (B.26) 

U Vf| = [(U v siny) 2 + (U 7 V cosy) j (B.27) 

1 

U v ^ =[(u v cosysina)' + (-U 7 Vsinysina) + (-U a V cosy cosa) j (B.28) 

The final step in the five hole probe data analysis is to convert the velocities 
from the five hole probe coordinate system to the turbine coordinate system. This 
is done as follows, 

V x = [V x cosa- V 6| sina]cos8- V r sin 5 (B.29) 

V r = [V X| cosa- V e sinaJsinS- V r cosS (B.30) 

V e = V 9i cosa - V Xi sin 8 (B .3 1 ) 


where a and 8 are the yaw and pitch angle offset, respectively, of the five hole 
probe relative to the turbine coordinate system. The uncertainty equations for these 
components are then calculated as 
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Using equations (B.12), (B.13), (B.18), (B.19), (B.22), (B.32), (B.33). and (B.34), 
the uncertainties for the properties measured by the five hole probe are calculated 
to be 


U a = 0.5 degrees 
U y = 0.5 degrees 
U„ =30 Pa 

r o 

U p = 33 Pa 


U v = 1 m/s 
U v =0.78 m/s 
U v = 0.45 m/s 

v e 

U v =0.78 m/s 

v r 

Using typical values downstream of the nozzle, the percentage uncertainty for the 
velocities are: 

U 

= 0.6 % 


u 


V, 

u. 


- = 2.5 % 

- = 0.5 % 


u 


- = 22 % 
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APPENDIX C 


HOT WIRE UNCERTAINTY ANALYS 


The accuracy of hot wire measurements can easily be quantified and are 
limited only by the accuracy of the calibration, the anemometer and voltmeter error 
and by ambient temperature drift. Each of these factors are discussed below. 


C.l Probe Calibration 

The single sensor hot wire was calibrated in a low turbulence (< 1.0 %) 

calibration jet. The sensor was placed perpendicular to the flow. The probe was 

calibrated over a range from 0 to 0.6U m . A pitot-static probe connected to a 

manometer was used to measure the calibration jet velocity. The air velocity 

measured by the pitot-static probe is given by 

2RTp HjQ gh 

P 



(C.l) 


where h is the manometer reading in inches of water. The manometer is used to 
measure the difference between the static and total pressures for the pitot-static 
probe. The error in reading the manometer is 

U hi = 0.025 inches H 2 0 (C.2) 

The uncertainty for the thermometer and barometer that are given by the 
manufacturer are respectively. 


U T) = 0. 1 degrees K 


(C.3) 
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U P| = 0.05 inches Hg (C.4) 

The significant source of error results from the pitot-static tube installation. 
According to Holman (1994) and Bradshaw (1970), the major source of error in a 
low turbulence flow away from solid surfaces is due to probe misalignment. The 
error due to the pitot-static installation is then 

U h2 = O.Olh = 0.02 inches H 2 0 (C.5) 

since h=2.0 inches of water is the h corresponding to the ma ximum velocity 
measured by the hot wire. 

The uncertainties in each of the measured variables is 



(C.6) 

u t =K+uU ! 

(C.7) 

u„=(u‘ +u’) ! 

(C.8) 


and the uncertainty in the calibration velocity measured by the pitot-static probe is 



ifUr 

H T 





(C.9) 
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Substituting the nominal values for T, h and P and equations (C.8), (C.9), and 
(C.10) into (C. 11) gives 

^ = 0.8% (C.10) 


C.2 Measurement Errors 

The velocity measured by the hot wire sensor is obtained from King's law 
which is 


E 2 = A + BV 045 


(C.ll) 


where E is the anemometer output voltage taken across the Wheatstone bridge in 
the anemometer, V is the fluid velocity and A and B are the calibration constants. 
The uncertainty in the measured velocity is thus 


U v = 


dv 

9E 


U t 


+ 


&- 


Ms"' 


(C.12) 


To find U v , the uncertainties in E, A and B need to be found first. The error stated 
by the manufacturer for the anemometer bridge output voltage and the voltmeter 
output voltage are 

U Ei = 0.0002V (C.13) 

Up = 0.002V (C.14) 

fc 2 


and the error due to reading the voltmeter is 
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U Ej = 0.0005V (C.15) 

Another source of error that can change the anemometer output voltage is ambient 
temperature drift. A change in ambient temperature causes a change in the heat 
transfer rate of the sensors. This results in a corresponding change in the measured 
output voltages from the anemometer. A correction scheme by Kristensen (1973) is 
used to correct the variation in anemometer output voltage due to temperature. 

Thus this error is neglected. The total uncertainty in E is then found from 

(C.16) 

The calibration constants A and B are obtained by fitting a straight line in 
the form 


InA = Ain V + B 


(C-17) 


through the calibration data. The uncertainty in A and B are then obtained from the 
equations 


U A = 


(IN' 


av v 


(C.18) 


u B = 
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aE Ue J + [av Uv 


(C.19) 
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and are assumed to be zero.) Thus, substituting the uncertainties in 

A, B and E into equation (C. 12) along with the nominal values for A, B and E 
yields 


U 


f =l9% 


(C.20) 


which is the total uncertainty in the measured velocity. Yavuzkurt (1984) states 
that the mean and rms components of velocity have the same percent error, thus 


U 


Tu 


Tu 


= 1.9 % 


(C.21) 


which is the total uncertainty in the measured turbulence intensity. 
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APPENDIX D 
LDV ERRORS 


LDV measurements are subject to numerous errors, most of which can be 
qualified. The discussion here, which is based on the error analysis of Patrick 
(1987) uses the ASME measurement uncertainty analysis described in Appendix A 
to calculate the uncertainty. The errors are separated into both precision and bias 
errors. Bias errors include errors from laser beam geometry, counter processor 
errors and seeding bias errors. The laser beam geometrical errors consist of finite 
probe volume bias, beam location bias, beam orientation bias, fringe spacing 
uncertainty bias, negative velocity beam bias, angle bias and frequency broadening 
bias. The processor bias errors are made up of errors due to comparison accuracy, 
clock synchronization, quantizing, threshold limit, electronic noise and pedestal 
filter removal. Finally, the seeding bias errors include errors resulting from the 
flow distortion, particle lag, statistical or velocity bias and Bragg bias. Most of the 
bias errors are very small compared to the precision errors (discussed below), and 
thus are neglected. The bias errors which can be on the order of magnitude of the 
precision errors are discussed in detail below. 

Statistical or velocity biasing was first mentioned by McLaughlin and 
Tiederman (1972). It occurs as a result of two reasons. The first is that the velocity 
magnitude varies with time. The second is that in a uniformly seeded flow, more 
particles pass through the probe volume per unit time during periods when the 
velocity is faster than the mean velocity than when the velocity is slower than the 
mean velocity. Thus a high data rate causes the measured values to be biased 
toward a higher value than the true mean. In the nozzle, where the turbulence 
intensity was less than 5 %, the error due to statistical bias is less than 0.5% based 
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on the analysis of Strazisar and Powell (1980). In the rotor, based on the same 
analysis, this error is about 1% outside the rotor wake and 5% inside the rotor wake 
(where the highest turbulence intensities occur). But in the rotor the need for 
statistical biasing corrections is questionable. According to Hathaway (1993) and 
Bell et al. (1993), measurements taken using conditional sampling techniques (such 
as with an optical shaft encoder) should not need this correction, since the velocity 
data is grouped into individual time windows and averaged for each window, and 
just because there are more measurements at higher velocities means there are 
more measurements at that particular window, but the average velocity is unbiased. 
Based on this, the velocity bias error is neglected. 

Angle bias occurs when the flow is not parallel to the plane containing the 
laser beams. The factor controlling the angle bias is the ratio N/Nfj- , where N is the 
minimum number of cycles required by the signal processor and Nfj- is the number 
of measurable fringes. The angle bias can be minimized by reducing the N/Nfj- 
ratio. Frequency shifting was used to minimize the angle bias, and thus this error is 
negligible also. 

The precision errors in LDV measurement are data processing errors which 
result from averaging a finite number of data samples per data point. In LDV 
measurements, the velocity being measured does not remain constant during the 
sampling period but fluctuates due to turbulence. Thus the precision error is 


where 


N 2 
S(v,-v) 

i=l 


S = \^ = v' 

(C.l) 

V N - 1 

i N 

v=-yv, 
n£ 1 

(C-2) 
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where V s is the velocity of the i’th sample and V is the sample mean velocity. The 
precision error calculated by equation (C.l) is an estimate of the rms turbulence 
level. For an infinite number of samples, S becomes an exact measurement of the 
turbulence level. For a finite number of samples, both the rms turbulence level, v', 
and the mean sample velocity V, will deviate from the true turbulence level and 
mean velocity of the flow field by precision errors S v . and S v , respectively. 

Patterson (1982) states that the mean square turbulence intensity has a Chi-square 
distribution. Thus for a large sample size (N>50) the precision error of the 
turbulence intensity measurement can be found from 


Sy- 1 

v # V2N 


(C.3) 


Relative to the mean velocity measurements, the sampling distribution of V is 
normal about the population mean V p (true mean) as a mean with a standard 

deviation of v p /Vn . Thus the precision error in the mean velocity measurement 
can be estimated as 



which is a function of the true turbulence intensity,(v'/V p ). The quantity, (v'/V p ), 

is unknown but can be approximated by the measured ratio of the turbulence to the 
mean velocity, (v'/v). 

For the LDV measurements in the rotor, 120,000 axial and tangential 
coincident samples were taken at each measurement position for the entire rotor 
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revolution. Since the data were phase-locked averaged to obtain an "average" rotor 
blade passage with 50 measurement windows, there were on average 2400 samples 
per bin. Thus outside of the rotor wake, the precision errors in the mean velocity 
and turbulence intensity are less than 0.2% and 1.4% respectively. In the rotor 
wake, where the samples per bin are much smaller and the turbulence intensity is 
much higher, the precision errors in the mean velocity and turbulence intensity are 
2% and 7%, respectively. For the LDV measurements in the nozzle, where the 
sample size was 200, the precision errors in the mean velocity and turbulence 
intensity are 0.1% and 5.0%, respectively. Thus using equation (A.6), the total 
uncertainty for a 95% confidence level is as follows. 

Nozzle flow field: = 0.2 % 

S v//, = 10 % 


Rotor flow field: 

Outside of rotor wakes: 


s y ^ = 0.4 % 
s = 2.8 % 
s y ^ = 4.o % 
s v^/, = 14.8 % 


Inside rotor wakes: 
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